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Conventions 


The  symbol  |  |  means  the  magnitude  of  the  terms  contained  within 
The  symbol  Re{  }  means  the  "real  part"  £.  Im(  }  means  "imaginary  part" 

The  symbol  < - >  indicates  a  Fcurier  transform  pair 

The  symbol  oc  refers  to  "proportional  to" 

The  use  of  *  over  a  symbol  means  the  estimate  of  an  unknown  parameter 
Random  variables  are  uppercase;  random  variable  values  are  lowercase 
A  /  within  parantheses  {/A}  means  "given  that"  ("given  that"  A  is  true) 
The  symbol  ln[  J  denotes  the  natural  logarithm  of  the  enclosed  quantity 
The  symbol  sqrt(  )  means  th*»  square  root  operation 
The  symbol  E(  )  means  the  expected  value 
The  symbol  var{  }  means  the  variance 
The  symbol  cov(  }  means  the  covariance 

The  symbol  p{  }  means  the  probability  of  the  enclosed  random  quantity 
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ABSTRACT 

This  investigation  addresses  the  problem  of  estimating  the  target 
angle  off  beam-axis  of  an  amplitude-comparison  monopulse  radar  to 
generate  an  error  signal  that  can  be  applied  to  a  servo-cdntrol  system  . 
to  reposition  the  beam-axis  on  target,  thus  providing  target  tracking. 
Cramer-Rao  (CR)  and  Ziv-Zakai  (ZZ)  bounds  are  derived  to  indicate  system 
performance  under  varying  signal-to-noise  ratio  (SNR)  conditions. 

Actual  tracking  error  is  approximated  from  a  computer  simulated  tracking 
loop,  and  then  compared  to  the  CR  and  ZZ  bounds  for  varying  SNR  levels, 
for  tracking  in  "slow"  and  "rapid"  target  fluctuation  environments,  and 
for  tracking  in  the  presence  of  specular  and  diffuse  multipath.  A  two- 
dimensional  tracking  model  and  associated  ZZ  performance  bound  are  also 
presented. 

At  high  SNR  levels,  the.CR  bound  results  lower  bounded  the  mean- 
square  tracking  error,  but  for  low  SNR  conditions  the  CR  bound  exceeded 
the  mean-square  error.  The  ZZ  bound  results  indicate  a  tight  lower 
bound  for  the  mean-square  tracking  error  at  low  SNR  levels  and  in  both 
the  target  fluctuation  and  multipath  environments.  "Slow"  target 
fluctuations  and  diffuse  multipath  results  indicate  that  target  tracking 
capability  is  not  seriously  degraded  in  either  of  these  two 
surroundings.  Conversely,  "Rapid"  target  fluctuations  and  specular 
multipath  environment  results  indicate  serious  tracking  degradation  is 
introduced  in  the  amplitude-compariron  monopulse  tracker. 
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ZIV-ZAXAI  POUND  APPLIED  TO  AN 


AMPLITUDE-COMPARISON  MONOPULSE  RADAR 

I.  INTRODUCTION 

Background 

An  amplitude-comparison  monopulse  radar  is  a  tracking  radar  that 
determines  target  direction  by  comparing  the  reflected  signal  amplitudes 
received  simultaneously  on  two  identical  but  noncoincident  antenna 
patterns.  Target  tracking  is  accomplished  by  maintaining  the  axis  of 
symmetry  of  the  antenna  patterns  (beam-axis)  on  the  target.  As  the 
target  moves  off  beam-axis,  the  amplitude-comparison  yields  an  error 
voltage  proportional  to  the  angular  difference  between  the  present  beam- 
axis  location  and  the  target  position.  This  error  information  is 
applied  to  a  servo-control  system  to  reposition  the  beam-axis  on  target, 
thus  providing  target  tracking. 

Due  to  practical  considerations  such  as  channel  noise  and  servo 
time-constant  limitations,  the  radar  beam-axis  will  be  steered  in  the 
vicinity  of  the  actual  target  position..  Associated  with  this  vicinity 
is  the  tracking  error  from  the  true  target  location.  An  accurate  bound 
on  the  tracking  error  provides  a  direct  indication  of  the  quality  of  the 
estimr ted  target  position  and  thus  provides  a  measure  of  the  ability  of 
the  system  to  track  a  target. 

Traditionaly,  the  tracking  error  has  been  lower  bounded  by  the 
application  of  the  Cramer-Rao  bound  to  the  variance  of  the  distribution 


of  the  estimated  angle  off  beam-axis.  This  bound  limits  the  variance  of 
the  estimate  to  a  value  directly  proportional  to  the  radar's  beamwidth 
(BW)  and  inversely  proportional  to  the  available  signal -to-noise  ratio 
(SNR) . 

var( estimated  position)  ae  BW 

SNR 

Inspection  of  the  Cramer-Rao  bound  indicates  that  for  any  signal- 
to-noise  ratio,  (SNR),  reducing  the  beamwidth  decreases  the  variance  of 
the  estimate.  This  presents  a  shortcoming  of  the  Cramer-Rao  bound  as 
this  relationship  will  not  hold  for  all  BW  and  SNR  values.  Thus,  SNR 
levels  for  which  the  Cramer-Rao  bound  holds  (ie  the  variance  of  the 
estimated  position  is  directly  proportional  to  the  beamwidth)  must  be 
defined  if  the  bound  is  to  be  useful. 

Problem  and  Scope 

This  thesis  concentrates  upon  generating  a  maximum-likelihood 
estimate  of  the  target  angle  off  beam-axis,  determining  the  mean-square 
error  between  the  estimate  and  the  actual  target  position,  and  applying 
the  Ziv-Zakai  bound  to  verify  system  performance. 

The  maximum-likelihood  estimate  of  an  unknown  parameter  is  the 
parameter  value  that  would  most  likely  cause  a  given  observation  to 
occur,  in  terms  of  statistical  information,  the  maximum-likelihood 
estimate  would  be  the  value  of.  an  unknown  parameter  for  which  the- 
conditional  probability  distribution  function,  given  the  unknown 
parameter,  is  maximum,  Thit*  thesis  develops  a  maximum-likelihood 
estimate  of  the  target  angle  off  beam-axis  and  uses  the  estimate  to 


track  the  target 


I 


The  optimum  solution  of  the  estimation  of  an  unknown  quantity  from 
data  is,  in  general,  very  difficult  analytically.  The  solution, 
however,  can  be  considerably  simplified  if  the  permitted  analytical 
operations  are  optimum  in  the  mean-square  sense.  The  error  of  the 
target  position  in  the  mean-square  sense  is  the  average  value  of  the 

squared  difference  between  the  estimated  and  actual  target  position: 

2 

mean-square  error  =  E{( estimated  position  -  actual  position)  ) 

This  thesis  considers  the  difference  between  the  actual  target  position 
and  the  maximum-likelihood  estimated  position  in  the  mean-square  sense. 

The  Ziv-Zakai  bound  is  designed  to  be  tighter  to  actual  performance 
values  at  low  signal-to-noise  ratio  levels.  The  bounds  established  by 
the  Ziv-Zakai  method  are  derived  by  comparing  -he  estimation  problem 
with  the  known  results  from  detection  theory.  The  derived  bound 
compares  a  suboptimal  estimate  of  one  of  two  possible  values  of  a 
parameter  to  an  optimal  detection  scheme  between  the  two  possible 
values.  A  criteria  for  determination  of  which  unique  parameter  value 
was  estimated  then  leads  to  a  simple  comparison  between  the  probability 
of  etror  associated  with  optical  and  suboptimal  detection  methods.  A 
form  of  the  Ziv-Zakai  bound  is  used  in  this  thesis  to  derive  a  lower 
bound  for  the  mean-square  tracking  error. 

The  goal  of  this  thesis  is  to  solve  the  general  problem  of 
estimating  an  amplitude-comparison  monopulse  radar  target  parameter  and 
to  apply  the  Ziv-Zakai  bound  to  develop  an  accuracy  estimate  for  system 
performance  as  a  function  of  signal-to-no<se  ratio. 
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The  development  and  solution  of  this  thesis  depends  upon  the 
following  assumptions: 

1 )  The  channel  noise  is  assumed  to  be  additive  white  Gaussian 
that  is  independent  from  channel  to  channel  and  pulse  to 
pulse  with  zero  mean  and  variance  equal  to  HQ/2. 

2)  The  model  is  restricted  to  one  dimensional  tracking;  thus 
total  independence  from  a  second  dimension  is  assumed. 

3)  Any  velocity  changes  are  assumed  to  be  slowly  varying  with 
respect  to  the  time  between  radar  pulses. 

4)  The  antenna  patterns  are  derived  assuming  a  uniform 
aperture  illumination  and  the  target  is  in  the  "far- 
field" . 

5)  The  amplitude  of  the  return  signal  is  assumed  either 
unknown,  but  non-randcn,  (deterministic)  or  a  Rayleigh 
distributed  random  variable  with  completely  specified 
characteristics . 

6)  The, carrier  phase,  although  unknown,  is  assumed  to  be  the 
same  in  all  receiver  beams. 
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2)  Process  the  observations  to  quantify  a  maximum-likelihood 
angle  error  estimate  that  will  be  used  to  steer  the  radar 
beam-axis  to  the  predicted  target  location. 

3)  Derive  Cramer-Rao  and  Ziv-Zakai  bounds  to  indicate  system 
performance  under  varying  signal-to-noise  ratio 
conditions . 

4)  Design  and  implement  a  computer  simulated  tracking  locp 
that  provides  for  beam-axis  position  changes  and  will 
yield  the  error  between  actual  target  position  and  the 
beam-axis  position. 

5)  Convert  the  tracking  error  into  a  mean-square  error  and 
plot  for  different  signal-to-noise  ratio  levels. 

6)  Compare  the  mean-square  error  plots  to  the  derived  Cramer- 
Rao  and  Ziv-Zakai  bounds  to  indicate  which  bound  provides 
the  best  estimate  of  accuracy,  in  the  mean-square  sense, 
for  different  signal-to-noise  ratio  levels. 

7)  Observe  radar  performance  under  situations  of  slow  and 
rapid  fading  of  the  signal  amplitude,  and  in  the  presence 
of  a  return  embedded  in  multipath. 
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II.  RADAR  MODEL 


The  Antenna  Functions 

Amplitude-comparison  monopulse  radar  uses  two  overlapping  antenna 
patterns  to  obtain  ang’e  tracking  for  a  single  coordinate  direction. 

The  patterns  are  usually  mirror  images  about  the  boresight  axis  (beam- 
axis).  When  a  target  is  on  the  boresight  axis,  the  signals  received 
from  the  two  patterns  are  equal.  Target  tracking  can  be  accomplished  by 
bringing  the  boresight  axis  into  coincidence  with  the  target  bearing,  as 
indicated  by  equal  received  signals. 


Figure  1.  Simple  Block  Diagram  of  Antenna  Circu:  — * 

Processing  of  the  received  signals  is  accomplished  by  applying  the 
antenna  patterns,  A  and  B,  to  a  hybrid  junctio..  which  suss  and  subtracts 
the  two. patterns.  A  block  diagram  of  the  amplitude»coaiparison  antenna 
circuitry  is  shown  in  Figure  1. 
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The  sum  of  the  two  beam  patterns,  A  +  B,  is  represented  as  an 

antenna  function  V  .  Likewise,  the  difference,  A  -  B,  is  represented  as 
s 

V,.  To  derive  the  two  antenna  functions,  V  and  V_,  the  following  steps 
d  s  a 

will  be  taken: 

1 )  The  technique  for  taking  the  Fourier  transform  of  a  uniform 
source  will  be  developed. 

2)  A  general  expression  will  be  derived  for  the  magnitude  of  the 
electric  field  at  a  distance  far  from  the  source  (far-field) . 

3)  An  expression  will  be  developed  for  squinting  (offsetting)  the 
antenna  patterns  from  the  boresight  axis. 

4)  Two  antenna  patterns  will  be  combined  with  one  squinted  in  a 
positive  direction,  and  the  other  squinted  an  equal  amount  in 
the  negative  direction.  These  two  patterns  will  then  be  added 
and  subtracted  to  form  and  V'd>  respectively. 

5)  Approximations  will  be  applied  to  the  general  expressions  Of  4 
to  develop  simple  expressions  for  the  antenna  functions. 

Antenna  Patterns 

Each  individual  antenna  pattern,  A  or  B,  is  assumed  to  be  generated 
from  a  uniformly  illuminated  rectangular  aperture.  A  uniform 
rectangular  aperture  has  aperture  fields  which  are  uniform  in  phase  and 
amplitude  across  the  physical  aperture.  The  described  situation  for  . 
uniform  illumination  in  the  x~y  plane  is  shown  in  Figure  2. 
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Figure  2.  Uniformly  Illuminated  Rectangular  Aperture 
in  the  x-y  Plane 


To  specify  antenna  patterns,  Fourier  transform  techniques  are  used. 
In  general,  the  Fourier  transform  is  used  to  convert  from  time  domain  to 
frequency  domain.  As  applied  to  antennas,  the  Fourier  transform  can  be 
used  to  convert  from  direct  space  (aperture  illumination:  to  K-space 
(antenna  far-field  pattern)  (1:523].  Three  Fourier  transform  tools  will 
be  utilized  in  the  far-field  pattern  development.  These  tools  are: 

1)  f(t)  < - >  F(u) 

2)  f'(t)  < - > juF(u) 

3)  f(t)  ■  f(t+a)  <—-->?( u)  •  exp(jua) 

where: 

f(  )  refers,  to  the  illumination  of  the  source 
F(  )  refers  to  the  f*r-field  antenna  pattern 

designates  the  Fourier  transform  operation 
primed  variables  refer  to  derivatives  of  that  variable 
}(  )  refers  to  the  Dirac  delta  function  { 2 : 273) 
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Figure  3.  Graphical  Determination  of  f'(t) 


Consider  f(t)  to  be  a  uniform  source  of  length  2a,  the  derivative, 
f'(t),  can  be  determined  graphically  as  shown  in  Figure  3,  Using  the 
Fourier  transform  tools  previously  mentioned,  the  derivative  of  the 


uniform  source  can  be  utilized  to  determine  the  K-apace  representation 
of  the  far-field  pattern. 

f'(t)  <— ~«> juF(u)«(2a)  1 [exp( jua)-exp(-jua)] 

Solving  for  F(u)  yield3: 

F(u)  *  (2aju)  1 )exp( jua)-exp(-jua)) 

*  (2aju)  1 (2jsin(ua)] 

*  (ua)  ^sin(ua)  *  Sa(ua)  (1) 

where  t 

Sa(  )  *  sampling  function  13:24) 


•j. 


Figure  4.  The  Rectangular  Aperture 


To  apply  the  Fourier  expression  for  the  far-field  pattern,  suppose 

the  aperture  electric  field,  E  ,  is  known  and  x  directed  as  shown  in  the 

& 

1  by  1  rectangular  aperture  of  Figure  4.  For  the  one-dimensional 
x  y 

problem,  the  source  geometry  takes  the  form  of  a  uniform  line  source 
along  the  x-axis.  The  pattern  of  concern  is  the  electric  field  in  the 
x-z  plane  that  will  be  generated  by  this  uniform  source.  For  the  x 
directed  source,  the  variable  u  will  describe  the  x  component,  of  the 
phase  constant  for  a  plane  wave.  This  phase  constant  is  commonly 
refere4  to  as  k  and  equals  2 (pi) /wave length.  Making  a  simple 
rectangular  to  spherical  coordinates  conversion,  k^  « 

(sin(d)cos(f )] 2 (pi) /wavelength.  Substituting  this  value  into  equation 
1,  and  noting  that  the  distance  a  <*  1^/2,  yields  an  expression  for  the 
magnitude  of  the  electric  field  due  to  a  uniform  source. 

E  *  Sa(2(pi)/wavelength) [sin(9)cos(» ) ) 1^/2 


(1A) 


For  the  one-dimensional  radar  model,  further  simplification  of 
equation  1A  can  be  made  by  noting  that  only  the  magnitude  of  the 
electric  field  in  the  x-z  plane  is  of  interest  (ie  electric  field 
changes  as  a  function  of  azimuth  only).  In  the  x-z  plane: 

»  -  0 

cos(D)  =  1 

Simplifying  equation  1A: 

E  *  Sa[( (pi)lx/wavelength)sin(8)]  (2) 

For  large  apertures  (ie  lx  >>  wavelength)  the  sin(8)  factor  is 
approximately  equal  to  8  [1:387].  Further,  the  factor  ( pi ) Wwavelength 
can  be  expressev.  as  [4:267]: 

( pi ) 1^/wavelength  «  2.78/BW  (3) 

where: 

BW  *  3dB  beamwidth 

Substituting  equation  3  into  equation  2,  a  final  form  for. the  normalized 
magnitude  of  the  electric  field  for  the  uniform  rectangular  aperture  can 
be  expressed  as: 

E  *  Sa(2.788/BW) 

A  plot  of  equation  4  for  varying  8/BW  ratios  is  provided  in  Figure  5. 


(4) 


Figure  5.  Antenna  Pattern  for  Uniformly  Illuminated  Line  Source 


Squinting  Beams  Off  Boresight 

Amplitude-comparison  monopulse  uses  two  overlapping  patterns,  of 
the  nature  of  equation  4.  that  are  squinted  (offset)  from  the  boresight 
(tracking)  axis.  The  boresight  axis  is  an  imaginary  line  drawn 
perpendicular,  and  centered  on,  the  plane  of  the  rectangular  apertures. 
This  squint  relationship  is  shown  in  Figure  6,  where  one  beam  is 
squinted  in  a  positive  direction  from  boresight  by  the  angle  9^,  and, the 
other  beam  is  squinted  negatively  from  boresight  by  9#. 

The  squinting  of  the  two  beams  from  boresight  will  cause  two 
distinct  antenna  voltage  gains  to  be  received  from  a  target  that  is 
arbitrarily  placed  at  some  angle  9  . 
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Figure  6 -  Squinting  of  Beams 


Figure  7  depicts  the  geometry  of  the  described  situation.  Beam  number 

one  would  locate  the  target  at  an  angle  »t  -  Similarly,  beam  number 

two  locates  the  target  at  an  angle  of  •  ♦  «  . 

t  s 
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In  order  to  locate  a  target  position  relative  to  the  boresight 
axis,  the  monopulse  radar  sums  and  subtracts  the  voltage  gain  from  each 
of  the  two  squinted  beams.  The  antenna  functions  Vg  (sum  of  the 
patterns)  and  (difference  of  the  two  patterns)  depend  upon  completely 
specifying  the  squinted  beam  patterns.  For  a  target  located  at  an  angle 
0 : 


V,(0)  »  £(0  -  0  ) 

1  s 

(5) 

V  ( 0 )  -  E(0  +  0  ) 

2  S 

(6) 

where : 

V  ( .8 )  *  received  voltage  due  to  positively  squinted  beam 
and  is  a  function  of  the  target  angle  0 
V2(0)  '*  received  voltage  due  to  negatively  squinted  'beam 
and  is  a  function  of  the  target  angle  0 

Using  equations  5  and  6  with  the  patterns  described  by  equation  4: 

V  (0)  -  Sa((2.78/BW)(0  -  0  ))  (5A) 

•  s 

V„(0)  «  Sa((2.78/BW)(0  +  0  ))  (6A) 

2  3 

At  this  point  of  the  antenna  function  development,  the  squint  angle 
is  established  as  one-half  the  be&mwidth.  With  this  assumption, 
equations  5a  and  6A  becomes 

Vn(0)  «  Sa{2 .780/BW  -  1.39)  (SB) 


V2<0)  *  Sa( 2 .780/BW  +  1.39) 


(6B) 


Figure  8.  Voltage  Sain  of  Two  Beams 

Squinted  Off  Bores ight  by  BW/2 


A  combined  plot  of  equations  SB  and  6B  are  shown  in  Figure  8.  The 


antenna  functions  then  follow  as  V  (0)  »  V_(3)  ♦  V_(9)  and  V.(0)  *  V,(0) 

3  i  2  a  l 

-  V2(0) .  The  development  for  the  sum  pattern  is  as  follows: 


V  (9) 
s 


v^e)  +  v2(9) 

Sa(2.788/BW  -  t.39)  +  Sa(2.780/BW  ♦  1.39) 

[2.780/BW 


-  Sa(2.780/BW  -  1.39) 

♦  Sa( 2.789/BW  +  1.39) 


[2.788/BW 
! 2.788/BW 


+  1 .39]" 
+  1.391> 


[2.780/BW 


-  1.39]1 
•  1 . 39] j 


let  a  »  1.39s  b  «  2.780/BW  - 

( 9 )  «  (a2  -  b2)  Vlsinta  -  b)  *  sin<a  ♦  b))a  VI 

[[sin(a  •  b)  +  sin(a  ♦  b)]b  J 

*  (a2  -  b2)  ' f I-co8(a)sin(b)  *  cos<a)sin<b)]b  +1 

[ lsin(a)cos(b)  ♦  sin(a)cos(b)]a  J 

•  (a2  *  b2)  1 l2asin(a)coa(b)  *  2bcos(a)sin(b) ] 
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2 


[  .5<8/BW)sin(2. ?80/BW)-1 .37cos(2.738/BW)]  (7) 


7.73 [0/BW] 2  -  1.93 

A  similar  development  for  the  difference  pattern  yields: 

Vd(0)  _ 2 _ 1 .25sin( 2 . 780/BW) -2 . 74( 0/BW)cos(2.780/BW) ]  (8) 

7. 73 [0/BW] 2  -  1.93 


In  order  to  make  subsequent  calculations  more  tractable,  and  to 

remove  the  function  discontinuities  sho«n  in  Figure  9,  equations  7  «nd  8 

need  to  be  simplified.  Observation  of  Figure  9  indicates  that  V  is  an 

■  8 

even  function  and  is  an  odd  function.  Consider  approximating  the  sum 
pattern  as  a  scaled  cosine  function  and  approximating  the  difference 
pattern  as  a  scaled  sine  function: 

Vg(approx.)  *  sqrt(2)cosn{«)  k  7A) 

V^tapprox.)  *  sinr{d)  (8A) 


Estimation  Model 


The  general  framework  of  the  estimation  problem  consists  of 
determination  (estimation)  of  the  value  of  an  unknown  parameter, 
possibly  random,  from  an  observation.  Normally  observations  are  made 
through  a  noisy  channel,  and  in  this  case,  estimation  techniques  must  be 
used  to  determine  the  desired  parameter  value.  Parameter  estimation  is 
then  a  mapping  of  the  observation  space  into  a  parameter  space. 


Figure  11  shows  the  basic  components  of  the  estimation  model.  Once 
a  suitable  rule  for  mapping  into  a  parameter  space  has  been  developed, 
it  is  of  interest  to  examine  measures  of  quality  of  the  estimation 
procedures,  possible  quality  measures  might  include  the  mean  and  the 
error  (in  the  mean  square  sense)  of  the  estimate.  Frequently,  quality 
measures  of  the  estimation  procedure  are  quite  difficult,  and  for  many 
cases  rather  than  approach  the  error  of  the  estimate  directly  it  is  of 
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more  concern  to  derive  a  lower  bound  on  the  variance  of  any  estimate. 

A  comparison  of  the  chosen  estimation  procedure  to  the  derived  lower 
bound  will  provide  the  needed  measure  of  quality  of  the  estimation 
procedure.  To  complete  the  estimation  model  for  the  amplitude- 
comparison  monopulse  radar,  the  following  steps  will  be  taken: 

1 )  The  observations  that  the  estimator  will  have  to  work  with  are 
derived  from  the  amplitude  weighting  provided  by  the  antenna 
functions  in  the  presence  of  additive  white  Gaussian  noise. 

2)  A  maximum-likelihood  estimation  procedure  will  be  applied  to 
the  received  observations  to  estimate  the  amplitude  that  most 
likely  caused  the  observation.  Subsequently,  the  amplitude 
estimate  will  be  used  to  develop  an  error  equation  based  upon 
the  angle  of  arrival  of  the  received  signal  from  a  target. 

3)  The  Cramer-Rao  lower  bound  will  be  developed  which  considers 
the  conditional  variance  of  any  unbiased  estimate  for  the 
target  angle  of  arrival. 

4)  An  alternate  bound,  the  Ziv-Zakai  hound,  will  be  developed  for 
the  specific  case  where  the  procedure  of  estimation  is  maximum- 
likelihood. 
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Observation  Formulation 


Figure  12.  Monopulse  Observation  Model 

> 


For  the  specific  case  of  additive  white  noise  (See  Figure  12)  the 


received  waveforms  are: 

r  (t)  -  As ( t ) V  +  n; (t) 

*  SI 

r  (t)  *  As( t )V  +  n^(t) 
2  <3  2 


(11) 

(12) 


r(t)  -  fr  (t)l 

[r2(t)J 

where; 

A  *  unknown  signal  amplitude  (nonrandom) 
s(t)  ■  transmited  signal 
V  »  sum  voltage  antenna  gain  *»  V  (0) 

8  S 

Vd  ■  difference  voltage  antenna  gain  ■  V^',8) 
n(t)  -  additive  white  Gaussian  noise 

In  general,  the  transmitted  signal  <s(t));  Me  disturbing  noise 
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(n(t));  and  the  received  signal  are  randon  processes.  The  key  to 
analyzing  the  random  observation  is  to  find  a  way  to  replace  all 
waveforms  by  finite  dimensional  vectors,  for  which  characterization  of 
the  random  process,  by  way  of  a  joint  density  function,  can  be  made. 

One  implementation  of  this  transformation  is  referred  to  as  Gram-Schmidt 
orthogonalization  [6:266].  The  orthogona-ization  procedure  permits  the 
representation  of  any  K  finite  energy  time  functions  as  linear 
combinations  of  M  £  K  orthonormal  basis  functions. 

M 

s,(t)  *£(s.  )»  (t)  k=1,2, .  ..K  (13) 

K .  K,m  m 

where : 

integration  is  taken  over  a  period 
{01(t),02(t)< ...^M(t))  are  orthonormal 


Figure  13.  General  Orthogonalization  Implementation 
Source:  [6:228] 
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Figure  13  illustrates  the  implementation  of  the  orthogonalization 

process.  Each  individual  observation  is  developed  by  the  basis 

function,  0  (t),  and  t  12  total  dimension  of  the  vector  observation 
m 

corresponds,  at  most,  ‘;o  the  number  of  signals  at  the  source.  For  the 
estimation  problem  at  hand,  the  nature  of  the  source  is  known  and  one- 
dimensional  (singular  valued).  The  random  variable  R^  will  be  the  only 
output  and  is  a  sufficient  representation  for  the  energy  time  function, 
r^(t),  as  all  other  R's  are  independent  of  R^  (with  white  Gaussian  noise 
assumed).  Further,  the  basis  function,  0^(t)  is  [6:267]: 

^(t)  =  s(  t)/sqrt(E)  (14) 

where: 

E  *  /s2(t)dt  *  signal  energy 
sqrt(  )  “  square  root  operation 
integration . is  taken  over  a  period 


RANDOM 

.  RANDOM 

PROCESS 

VARIABLE 

(t) - <X)— » 

T 

j()dt 

0 

— *1 

T 

s(U 

.  t 

Figure  14.  Monopulse  Orthogonalization 


*  A  «V. 


A-V 
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Substitution  of  equations  11  and  14  into  equation  13  yields: 

R1  =  Jr1(t)»1(t)dt: 

*  HAV  s(t)  +  n  ( t )  ]  p  ( t ) dt 

J  s  11 

*  fl AV  s2(t)/sqrt(E)]dt  +  N. 

J  3  i 

-  AV  sqrt(E)  +  N.  (15) 

s  1 

where: 

N1  -  J s(t)n1(t)/sqrt(E)dt 
integrations  are  taken  over  a  period 

A  similar  development  for  R^  yields: 

R2  »  AVdsqrt(E)  ♦  N2  (16) 

To  form  the  joint  density  function  required  to  characterize  the 

received  waveforms,  the  first  and  second  moments  of  tls  vectors  R^  and 

R2  must  be  determined.  Por  the  R1  vector,  and  conditioning  on  the 

amplitude  and  angle  of  arrival  (conditioning  is  performed  on  the  angle 

of  arrival  rather  than  the  received  sum  channel  voltage  because  the  sum 

voltage  antenna  pattern  is  a  function  of  9,  therefore  E(R,/A,v  )  is 

i  s 

proportional  to  E(Rj/A,9}): 

EtR^A.ej  -  AVssqrt(E)  ♦  E(  /s(t)n1  (t)/sqrt(E)dt) 

-  AVasqrt(£)  ♦  ^E(s(t) }E(n^(t)}/sqrt(E)dt 
•  »  AV  sqrt(E)  (17) 

S" 

where: 

signal  and  noise  are  independent, 
mean  of  the  noise  is  assumed  zero 
E(  )  refers  to  the  expectation  operation 
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where:  (continued) 


integrations  are  taken  over  a  period 

A  similar  development  for  the  vector  yields: 

E{R2/A,8)  *  AVdsqrt(E) 
where : 

E(R2/A,Vd)  is  proportional  to  E{R2/A,8) 

varCR^A.8)  -  E([R1  -  AVssqrt(E)]2/A,8) 

-  EHR^/A.S)  -  lAVgsqrt(E)  ]  2 

-  E{ ( AV  sqrt(E)  +N  )2/A,8)  -  [AV  sqrt(E)]2 

S  1  s 

Expanding  the  first  term  of  equation  19: 

E{ [AV  sqrt(E) ] 2/A, 8J  +  E{2AV  sqrt(E)N./A,8)  ♦  E{N  2/A,8) 
s  si  1 

-  (AV^sqrttE) ]  2  +  E{{  (u)s(t)s(u)/E)dtdu/A,8} 

Substituting  equation  19A  back  into  equation  19: 

var(R^/A, 8}  -  1/E  (t)^ (u) )E(s( t)s(u) Jdtdu 

-  V“//*“  -u)s(t)s(u)dtdu 
■  Nq/2E  Ja(n)alu)A\x 

-  NqE/2E  -  Nq/2 

where: 

j" i (tTU)dt  •  1»  iff  t  -  u  (sifting  property) 

■  Oi  otherwise 
integrations  are  taken  over  a  period 
each  noise  sample  is  independent,  Gaussian,  zero  mean 
autocorrelation  of  the  noise  <■  <N  /2)4(r) 
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(18) 


(19) 


(19A) 


(20) 


(21) 


A  similar  development  for  yields: 
var(R2/A,0}  *  NQ/2 

COV(R  , R  /A, 0 )  =  E{(R.R,)2)/A.0}  -  A2EV  V.  (22) 

12  i  2  s  d 

Expanding  the  first  term  of  equation  22: 

E{ (AVgsqrt(E)  +  (AVdsqrt(E)  +  N2J/A,0) 

-  A2V  ♦  AV  sqrt(E)E{N_/A, 0}  +  AV^sqrttEjECH./A, 9}  +  E{N,N,/A,0) 

sd  S  2  d  2  12 

-  A2v  V..E  (22A) 

s  d 

Substituting  equation  22A  for  the  first  term  of  equation  22: 

cov(R1#R2/A,0}  «  0  (23) 

Equation  23  shows  the  vector  observations  to  be  uncorrelated. 
Further,  with  Gaussian  noise  assumed  and  linear  operations  to  develop 
the  vector  observations,  R^  and  r2  will  be  Gaussian  and  uncorrelated, 
therefore  independent  [7:199].  The  receiver  observation  is  then  the 
joint  density  of  two  independent  Gaussian  random  variables. 

Independence  implies  that  the  joint  density  of  R^  and  R2  is  the  product 
of  the  marginal  R^  and  R2  distributions.  With  Gaussian  distributions 
and  statistics  given  by  equations  17,  18,  20,  21,  and  23  the  joint 
conditional  density  of  the  receiver  observation  is: 
pR(r/A,0)  -  p{ry/A.9)p(r2/\,9) 

*  [piN0]*1exp(-1/N0 [(r1  -  AVgsqrt(E))2  ♦  (r2  -  AVjSqrt(E) )2] )  (24) 
where: 

r  - 
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Maximum-Likelihood  (ml)  estimation 

The  maximum-likelihood  approach  is  to  choose  the  estimate  for  the 
unknown  parameter  to  be  the  parameter  value  that  most  likely  caused  a 
given  observation  to  occur.  In  the  general  case,  the  joint  conditional 
density  of  the  receiver  observation  is  denoted  as  the  likelihood 
function  £5 : 65]  .  Another  useful  function  for  monotonic  functions  (1  to 
1  mapping  between  the  function  and  the  likelihood  function)  is  the 
logarithm  of  the  likelihood  function.  With  a  1  to  1  mapping,  working 
with  the  log-likelihood  function  is  equivalent  to  working  with  the 
likelihood  function  [8:181].  The  maximum- likelihood  (ml)  estimate  is 
then  that  value  of  the  unknown  parameter  for  which  the  likelihood 
function  is  a  maximum.  Or  equivalently,  a  necessary  condition  for  the 
ml  estimate  is  obtained  by  differentiating  the  log-likelihood  function 
with  respect  to  the  unknown  parameter  and  setting  the  result  equal  to 
zero  [5:65] . 

Calling  this  condition  the  log-likelihood  equation,  it  can  be 
expressed  as : 

A* [ln[pR(r/A,e)]J |  -  0  (25) 

A  -  A  . 
ml 

where : 

A*  [  J  represents  the  partial  wrt  the  parameter  A 

Maximum- likelihood  estimates  can  then  be  obtained  from  unique 
solutions  of  equation  25.  To  apply  the  procedure  described  by  equation 
25  to  the  conditional  density  of  equation  24,  the  log-likelihood 
equation  must  first  be  established.  Using  equation  24 i 
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In [pR(  r/A, 9 )  ]  =- 1 /NQ l ( r 1 -AVgsqrt { E) ) 2+ ( r2~AVdsqrt ( E ) ) 2 ] -In [ ( pi  )NQ  ] 

*  -1/N  [ ( r  ) 2  -  2Ar  V  sqrt(E)  +  (AV  sqrt(E))2] 
Oils  s 

-1/N0l(r2)2  -  2Ar2Vdsqrt(E)  +  (AVdsqrt(E) )2J 
-ln[{pi)Ng] 


Takin-,  the  partial  wrt  the  parameter  A  and  setting  the  result  equal  to 

zero  will  provide  the  log-likelihood  equation  described  by  equation  25: 

0  -  2r  V  sqrt(E)/N.  -  2A(V  sqrt(E))2/H.  + 

IS  OS  0 

2r2Vdsqrt(E)/N0  -  2A(Vdsqrt(E) )2/NQ  (25A) 

*  r.V  sqrt(E)  -  A(V  )2E  +  r.,V.sqrt(E)  -  A(V12E 
is  s  2d  d 


The  other  unknown  parameter  of  interest  is  the  angle  of  arrival. 

An  estimate  of  9  allows  boresight  movement  to  keep  the  tracking  axis 
(boresight)  upon  the  target.  The  log-likelihood  function  terms  that  are 
functions  of  the  angle  9  are  V  and  V  which  are  described  by  equations 
9  and  10,  respectively.  Taking  the  partial  with  respect  to  the 
parameter  6  and  setting  the  result  equal  to  zero  yields: 

0  *  2AT  V  'sqrt(E)/N  -  2AV  sqrt(E)AV  'sqrt(E)/M  ♦. 

1  S  OS  S  v 

2Ar2Vd'sqrt(E)/N0  -  2AVdsqrt(E)AVd,sqrt(E)/M0  (27, 

-  r.V  *  -  AV  V  'sqrt(B)  ♦r.V.'  -AV.V  ’sqrt(E)  (27a) 

IS  SS  20  a  a 
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where 


V  V  '  represents  the  partial  wrt  the  parameter  0 

S  Q 


Substituting  equation  26  for  A  of  equation  27A  and  simplifying: 

0  «  r.V  ‘V*  -  r„V  V  .V  '+  r.V  'V  2  -  r.V  V  V  1 
Isd  2  s  d  s  2d  s  Isdd 


or  equivalently: 


r,v, [v  'V.  -vv;]  -  r.V  [V  'V  -  V  V  *] 
Id  sd  s  d  2  s  s  d  s  d 


equating  coefficients,  and  evaluating  at  the  position  estimate: 

A  A 

r.V.(8)  -  r  V  (0)  »  0  (28) 

1  a  2  S 

The  maximum-likelihood  estimation  procedure  has  generated  two 
estimates.  The  first  estimate,  equation  26,  describes  the  most  likely 
amplitude  for  a  given  observation.  The  second  estimate,  equation  28, 
was  derived  from  the  amplitude  estimate  and  gives  a  relationship  that 
contains  both  observation  terms  and  amplitude  weighted  angle  terms.  In 
order  to  get  a  better  feel  for  the  relationship  between  equation  28  and 
the  estimate  for  target  position,  it  is  necessary  to  look  more  closely 
at  equation  28. 

For  the  second  estimate  (equation  28)  the  random  quantities  will  be 
the  noise  corrupted  observations,  R^  and  R^.  In  a  no-noise  environment, 
the  observations  R^  and  R^  would  equal  their  respective  mean  values  and 
equation  28  restated  for  the  no^noise  case  is: 

^  *> 

0  -  [AV  (0)sqrt(E)]V.(0r  -  lAV.(0)sqrt(E))V  (0) 
s  a  a  s 


(29) 


•  ^ 


where : 

A  A  A 

V  (6),  V .  ( 8 )  are  functions  of  an  estimated  target  position,  8 

s  a 

V  (8),  V.(8)  are  functions  of  the  true  target  position,  8 

s  a 

Substitution  of  equations  9  and  10  into  equation  29  yields: 

A 

0  *  Asqrt( 2E)/2 [ ( 1  +  cos(2. 18/BW) )sin( 1 .888/BW)  - 

A 

(1  *  cos(2. 18/BW) )sin( * .888/BW) ]  (29A) 

Inspection  of  equation  29A  indicates  that  the  equality  will  be 
satisfied,  ir.  the  no-noise  case,  when  the  estimate  for  theta  equals  the 
parameter  value.  For  all  other  estimated  values  of  8,  unequal  to  the 
true  target  position,  the  equation  will  not  be  zero  but  will  be  an  error 
that  is  proportional  to  the  difference  of  the  estimate  from  the  actual 
value.  Further,  with  no-noise  assumed,  the  target  position  could  be 
determined  exactly  with  one  observed  pulse  (ie  the  required  estimate 
for  target  position  to  force  equation  29A  to  zero)  and  the  boresight 
could  ther.  move  to  the  exact  target  location.  With  additive  noise, 
equation  23  will  not  be  equal  to  zero  and  the  maximum-likelihood 
estimate  for  target  position  will  be  in  error.  In  this  light,  the 
second  estimate  (equation  28),  is  the  error  equation  for  the  amplitude- 
comparison  monopulse  radar.  The  derived  second  estimate  will  be  used  to 
wove  the  boresight  to  the  vicinity  of  the  target. 

In  order  to  bound  the  estimation  error  for  the  target  position,  it 
is  first  necessary  to  compute  the  statistics  for  the  generated  estimate 
for  the  amplitude  of  the  return  signal. 


(30) 
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Cramer-Rao  Lower  Bound  (CRB ) 


For  any  unbiased  estimate  A  of  a  scalar  A,  the  conditional  variance 
is  bounded  by  [8:232]: 

var{A/A)  _>  [E{(A*  [lntp^r/A))])2)]'1 

*  -1 

or  equivalently,  var{A/A)  >  [-E{A"[ln(p  (r/A))}}]  (32) 

*“  A 

where : 

A'  represents  the  partial  wrt  the  parameter  A 
A"  represents  the  second  partial  wrt  the  parameter  A 


Equation  32  is  usually  referred  to  as  the  Cramer-Rac  inequality. 
Consider  applying  the  steps  described  by  equation  32  to  the  amplitude 
estimate: 


A'[ln(pR(r/A))J 


2r,V  sqrt(E)/N  -  2AV  E/N.  + 
IS  0  so 


-A"[ln(pR(r/A))] 

E(-A"[ln(pR(r/A))]} 


-1 


2r2Vdsqrt(E)/H0  -  2AVfl  E/NQ 
2Vs22^0  +  2Vd2E/N0 


H. 


2  2 
2E(V  +  V  _) 
8  ct 


(33) 
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Using  the  inequality  given  by  equation  32: 
- 


vartA^/A,®}  * 


2E(V  2  +  V.2) 
S  Q 


(34) 


Equation  31  previously  defined  the  conditional  variance  of  the 
amplitude  estimate.  Comparison  of  equation  31  with  th*-  conditional 
inequality  of  equation  34  indicates  that  the  conditional  variance 
satisfies  the  CRB  with  equality.  Any  unbiased  estimate  that  satisfies 


the  inequality  of  equation  32  with  an  equality  is  called  an  efficient 
estimate  [5:66] .  There  is  no  guarantee  that  an  efficient  estimator 
exists  for  a  given  problem,  however,  if  one  doe  exist,  the  maximum- 
likelihood  estimate  will  be  efficient  [8:233].  Efficiency  is  of 
importance  because  equality  with  the  lower  bound  described  by  equation 
32  assures  minimization  of  the  conditional  variance  of  the  estimate,  and 
therefore,  provides  the  best  estimate  available.  Applying  equation  32 

/v 


to  an  estimate  of  target  location,  8: 

8 *  ( lnp  ( r/A, 0 ) ]  *  2Ar,V  ’sqrt(E)/N  -  2AV  sqrt(B)AV  *sqrt(E)/N  + 

R  IS  U  8  8  U 

2Ar2V/sqrt(E)/Nrt  -  2AVdsqrt{B)AVd’sqrt(E)/N0 

0"  [lnp_(r/A,8)]  -  2Ar  V  "sqrt(E)/Nft  -2A2(V  ’  )2E/N„  -  2A2V  V  "E/N  + 
R  IS  0  8  0  SSO 

2Ar2Vd"sqrt(E)/N0  -2A2  (V^  )2E/NQ  -  2A2VdVd"E/NQ 
*  2Asqrt(E)/N0[Va-(r1-AVasqrt(E))  +  V’(r2'AVd®qtt(E>  )J  - 
2A2E/N0[(Va')2  +  (Vd*)2] 

-  2Asqrt(E)/N(J[Va"N1  +Vd*N2  -Asqrt(E)  (V#*  )2  -Asqrt(E)  (V^  )2]  (35) 

-E(8"  [lnp  (r/A,  8 )  ]  }  «  2A2E/Nn[(V  ’)2  +  (V  ’)2]  (35A) 

R  us  a 


where: 

E(V  "N  )  «  E(V  }  »  0 
si  a  2 


C«E{8" [lnpR(r/A,9)J )] 


N. 


2A2E 


(v#*)2 


1 

*  'V>2 


(35B) 


Substitution  of  equations  9  and  10  into  equation  35B  and  performing 
the  indicated  first  partial  squared,  and  then  simplifying  yields  the 
final  form  for  the  Cramer-Sao  bound  (CBS)  for  the  conditional  variance 

/V 

of  the  unbiased  estimate  0t 


var(0/A,0)  > 


NQ [BW] 


2A2E 


[2.87  +  1 .8cos( 3 .80/BW)  -  1 . 1cos(4.20/BW)J 


(36) 


Equation  36  is  one  of  the  bounds  for  the  tracking  error  of  the 
monopulse  radar.  Application  of  the  conditional  variance  as  a  mean 
square  error  implies  a  priori  knowledge  about  the  average  error. 
Specifically,  a  zero  average  error  is  implied.  This  will  clearly  not  be 
the  case  if  the  estimate  is  no  longer  unbiased,  and  at  that  point  this 
form  of  the  CRB  loses  usefulness.  Before  leaving  the  CRB,  a  few 
comments  about  equation  36  are  in  order.  Inspection  of  the  leading 
coefficient  indicates  that  the  bound  is  directly  proportional  to  the 
beamwidth  and  inversely  proportional  to  the  signal-to-noise  ratio.  This 
result  is  intuitively  appealing,  as  the  bound  predicts  a  small  tracking 
error  in  a  low-noise  (high  SNR,  low  BW).  environment.  Nota  also  that  if 
the  converse  is  true  (low  SNR,  high  BW),  then  the  variance  of  the 
estimated  position  grows  quite  large.  For  low  SNR,  the  CRB  does  not 
adequately  approximate  the  estimation  error  because  the  likelihood 
function  can  have  several  peaks  and  maximization  of  the  likelihood 
function  may  only  produce  a  local,  instead  of  absolute,  maximum.  Errors 
are  then  made  on  the  likelihood  function  sidelobe  peaks  (9:148). 
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Ziv-Zakai  Lower  Bound  (ZZB ) 


The  Ziv-Zakai  bound  is  derived  by  comparing  the  estimation  problem 

with  the  optimal  detection  problem.  The  resulting  bounds  are 

independent  of  bias  and  explicitly  include  the  dependence  on  the  a 

priori  interval  of  the  desired  parameter  [10:386].  The  estimation 

decision  rule  will  have  an  associated  probability  of  error  (P  )  that 

e 

will  be  lower  bounded  by  the  Pq  associated  with  an  optimal  detection 
scheme.  This  lower  bound  (estimation  decision  vs  optimal  decision) 
enables  the  derivation  of  performance  bounds  that  are  based  upon 
detection  theory., 

Consider  an  estimation  technique  for  the  target  angle  off  beam- 
axis.  et,  of  a  received  message  when  it  is  known  that  the  angle  is 
either  9Q  or  .  The  Ziv-Zakai  approach  is  to  compare  the  angle 

A  A 

estimate.  9,  with  an  average  value  (ie  compare  9  with  (9Q  +  5.^/2}  and 
then  decide  between  3  and  9.. 


Figure  15.  Estimation  Decision  Space 
Source:  (10:387) 


Fron  the  law  of  total  probability: 

A  A 

Pe(est)  -  p{e0)p{<e  >  aQ  +  a,)/®,,!  +  Pte^pue  <  eQ  +  e^/e^ 

2  2 

A  A 

-  1/2P{{9  >  9q/2  +  91/2)/9q)  +  1/2P{(9  <  9Q/2  +  9^/2)/*^ 

where : 

P(9q)  -  P{91) 

conditioning  on  9^  assumes  s(9Q)  was  transmitted 
Let:  d  -  91  -  9 

A  A 

P  test)  -  1/2P{ (9  >  9q  +  d/2 )/9q}  +  1/2P{(9  <  »,  “  &/2)/9y) 

A  A 

•  1/2PU9  -  *0  >  d/2)/90J  +  1/2P{ (9  -  9^  <  - d/2)/9y ) 

A  A 

-  1/2P{  ( 1 9  -  9q|  >  d/2)/«0)  ♦  1/2P{(|9t  -  9 1  >_  A/2)/9y) 

Recall  that  P  must  be  greater  than  an  optimal  detection  scheme,  (ie 
for  equiprobable  binary  signaling): 

Pe  1  V2P{  ( 1 9  -  9Q1  >  d/2  )/9fl)  ♦  1/2P{(|9  -  9j  >  &/2)/9y)  (37) 

where 

I ®  i I  -  -  «l 

To  proceed,  the  right  side  of  equation  37  must  be  simplified. 
Reference  9  applied  the  Tchebychef f  inequality  to  the  right  side  of 
equation  37.  The  Tchebychef f  inequality  is  a  specific  case  of  a  more 
general  rule  which  is  termed  the  inequality  of  Bienamye’  (11:151): 

Let  Y  "  |x  •'  aj”;  since  Y  is  always  positive: 


Applying  equation  38  to  equation  37: 


Let:  x  *  8}  d/2  *  b;  a  *  8  or  8^j  n  *  2 

P{|8  "  e0l  1  d/2}  £  E{l®  '  80|2}/<d/2>2  ( 38A) 

PC  i  0  -  ej  _>  d/2)  £  EC  |  e  -  81|2}/{d/2)2  { 38B) 

Using  equations  38A  and  38B  in  the  inequality  of  equation  37: 

Pa  <_  1/2EC | 0  -  9o|2/0o)/(d/2)2  ♦  1/2E{ | 8  -  0 , | 2/9^/W2)2 
ld2/4]Pe  _<  1/2tE{{0  -  80)2/80)  ♦  E( ( 8  -  •  j)2/^)! 

[d2/4]Pe  <  1/2te2{80)  +  e2 ( 8 1 ) J  (39) 

where : 

*  A- 

2  2 

E{(9  -  0)  /8)  «  e  (0)  *  mean-square  estimation  error 
To  simplify,  assume  0Q  *  -8  : 

[(28)2/4]P  (0.-0)  <  1/2  (E(  ( 0  -  8 )2/0)  ♦  E{ (9  -  9)2/-9}J  (40) 

e  — 

where: 

P  (9,-9)  «  error  probability  of  the  best  procedure  for 
© 

deciding  whether  a  target  is  at  9,  or  -9,  when  it 
is  known  to  be  at  one  of  these  positions  with  equal 
probability 

02Pb(0,-9)  £  1/2e2(0)  ♦  1/2e2(-0)  (40A) 

The  right  hand  side  of  the  inequality  of  equation  40A  is  a  lower, 

A5 

bound  to  the  arithmetic  average  of  the  mean  square  estimation  error  (e  ) 
for  any  pair  of  the  values  of  the  parameter  8  which  are  20  units  apart 
(10:387).  imposing  the  symmetry  of  the  mean-square  errors  equation  40A 
then  becomes: 

-36- 


>%rf.V.V.V 


2  2  2 
8  P  ( 8, -€  )  <  1/2e  (8)  +  1/2e  (8) 
e  — 

“2  2 

e  >  8  P  (8,-8) 

—  e 


Many  bounds  on  the  estimation  error  can  be  derived  from  this  basic 
result.  Clearly,  the  worst  case  error  will  be  lower  bounded  by  the 
M»ima  in  the  a  priori  interval  of  the  parameter  8^  [12:650).  In  this 
case,  the  estimation  error  is  bounded  from  below  by  the  scaled  maximum 
of  the  product  of  the  square  of  the  difference  of  the  two  possible 
angles  and  the  minimum  achievable  detection  error  for  equiprobable 
binary  signaling  over  the  a  priori  interval. 


0  <  8  <  8 
t  o 


•t  VV-V 


(40B) 


Reference  9  developed  an  extension  to  the  bound  on  the  mean-square  error 
for  bearing  estimation  given  in  equation  40B  for  the  specific  case  of 
maximum-likelihood  estimation.  The  development  is  as  follows  [9:156): 


A  tighter  bound  can  be  obtained  for  the  catse  where  the  probability 
density  function  is  a  function  only  of  the  difference  between 
estimated  and  actual  parameter  values.  This  constraint  is 
satisfied  by  a  maximum-likelihood  estimator  u  *  sin(S).  Realizing 
that  |sin(8Q)  -  sinfB^I  £  1 8^  -  8^J,  the  maximum-likelihood 
estimator  can  be  applied  to  give  the  appropriate  bounds  on  the 
mean-square  error: 
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(40C) 


2  2 
e  >  max  (sin(8.)]  P  ( 0. , -9  ) 

“  0  <  0  <  0  t  e  t  t 

-  t  -  p 

where: 

sin{8  )  *  lsin(8  )  +  sin(8)]/2 
p  max 

0  =  maximum  possible  anglt  (a  priori  interval) 

max 

Equation  40C,  [3:156],  will  be  the  form  of  the  Ziv-Zakai  bound  used 
to  predict  the  mean-square  estimation  error  for  a  maximum-likelihood 
estimate  of  azimuth  angle.  To  apply  equation  40C,  the  error 
probabiltiy,  P^,  must  first  be  computed. 

HQ  and  H1  are  the  two  hypotheses  for  the  binary  detection  problem 
of  whether  the  azimuth  angle  of  the  target  is  at  0Q  or  0^,  when  it  is 
known  that  a  target  is  at  one  of  these  two  positions.  The  two 
hypotheses  are: 

HQ :  R  “  E{R/Hg)  +  N 

H  j :  R  -  E(R/H^)  +  N 

The  total  probability  of  error  is: 

P@  P(making  an  incorrect  decision) 

*  P{ deciding  HQ  but  H1  is  true  or  deciding  H1  but  HQ  is  true) 

-  P(H^)  ^pR(r/ H0)dr  +  P{H0)  ^p^r/H^dr 

where: 

H1  integral  is  area  under  "0"  density  in  the  "1*  region 
Hg  integral  is  area  under  "1"  density  in  the  "0*  region 


The  decision  rule  is  then: 

I  >"’  ><v 

r  /  n  \  w 


where : 


L(R)  is  the  likelihood  ratio  (5:26J  -  ptr/H^ 

P(r/H0) 

equal  a  priori  probabilities  is  assumed 
equivalently: 

1(R>  0 

_ HO^ _ 

where : 

1(R)  is  the  natural  logarithm  of  the  likelihood  ratio 


ZERO 

THRESH 

♦  ♦  ♦ 

♦  Pe<H . >  ♦ 

-  - Pe(H0)  * 

♦  ♦  +  + 

♦  ♦  ♦  ♦  ■ 

*  *  *  * 

♦  ♦  • .  ♦  ♦ 

•  •  •  m  1  • 

♦  ♦  ♦  ♦ 

•  •  «*  *  0 

Figure  16.  Decision  Space  for  1(R) 


From  Figure  16,  the  total  probability  of  error  is: 


P  =*  1/2P(1(R)  >  0/Hn)  +  1/2P{  1(R)  <  0/H.) 
e  u  1 


(41) 


The  log- likelihood  ratio  is  a  random  variable.  This  is  easily  seen 
from  the  derivation  for  the  likelihood  ratio  which  involves  the  ratio  of 
two  functions  of  a  random  variable;  therefore,  X,(R)  is  a  one-dimensional 
random  variable  [5:26],  and  the  logarithm  of  a  random  variable  is 
random.  With  Gaussian  statistics  for  each  of  the  joint  conditional 
densities,  1(R)  is  the  stun  of  independent  Gaussian  random  variables; 
therefore  Gaussian.  Computing  the  statistics  of  the  one-dimensional 
Gaussian  random  variable  1(R)  will  allow  the  development  of  the 
probability  of  error  as  described  by  equation  41.  From  equation  24, 
the  joint  conditional  density  when  HQ  is  true  is: 

Pr(  r/H0 )  *  equation  24  evaluated  at  8  -  0Q  (42) 

-  (1)exp(-1/N0[(r1-AVa(90)sqrt(E))2  +  (r2-AVd(®0)sqrt(E) )2)  )  (42A) 

(pi)N0 

Similarly,  the  joint  conditional  density  when  is  true  is: 

-  ( 1)exp(-1/M/.[(r.-AV  (9.)sqrt(E))2  +  (r^-AV^S.  )sqrt(E)  )2]  )  (42B) 

0  l  e  1  2  ai 

(pr)N0 


Then  1(R)  is  the  logarithm  of  equation  42B  divided  by  equation  42A: 

1  ( R)  *  -l/NQ(r1-a)2  -l/NQ(r2-b)2  tl/N^r^c)2  +1/n0(r2-d)2  (43). 

where: 


a  -  E{R,/H.)  «  AV  (9,)sqrt(E) 
ll  si 

b  -  E(R_/H . )  *  AV  (0. )sqrt(E) 
2  1  a  1 

C  »  E(R,/H. }  «  AV  (9rt)sqrt(E) 
10  s  o 

«  AVd(90)sqrt(E) 


d  -  E(R2/H0) 


1<R)  --1/H  (r,2  -  2r^  +  a2)  -1/NQ(r22  -  2i£>  +  b2) 

+1/UQ(r12  -  2r1c  +  c2)  +1/NQ(r22  -  2r2<3  +  d2) 

=  1/N0[2r1(AVs(81)sqrt(E)  -  AVg(0Q)sqrt(E) )  + 

2r  ( AV  ( 0 . ) sqrt ( E )  -  AV.(8.)sqrt(E) )  - 
2  a  1  a  u 

A2EVg2(81)  +  A2EVg2(80)  -  A2EVd2(81)  +  A2EVd2<0o)] 

-  2Asqrt{2)/No[r1(Vg(01)  -  Vg(80))  +  r^V^)  -  V.^))  + 
Asqrt(E)/2(Vg2 (8g )  -  Vg2(01)  +  Vd2(0Q)  -  V^2 ( 0 , ) 1  (43A) 

Introducing  a  new  version  of  1(R): 

1*  (R)  »  1  ( R)  /  ( 2Asqrt ( E ) /tIQ )  (43B) 

The  statistics  of  l'(R)  for  the  hypothesis  HQ  are: 

ECr(R)/H0>  -  E{R1(Vg(01)  -Vs(0o)/Ho)J  +  E{R2(Vd(01)  ♦ 

Asqrt(E)/2(Vg2(0Q)  *  Vg2(01)  +  Vd2(0o) 
-Asqrt(E)[Vs(0o)(Vg(01)-Vs(0o))+Vd(0o){Vd(01)-Vd(0o))]+ 
Asqrt(E)/2(Vg2(0o)  -  Vg2(81)  +  V^COq)  -V^tS,)) 

-  Asqrt(E)[Vg(0o)Vg(01)  +  Vd(0o)Vd(01)  - 

Asqrt(E)/2 ( (Vg2( 0Q)  +  V#2(0t)  ♦  Vd2(0Q)  ♦  Vfl2 { 0 n ) ] 

«  Asqrt(E)W/2 Ip  -  1.0]  (44) 

where: 

W  -  V  2 { 0* )  +  V  2(0.)  +  V  2(0ft)  ♦  V  2{0.)  (45) 

s0  si  d  o  ai 

p  -  2/WtVg(0o)Vg(81)  +  VVW1  U6> 

var(I'(R)/H0)  *  var(R1/H(J)  [v^e^  -  Vg(0o)]2  + 
var(R2/Ho)[Vd(01)  -  Vd(0o)]2 

-  V2lvs2(9o}  +  vg2(V  *  V l#o>  +  vd2(ein  * 
VVVVV  +  VVVV1 
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(47) 


var(l'(R)/H0}  =  NqW/2 (1.0  -  p] 

A  similar  development  for  hypothesis  yields: 
E(1'(R)/Ht)  =  Asqrt(E)W/2 [ 1 .0  -  p] 
varU'tRJ/H^  =  NQW/2[1.0  -  p] 


The  conditional  distribution  under  the  hypothesis  is: 

pdMW/H^  =  kexp[-(l'  (R)-E(X'  (R)/H1 )  )2/2var{l*  (R)/H^)] 
where: 

k*  [sqrt(2(pi)var(l*  (RJ/H^  )J_1 
s;i'(R)/H  )  *  Asqrt(E)W/2 [1.0  -  p] 

var{l'  (R)/H1 }  *=  NQW/2[1.0  -  pj 

Similarly,  under  the  hypothesis  HQ: 

p(l'(R)/H0)  *  kexpi-(l’(R)+E{l,(R)/H1})2/2var(l'{R)/H1)] 
where : 

E(1,(R)/H0)  *  -  K(l'(ft)/tt  } 

var(l' (R)/Hq}  *  var{l*  (RJ/H^ 


Equation  41  then  becomes: 

«  1/2P(1'(R)  >  0/Hq)  +  1/2P{1'(R)  <  0/H^ 

-  1/2P{1'(R)  >  0/Hq)  +  1/2P(1*(R)  >  0/Hq} 

»  P{1'(R)  >  0/Hq} 
where : 

Due  to  the  symmetry  of  the  shifted  distributions: 


(48) 

(49) 


Pe  -  Jp( l*(R)/HQ)dl'(R) 

-  Jkexpt-(l,(R)+E{l,(R)/H1J)2/2var{l,(R)/H1))dl'(R)  (50) 

where 

integrations  range  from  zero  to  infinity 

Let  u  -  (l,(R)+E{l,(R)/H1})/sqrt{var{l'(R>/H1}) 

du  *  dl'  (R)/sqrt(var{l'  (RJ/H^  ) }  dl*(R)  *  sqrt(var(l*  (RJ/H.j)  )du 
when  1*(R)  =  0  ;  u  *  E(l* (RJ/H^J/sqrtfvartl* (RJ/H^J ) 

Equation  50  then  becomes: 

co 

Pe  *  l’(R)*0 J  •q*t(2(pi)>"“1«icp[-a2/2]du 

*  Q(u  when  l'(R)  *  0| 
where: 

Q(al  «  sqrt(2(pi) ) 

Pe  -  Q[E(1'  (RJ/H^/sqrttvartl'  (RJ/H^J)) 

»  Q [Asqrt( E)W/2 (1.0  *  p) /sqrt(NQW/2 ( 1 .0  -  p] ) 

«  &[Asqrt( [EW/2Nq] (1.0  -  p))]  (51) 

Equation  51  describes  the  necessary  P@  of  the  best  procedure  for 
deciding  whether  a  target  is  a  9Q  or  9^,  when  it  is  known  to  be  at,  one 
of  these  two  positions  with  equal  probability.  For  the  case  at  hand,  *Q 
**  -9 1  and  substitution  of  equations  9  and  10  into  equations  45  and  46 
the  factors  w  and  p  become: 


-  (vs2‘V  *  vs2‘V  +  va2(9o’  *  va^ * ® i * 

»  1/2 [ 1+2cos( 2 .18/BW)  +  cos2 (2. 18/BW)]  +  sin2( 1 .90/BW)  + 

1/2  [  1+2cos(  -2 . 18/BW)  +  cos2 (-2. 18/BW)]  *  sin2(  -1 .90/BW) 

=■  5/2  +  2cos( 2 . 18/BW)  +  1/2cos( 4 .29/BW)  -  cos(3.89/BWj  (52) 


Figure  17  provides  a  plot  of  the  factor  W  out  to  three  beamwidths. 


Figure  17.  W  vs  Normalized  Beamwidth 


p  *  2/W[Vs(80)V3(81)  ♦  v<J(®0)Vd(®1)l 

*  [1/2  +2cos(2. 18/BW)  ♦  1  /2  CQS( 4 . 2 9/BW )  +COS(3.88/BW)]/W  (53) 


A  plot  of  the  parameter  p  is  provided  in  Figure  18.  Inspection  of 
Figure  18  shows  the  parameter  p  to  be  periodic  in  normalized  beamwidth 
and  approximately  equal  to  cos(2.228/BW) . ,  F.  /are  18  displays  only  a 
small  window  of  the  function  p  because  the  particular  form  of  equation 
52  has  discontinuities  at  multiples  of  1.7  times  the  normalized 
beamwidth.  No  approximation  for  p  is  utilized  because  all  angles  of 


interest  are  within  plus  or  minus  one-half  the  beamwidth.  With  the 
magnitude  of  p  always  less  than  or  equal  to  1,  p  satisfies  the  condition 
of  a  correlation  coefficient. 


Figure  18.  p  vs  Normalized  Beamwidth 


Inspection  of  Figures  17  and  18  reveal  that  when  the  parameter  8  is 
zero.  W  equals  4  and  p  is  unity.  When  p  is  unity,  equation  51  is 
maximum;  and  therefore  P^  is  maximum.  This  result  is  intuitively 
appealing  as  targets  with  small  seperation  angles  will  be  difficult  to 
distinguish.  Substitution  of  equation  Si  into  equation  40C,  along  with 
equations  52  and  53,  will  provide  the  Ziv-Zakai  lower  bound  for  the 
estimation  error  of  a  maximum-likelihood  estimate  of  the  target  azimuth 
angle. 
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A  '  A 

e  *  R,(0)V,<0)  -  R,<0)V  (0)  (28) 

la  2  s 

To  provide  automatic  tracking  of  targets,  the  error  signal  actuates 
a  servo-control  system  to  steer  the  beam-axis  on  target.  If  the  error 
signal  does  not  change  rapidly  with  respect  to  the  operating  frequency 
of  the  radar,  the  error  will  be  approximately  a  dc,  or  step,  input  to 
the  antenna  servo  system.  The  time  domain  response  of  the  servo  is  then 
approximately  the  response  that  would  be  obtained  due  to  a  scaled  step 
function  proportional  to  the  error  signal.  In  that  light,  the 
development  for  the  tracking  loop  will  be  as  follows: 

1)  A  relationship  will  be  derived  between  the  error  signal  and  a 
step  function  that  is  proportional  to  the  difference  of  the 
estimated  target  location  and  actual  target  position.  This 
conversion  will  be  in  the  form  of  a  discriminator  curve  that 
will,  within  limits,  provide  a  1  to  1  mapping  between  error  and 
a  step  function  with  amplitude  proportional  to  required 
boresight  movement. 

2) .  Since  the  error  signal  will  be  in  the  form  of  a  step  function, 

an  improved  Type  I  servo  system  will  be  developed  to  provide 
the  desired  response  time  and  steady  state  error  for  the 
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monopulse  tracker.  A  Type  I  servo  gives  zero  steady  state 
error  for  a  step  input,  and  an  improved  Type  I  provides  the 
desired  response  time. 

3)  Relationships  will  be  established  for  tracking  moving  targets. 

Discriminator 

With  additive  noise,  the  error  equation  is  a  random  quantity  as  R^ 
and  R^  are  independent,  Gaussian,  distributed  random  variables.  As 
stated  earlier,  the  discriminator  curve  is  established  assuming  a 
perfect  mapping  from  error  signal  to  a  step  function  required  to  align 
the  beam-axis  with  the  target  position  (ie  the  mapping  is  accomplished 
in  the  no-noise  environment).  Zn  this  case,  the  observations  are  given 
by  the  conditional  means: 

R  «  E{R( 9 )/A, 8} 


*  AVg(8)sqrt(E)  for  observation  R^  (54) 

*  AV  ( 8 ) sqrt ( E )  for  observation  R  (55) 

Q  2 

Substitution  of  equations  54  and  55  into  equation  28  yields: 

* 

e  •  AV  (8)sqrt(E)V,(8)  -  AVJ(8)sqrt(E)V  (8)  (56) 

s  a  a  s 

Recall  that  V  (8)  and  V. (8)  were  described  by  equations  9  and  10: 
s  d 

V  (8)  -  sqrt(2)/2 ( 1  +  cos(2.O940/BW)J 
s 

V^te)  *  sin(1.88e/BWl 
a 


u 


¥ 


As  previously  stated,  it  is  necessary  to  express  the  error  equation 
(equation  56)  as  a  function  of  the  difference  between  estimated  and 
actual  target  location*  Such  a  function  allows  boresight  movement 


•  *  * *s" •  * *N •  •  ’**..** •  ■  •-< v 


*  /  w 


proportional  to  this  difference.  If  all  anqles  are  considered  relative 
to  the  boresight  axis  the  radar  will  estimate  a  target  location,  the 
target  may  move  relative  to  that  location,  and  the  boresight  will  then 
move  to  the  estimated  position  where  the  whole  procedure  is  repeated 

A 

over  again.  In  this  light,  each  angular  difference  between  9  and  8  is 

A 

approximately  equal  to  the  difference  between  0  and  the  boresight  (ie  0 
is  approximately  equal  to  zero) . 

A  A  A 

Let  d9  *  9  0»  0  *  d0  +  8;  For  small  0: 

A 

V  (0)  ■  sqrt(2)  (56A) 

s 

A  A 

V.(8)  *  1.880/BW  which  is  approximately  0  (56B) 

d 


Substitution  of  equations  56A  and  56B  into  equation  56: 
e  -  -  Asqrt ( 2 ) sqrt ( E) [sin( 1 . 886/BW) ] 

-  -  Asqrt (2E) [sin( 1 ,88d0/BW) ] 
d0  -  -  BW/1 .88arcsin[e/Asqrt(2E)J 


(57) 

(57A) 


Equation  57a  is  the  transformation  required  from  an  error  signal  to 
a  step  function  with  amplitude  proportional  to  the  required  boresight 
movement  for  target  tracking. 

Figure  19  provides  a  plot  of  equation  57.  Due  to  thfr  inverse  sine 
function  associated  with  the  discriminator  curve,  all  values  of  error  do 
not  map  to  a  unique  value  of  d9.  Limits  must  be  place  on  d0  such  that 
within  the  defined  interval,  there  is  approximately  a  one  to  one  mapping 
between  d9  and  the  error  signal.  For  large  errors,  d0  can  then  only 
move  to  an  established  maximum,  and  no  further.  This  constraint  doesn't 
limit  the  usefulness  of  the  discriminator  curve,  because  in  a  realistic 
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sense  the  mechanical  limitations  on  how  far  the  bores ight  can  move  for 
large  inputs  may  be  the  final  performance  limitation.  The  limits  used 
for  d8  are  noted  on  Figure  19.  Equation  57A  will  be  then  be  used, 
within  the  specified  limits,  to  actuate  a  servo-control  system  to  steer 
the  beam-axis  on  target. 


Figure  19.  Discriminator  Curve 


Angle  Servo 

A  tracking  radar  can  be  divided  into  two  pirts.  The  first  part 
estimates  the  position  of  the  target  within  the  resolution  cell.  The 
second  part  centers  the  resolution  cell  on  the  target,  typically  with 
the  aid  of  servomechanisms  (4:3121.  If  gear  backlash,  and  other  non- 
linearities,  are  neglected,  the  servo  system  can  be  modeled  as  a  linear 


system.  Analysis  of  this  linear  system  can  be  performed  with  the  aid  of 
Laplace  transforms. 
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C(s)  *  - ° - 

*Is+Oj) 

H(s)  *  1 
K  *  *0** 

Figure  20.  Type  I  Angle  Servo 
Source:  [4:313] 


Figure  20  shows  a  Type  1  angle  servo  and  associated  parameters. 

The  corner  frequency,  w  ,  is  fixed  by  antenna  inertia  and  motor  torque 
[4:3141.  Combining  the  gains,  Xp  and  Km  of  Figure  20,  results  in  an 
open  loop  transfer  function: 

G(s)  -  K/s [s  +  w^J 

Closing  the  loop,  the  closed  loop  transfer  function  is: 

®out/0in  *  H(s)  "  *  GCs)l 

■  K/(s(s  +,  wy)  ♦  KJ 

-  K/ls2  +  sw1  ♦  X]  (58) 

The  characteristic  equation  of  equation  58  is: 

2 

s  +  SWj  ♦  K  -  0  ( 58A) 

The  roots  of  equation  58A  determine  the  time  response 
characteristic*  of  the  linear  system.  Exciting  equation  58  by  the 


Laplace  transform  of  a  unit  step  function  (1/s),  separating  into  partial 
fractions,  taking  the  inverse  Laplace  transform,  and  then  simplifying 
the  result  yields  the  time  domain  output  response  due  to  a  unit  step 
input.  For  t  0,  this  response  can  be  shown  to  be: 

0  (t)  *  1  -  exp(-fw  t) [sin(w  sqrt(1  -  f2)t)  + 

out  n  n 

cos(w  sqrt( 1  -  f2)t)  /sqrt( 1  -  f2)]  (59) 

n 

where : 

f  »  w^/2sqrt(K)  *  damping  coefficient 
w^  ■  sqrt(K)  “  system  natural  frequency 

Since  w^  is  fixed  by  the  system  elements,  desired  system 
performance  (ie  allowable  settling  time,  steady  state  error,  and 
percent  overshoot  in  the  output  response  to  a  unit  step  function)  must 
be  established  by  the  adjustment  of  the  gain  parameter,  K.  Settling 
time  is  the  time  required  for  the  system  to  damp  out  all  transients.  In 
practice,  the  settling  time  is  established  where  the  error  (ripple  about 
the  desired  response)  is  reduced  below  5%  of  the  initial  error  [13:90]. 
For  approximation  purposes,  the  settling  time  is  within  4  time  constants 
of  the  envelope  of  the  damped  sinusoidal  oscillation.  The  ratio  of  the 
step  response  peak  value  to  the  steady  state  settling  value  is  termed 
percent  overshoot.  The  amount  of  overshoot  allowable  depends  upon  the 
particular  system,  and  is  a  function  of  the  damping  ratio,  but  ten 
percent  over  shoot  is  reasonable  (13:91).  For  the  mohopulse  radar,  the 
constraints  for  the  servo  are: 

1)  Thfi  servo  response  to  a  pulse  must  settle  to  steady  state 

before  the  arrival  of  the  next  pulse.  This  fixes  the  settling 


time  to  be  less  than  the  PRI,  which  in  turn  fixes  the  rise  time 
of  the  servo  1.0  be  less  than  one-fourth  the  PRI. 

2)  To  maintain  allowable  overshoot,  the  damping  ratio  must  be  set 
to  approximately  1/sqrt(2)  for  optimum  second  order  response  to 
a  step  input  [4:315]. 

Due  to  the  pulse  repetition  frequency  (PRF)  of  the  radar  (lOOpps), 
the  Type  1  servo,  as  described  by  equation  59,  could  not  settle  within 
the  pulse  repetition  interval  (PRI)  and  a  slightly  different  servo 
system  than  the  Type  I  servo  must  be  employed. 

Higher  gain  can  be  utilized  while  holding  the  damping  ratio  to 
approximately  one-half  with  the  addition  of  a  reciprieal  corner  (phase- 
lead  network)  in  the  system  closed  loop  transfer  function  of  a  Type  X 
servo  system  [14:251].  This  phase-lead  network  eliminates  the  overshoot 
problem  at  high  values  of  gain  while  reducing  the  step  response  rise¬ 
time  [4:320]  .  The  step  response  of  the  improved  Type  X  servo  is  given 


in  Reference  14  as  [14:254]: 

0  (t)  ■*  [1  -  exp(-fw  t)/sqrt{1  -  f2)]sqrt(1  -  w  /w  ) 

out  n  in 

[cos(w  sqrt( 1  -  f2)t  -  0)]  (60) 

n 

where: 

w^  ■  antenna  corner  frequency  -  10  rad/sec  [4:321]  (60A) 

w2  -  phase  lead  network  break  frequency 

wn  «  natural  frequency  *  sqrt(K)  (60B) 

f  »  damping  coefficient  •  K/(2w  w_)  +  w,/(2w  )  (60C) 

n  z  in 

0  *  arctanK  1/2sqrt(  1  -  f2))(w  /w  *  w/w)]  (60D) 

1  n  n  2 
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The  rise  time  is  provided  as  [14:254]: 

(pi)  -  arctan[( 1/2sqrt( i  - 

,  t  *  - 

r  w^sqrtd  - 

To  establish  performance  at  the  desired  level,  the  gain  (sqrt(K) 
and  the  phase-lead  break  frequency  (w^)  were  adjusted,  by  am  iterative 
technique,  with  constraints  of  a  fixed  damping  ratio  ( f  -  0.707)  and  a 
rise  time  that  must  be  slightly  less  than  one-fourth  the  PRI  (tf  < 
2.5msec).  For  the  iteration,  the  ratio  w^/w^  was  held  constant  at  1.4 
until  an  acceptable  value  of  rise  time  was  attained,  and  then  w^  was 
adjusted  until  f  *  0.707.  All  calculations  were  performed  on  a 
programmable  hand-held  calculator.  The  iterative  results  are: 


wr  •  900  rad/sec  (61A) 

w2  »  641  rad/sec  (61B) 

From  equations  6ilA  and  60B: 

K  «  w  2  »  810E+3  (62) 

n 

Substituting  equations  62,  61A,  60A,  and  618  into  equation  60C: 

f -  0.707  (63) 

Using  all  of  the  constants  and  solving  equation  61: 

tr  *  2.48E-3  sec  (64) 


Solving  for  the  step  response  (equation  60): 

*  „<t)  *  M  *  exp< -fw  t ) / sqrt ( 1  -  f2 ) ] sqrt(  1  -  w  /w  ) 

ouu  n  in 

[cos(wnsqrt( 1  -  f2)t  -  f)] 

«  1  -  1 .39exp( -6J6t) (cos(636t  -  0.778)) 


9  (t)  =  1  -  1.39exp(-636t) [cos(636t)cos( .78)  -  sin(636t)sin( .78) ] 

out 

=  1  -  .99exp( -636t) [cos(636t)  -  .99sin(636t) 1 
=  1  -  exp(-636t) [cos(636t)  -  sin(636t)J  (65) 

Equation  65  is  the  output,  time  domain,  servo  response 
characteristic  for  the  monopulse  tracker.  An  amplitude  scaled  version 
of  equation  65  is  used  to  steer  the  beam-axis  on  target  for  a  given 
error  signal  input.  A  plot  of  equation  65  is  provided  in  Figure  21  for 
a  unit  step  input. 


Figure  21.  Response  of  Improved  Type  1  Servo  to  a  Unit  Step  Function 


Moving  Target 

Several  moving  target  considerations  must  be  accounted  for  in  the 
radar  model. ,  Among  the  considerations  are  range  and  signal-to-noise 
ratio  (SNR)  changes  because  of  target  movement.  Consider  linear  target 
translation  as  shown  in  Figure  22. 


Figure  22.  Simple  Linear  Target  Translation 


From  Figure  22,  a  simple  relationship  can  be  established  for  the 
new  target  range  as  a  function  of  velocity,  time,  and  previous  range. 
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The  range  and  signal-to-noise  ratio  are  related  by  the  radar  range 
equation  [4:29] : 


P  G  2cr  X2 
r  r  t 


Wpill  Hlnt(S/S)p 


where : 


Rm  *  maximum  range  for  a  given  set  of  radar  parameters 
Pr  *  actual  power  deceived  by  the  radar 

G  *  antenna  gain 

X  »  wavelength 

Nint  «  receiver  internal  noise  power 

(S/N)p  «  signal-to-noise  ration  per  pulse 
<r  « 


target  cross-section 


Solving  equation  67  for  (S/N)  : 

P 

P  G  2cr  A2 

{  s/n  )  =  - - 2 

»i„A 


Using  equation  67A  and  after  cancellation  of  like  terms,  the  ratio  of 
new  SNR  to  the  previous  SNR  is: 


SNR(new)  M/R  (new) 


SNR(old)  I 1/R  (old) 
k  m 


R  (old) 


R  (new) 

RL 


or  equivalently: 


SNR(new)  =  SNR(old) fR  (old)/R  (new;) 

m  in 


where: 


R  (new)  *  sqrt [ R  2  +  (vt)2] 

m  oia 

R  (old)  is  established  by  the  initial  conditions 
m 


To  establish  initial  conditions,  some  constraints  must  be  placed  on 
the  simple  translation  described  by  Figure  22.  Let: 
aircraft  velocity  *  500  knots  »  257.4  m/ sec 
antenna  beamwidth  *  3  degrees 

aircraft  travel  one-half  the  beamwidth  in  50  pulses 
radar  PRI  *  10E-3  sec 


After  50  pulses,  the  target  will  travel  50  times  the  PRI  times  the 
velocity,  meters.  Plugging  in  the  constraints,  this  translation  is 
128.7  meters.  Figure  23  provides  the  geometry  of  the  described 
situation  that  now  allows  for  initial  range  computation. 


-aVk’.'.’.  V.'-  : 


Solving  for  the  intial  target  range  yields: 

Rinit  *  1128.7  ml/tan(1.5  degrees)  *  4914.9  meters 

Figure  24  depicts  the  complete  amplitude-comparison  monopulse  radar 
development  that  will  be  utilized  in  the  radar  simulation  program  that 
follows. 
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III.  RADAR  SIMULATION 


Computer  Program 

To  compute  the  mean-square  tracking  error  and  associated 
performance  bounds  of  the  amplitude-comparison  monopulse  radar,  a 
simulation  program  was  designed  for  the  radar  model  as  developed  in  the 
previous  chapter.  An  annotated  program  listing  of  the  program  is 
provided  in  Appendix  B.  The  input  variables  of  the  simulation  allow  for 
choice  of  a  stationary  or  moving  target  and  desired  signal-to-noise 
ratio  (SNR)  level.  The  output  parameters  are  the  mean-square  tracking 
error,  the  Cramer-Rao  bound  for  the  conditional  variance  of  the 
estimated  target  position,  and  the  Ziv-Zakai  lower  bound  for  the 
estimation  error, of  the  target  position.  Two  external  International 
Mathematical  and  statistical  Libraries  (IMSL)  routines  are  called  during 
the  program  execution.  The  first  routine  is  a  Gaussian  random  deviate 
generator,  GGNML,  that  creates  random,  independent,  noise  samples.  The 
second  routine,  HONOR,  computes  the  area  under  a  Gaussian  curve  and  is 
utilized  to  compute  the  Pa  for  the  Ziv-Zakai  bound. 

Average  tracking  error  as  conducted  in  the  simulation  program  is 
developed  by  observation  of  the  radars  ability  to  maintain  the  boresight 
in  the  vicinity  of  the  target  for  a  finite  interval  of  pulses,  and  then 
repeating  the  experiment  many  times  and  averaging  the  results. '  The 
observation  interval  is  established  at  50  pulses,  and  each  individual 
pulse  is  sampled  once.  From  that  one  sample,  an  error  signal  is 
generated  proportional  to  the  difference  between  the  boresight  and 


59- 


target  position.  The  error  signal  is  then  applied  to  the  angle  servo  to 
move  the  boresight  in  the  target  direction.  Each  set  of  50  pulses  is 
considered  a  run.  A  total  of  15  runs  are  performed  and  then  the  error 
is  summed  and  averaged  for  each  of  the  50  pulses  over  the  15  runs.  The 
simulation  program  contains  four  discrete  components: 

1 )  An  input  section 

2)  A  tracking  section 

3)  Computation  of  performance  and  bounds 

4)  An  output  section 

Input 

The  input  is  a  short  section  that  prompts  the  user  for  desired  SNR 
level  and  whether  the  target  is  moving  or  stationary,  establishes  the 
radar  and  program  parameters  (pulse  width,  pulse  repetition  frequency 
(PRF),  beamwidth,  etc.),  and  positions  the  target  in  the  beamwidth  at 
the  beginning  of  a  run.  A  few  comments  about  the  input  section: 

1 )  Two  separate  1X50  vectors  of  zero  mean,  Gaussian,  random 
variables  are  established.  The  two  seeds  are  utilized  to 
insure  independent  noise  samples. 

2)  Provisions  are  provided  for  the  testing  of  15  signal-to-noise 
ratio  levels  in  the  range  from  -20  to  +35  dB. 

3)  The  only  exit  for  the  simulation  is  the  entering  of  a  0  to  the 
prompt  as  to  whether  the  target  is  moving  or  stationary. 

4)  Target  placement  provides  for  translation  to  be  one-half  the 
beamwidth  at  the  end  of  the  50th  pulse  for  the  initial  range, 
velocity,  and  PRP  provided. 
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5)  Signal  amplitude  was  adjusted  while  monitoring  boresight 
movement  for  a  stationary  target  at  changing  SNR  levels.  A 
compromise  between  the  amount  of  pulses  necessary  to  move  to 
the  target  and  target  tracking  potential  over  the  established 
range  of  SNR  levels  fixed  the  amplitude  at  5. 

A  flow  chart  of  the  input  section  is  provided  in  Figure  25. 
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Figure  25:  Program  input  Section 


Tracking 

The  tracking  section  moves  the  boresight  to  the  vicinity  of  thr 
target  for  each  received  pulse.  The  radar  tracks  the  target  for  a  50 
pulse  interval  (run).  During  the  run.  noise  samples  are  generated  from 


the  stored  GGNML  vectors,  the  receiver  observations  are  formulated,  an 
error  signal  is  generated  from  the  observations,  and  the  boresight  is 
moved  in  a  direction  to  reduce  the  error.  Figure  26  gives  a  flow  chart 
of  the  tracking  section. 


Figure  26.  Tracking  Section 
A  few  comments  about  the  tracking  section: 

1)  Provisions  are  provided  for  the  changing  range  and  SNR  of 
moving  targets. 

2 )  Large  errors  are  limited  to  a  maximum  positive  or  negative 
value  to  model  the  limitations  on  how  far  the  boresight  can 
move  for  large  inputs. 

3)  The  time  parameter  for  the  servo  response  is  the  radar  pulse 
repetition  interval  (PRX). 
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Computation  of  Performance  and  Bounds 

This  section  computes  the  squared  tracking  error,  the  Cramer-Rao 
bound,  and  the  Ziv-Zakai  bound.  The  tracking  error  (angle  between 
boresight  and  target  after  observing  a  pulse)  was  computed  in  the  last 
section.  Squaring  this  result  (in  radians)  gives  the  desired  squared 
error  (SE).  The  Cramer-Rao  bound  is  computed  by  substitution  of  the 
angle  9,  after  pulse  observation,  into  equation  36.  The  Ziv-Zakai  bound 
requires  the  maximum  of  equation  40C  for  all  target  angles  from  0  out  to 

8  .  To  find  the  maximum,  the  following  routine  was  utilized: 

P 

1)  Compute  the  arithmetic  mean  angle  8  (8  varies  from  .75 

P  P 

degrees  to  1.5  degrees ) . 

2)  Segment  8^  into  angle  intervals  with  provisions  for  adjustment 
of  the  number  of  segments  (8  segments  when  8*0,  16  segments 
when  9«BW/2 )  to  provide  more  resolution  at  large  angles 
(approximately  equal  step  sizes)  and  then  evaluate  and  store 
equation  40C  for  each  seqmented  interval.  The  stored  values 
are  then  sorted,  and  the  maximum  is  picked  off. 
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Figure  27  provides  a  flow  chart  of  the  described  sequence 


Figure  27.  Computation  of  Performance  and  Bounds 


Output 

The  output  section  averages  the  squared  tracking  error  (SB).  and 
bounds  over  the  15  runs  and  prints  the  three  output  vectors  (the  average 
of  the  SB  will  be  designated  as  mean-square  error  (MSE)).  A  flow  chart 
of  the  output  section  is  shown  in  Figure  23. 
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Figure  28.  simulation  Output  Section 

Performance  Verification 

Stationary  Target 

The  outputs  of  the  simulation  program  are  plotted  in  the  following 
diagrams.  The  first  plot.  Figure  29,  shows  the  general  trend  of  the 
mean-square  tracking  error  for  a  stationary  target  at  different  SNR 
levels.  As  the  figure  indicates,  there  is  an  increasing  error  for 
decreasing  SNR  levels.  There  is  also  a  general  tendency  for  the  mean* 
square  error  to  be  spread  over  a  larger  range  for  the  smaller  values  of 
SNR.  This  is  a  direct  indication  of  the  increased  difficulty  of 
maintaining  the  boresight  within  the  vicinity  of  the  target. 


Figure  29.  General  Trend  of  the  Stationary  Target  Tracking  Error 

Figure  30  is  a  stationary  target  composite  plot  for  mean-square 
tracking  error  (MSE),  the  Cramer- Rao  bound  (CRB),  and  the  Ziv-Zakai 
bound  (ZZB)  for  SNR  levels  of  15,  7,  -5,  and  -20  dB.  As  the  -5dB  plot 
indicates,  the  CRB  begins  to  approach  the  MSE.  The  -20dB  plot  clearly 
shows. the  CRB  exceeding  the  MSE.  Since  by  definition  a  lower  bound  must 
always  be  less  than  or  equal  to  the  actual  error,  the  CRB  is  not  useful 
at  low  SNR  levels.  Where  the  actual  bound  loses  its  usefulness  depends 
upon  the  system  parameters,  but  for  the  purposes  of  this  simulation,  the 
CRB  for  a  stationary  target  will  only  be  valid  for  SNR  levels  greater 
than  -5dB.  In  contrast  to  the  CRB,  observation  of  Figure  30  shows  the 
ZZB  is  tight  at  the  lowest  measured  SNR  level.  The  -20dB  plot  shows  the 
zcu  to  be  approximately  constant  and  no  greater  than  the  average  MSS 
excursions. 
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Figure  31  provides  the  general  trend  of  the  mean-square  tracking 
error  for  a  moving  target  at  different  SNR  levels.  As  the  figure 
indicates,  there  is  an  increasing  error  for  decreasing  SNR  levels. 
Further,  comparison  with  the  general  trend  for  a  stationary  target 
(Figure  29)  shows  two  significant  results.  The  first  observation  is 
that  for  high  SNR  (15dB)  the  MSE  settles  out  to  the  same  value  as 
stationary  target  tracking.  The  second  observation  is  that  for  the 
lower  SNR  levels,  the  moving  MSE  is  slightly  higher  than  the 
corresponding  stationary  values  and  the  degree  of  MSE  spreading  for  a 
moving  target  is  noticeably  larger  (increased  tracking  difficulty). 


Figure  31.  General  Trend  of  the  Moving  Target  Tracking  Error 

Figure  32  provides  the  moving  target  composite  plot.  As  the  -2QdB 
plot  indicates,  the  CRB  again  far  exceeds  the  MSE.  In  contrast  to  the 
CRB,  the  ZZB  is  once  again  tight,  approximately  constant,  and  no  greater 
than  the  average  MSE  excursions. 
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IV.  TARGET  FLUCTUATIONS 

The  receiver  design  and  analysis  up  to  this  point  considered  only 
deterministic  amplitudes  of  the  received  signal.  Consequently,  the 
performance  estimates,  up  to  this  point,  are  based  upon  a  deterministic 
amplitude .  Of  interest  in  this  chapter  is  what  happens  to  performance 
estimates,  and  to  target  tracking  capability,  when  target  fluctuations 
are  present  that  were  not  designed  for. 

To  properly  account  for  target  cross-sections,  the  probability 
density  function  must  be  known  for  the  particular  type  of  target  that  is 
being  tracked.  The  model  of  the  radar  considered  thus  far  has  utilized 
deterministic  amplitudes  that  are  proportional  to  the  positive  square 
root  of  the  radar's  cross-section.  The  next  extension  is  to  assume  that 
there  are  many  individual  point  scatterers  with  each  individual 
scatterer  contributing  (in  an  additive  fashion)  to  the  cross-section  of 
the  target.  The  probabilistic  model  of  a  target  in  free  space  with  many 
individual  contributers  will  turn  out  appear  as  one  source  which  has 
an  Exponentially  distributed  cross-section  (15s 47].  The  corresponding 
target  amplitude  will  have  a  Rayleigh  distribution;  Inherent  in  this 
development  is  the  assumption  that  none  of  the  individual  scatterers  are 
much  larger  than  the  rest  (point  targets)  and  that  the  cross-section 
does  not  change  during  the  sampling  interval. 

To  account  for  fluctuating  cross-sections,  the  nature  of  cross- 
section  dependence  upon  radar  viewing  angle  must  be  known.  The  target 
cross-section  is  very  dependent  upon  aspect  angle,  and  the  aspect  angle 


r 
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changes  with  time  [4:38] .  If  the  cross-section  changes  significantly 
over  a  time  of  viewing  equal  to  the  pulse  repetition  interval  (PRI),  the 
fluctuations  are  assumed  independent  from  pulse  to  pulse  [15:47].  This 
paper  will  consider  such  fluctuations  to  be  "rapid"  fluctuations.  If 
the  cross-section  changes  are  significant  only  over  times  on  the  order 
of  a  burst  of  pulses,  many  times  the  PRI,  the  fluctuations  are  assumed 
to  be  independent  from  scan  to  scan  [15:46].  This  paper  considers  such 
fluctuations  to  be  "slow"  fluctuations.  Target  fluctuations  will  be 
investigated  in  the  following  steps:  . 

1)  A  variable  transformation  will  be  performed  on  an  Exponential 
(cross-section)  density  to  insure  that'  the  positive  square  root 
(amplitude)  will  be  Rayleigh  distributed. 

2)  The  effect  of  a  random  amplitude  on  the  receiver  estimator  will 
be  investigated. 

3)  The  simulation  program  developed  in  chapter  3  will  be  modified 
to  investigate  tracking  performance  with  "slow"  and  "rapid" 
fluctuations  on  the  signal  amplitude. 

4)  A  simple  method  to  reduce  the  mean-square  tracking  error  for 
target  tracking  in  the  presence  of  amplitude  fluctuations  will 
be  proposed  and  investigated. 

Distributions 

As  previously  stated,  the  probabilistic  model  of  a  target  in  free 
space  with  many  individual  contributors  appears  to  the  radar  to  be  one 
source  that  has  an  Exponential  distribution  on  the  radar  cross-section. 
Let  the  cross-section  be  represen" by  the  random  variable  X,  the 


probability  density  function,  pdf,  of  X  is  given  by  [7:45]: 

fx(x)  *  aexp[-ax]  x  ^  0  (69) 

where: 

a  *  VE(X) 

The  units  of  cross-section  are  units  of  power.  To  convert  to  a 
form  of  positive  amplitude,  power  must  be  converted  to  field  strength. 
The  appropriate  transformation  is: 

Y  «  +sqrt(X)  (70) 

The  transformation  requires  the  following  theorem  [7:118]: 

Let  X  be  a  continuous  random  variable  and  Y  «  g(X)  where  g(X) 
is  continuous  in  X  and  strictly  monotonic.  Then 
fy(y>  -  fxt9‘1(y)]i(g‘1(y)),|  (7t) 

where : 

|  |  denotes  absolute  value 
(  ) *  denotes  the  derivitive  wrt  y 

Utilizing  equations  69  and  70  in  equation  71  and  realizing  that  for 
X  restricted  to  positive  values  and  Y  restricted  to  the  positive  square 
root: 

fx(x)  *  aexp{-ax]  x  ^  0 
g(x)  «  fsqrt(X) 
g^ty)  *  y* 

I  (g'^y) )  *  |  -  2y 

2  2 

f„(y)  *  2ayexpl-ay  ]  y  >  0 


Let  a 


-  1/2< t2 

f^fy)  *  y/a2expt-y2/2<r2]  y  _>  0  (71A) 

Examination  of  equation  71A  shows  that  Y  is  a  Rayleigh  distributed 
random  variable  (11:195].  The  general  egression  for  the  moments  of  a 
Rayleigh  distributed  random  variable  are  given  by  ['.1:1481: 


E(Yn)  «  sqrt(  (pi )/2)[  1*3*5*.  .... nern]>  a  odd  (71B) 

*  2*^n(  .5n)  lorn ;  a  even  (71C) 

From  718  and  7 1C.  the  mean  and  variance  of  Y  are: 

E(YJ  -  sqrt((pi)/2)ff  (71D) 

var(Y)  -  E(Y2)  -  E(Y}2  v 

-  (2  -  (pi)/2)d2  (71E) 


The  parameter  a  of  the  Rayleigh  distributed  amplitude  will  be 
adjusted  in  the  fade  simulation  program  to  force  equation  710  to  be 
equal  to  the  deterministic  amplitude  previously  used  (A  *  5).  Note  that 
equation  7 IE  shows  that  the  variance  of  the  random  amplitude  will  change 
as  the  parameter  a  is  changed. 

Receiver  Structure 

Earlier  in  the  estimator  development,  the  necessary  condition  for 
the  maximum-likelihood  estimate  was  obtained  by  differentiating  the 
conditional  log-likelihood  function  with  respect  to  the  unknown 
parameter  and  setting  the  result  equal  to  zerd.  This  condition  was 
called  the  log-likelihood  equation  and  was  expressed  in  equation  25A. 

The  inherent  assumption  was  that  the  conditioning  due  to  the  amplitude 
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and  angle  of  arrival  resulted  from  deterministic,  but  unknown, 
parameters.  With  the  amplitude  now  known  to  be  a  Rayleigh  random 
variable,  the  log-likelihood  equation  must  be  expressed  in  terms  of  the 
a  priori  knowledge  of  the  amplitude  distribution.  To  account  for 
conditioning  with  a  random  variable,  an  average  log-likelihood  eq  ation 
can  be  used  to  remove  the  conditioning  by  averaging  over  all  possible 
amplitude  values  [16:310].  Averaging  equation  25A: 

E{0}  *  E(2r,V  sqrt(E)/N  -  2A(V  sqrt ( E ) ) 2/Nft  + 
is  US  0 

2r2VdSqrt(E)/N0  "  2A(Vdsqrt(E))2/M0) 

0  »  r,V  sqrt(E)  -  E{A)(V  )2E  +  r,V.sqrt(E)  -  E{A)(VJ2E 

13  S  2.  u  a 

A 

Solving  for  E(A  , } : 

A 

E(AnlJ  ■  r1Vgsqrt(E)  +  r2Vdsqrt(2) 

(v  )2e  +  (V.)2E 

S  a 

-  1  r,V  +  r.V  " 

sqrt(E)  - 5 — r - r  (72) 

L(VS)2  ♦  (Vd)2J 

Comparison  of  equation  72  and  the  previously  derived  estimate  for 
amplitude  (equation  26)  shows  that  on  the  average,  the  maximum- 
likelihood  estimate  for  amplitude  will  be  unaffected  by  its  random 
nature.  Recall  that  the  amplitude  estimate  is  later  utilized  to  develop 
the  monopulse  error  equation,  and  therefore  the  receiver  structure  that 
develops  an  error  signal  of  the  angle  off  boresight  is,  on  the  average, 
unaffected  by  the  Rayleigh  distributed  amplitude.  Equation  72  only 
considers  the  average  value  of  the  amplitude  estimate.  As  shown  earlier 
in  equation  7 IE,  the  random  amplitude  will  have  an  associated  variance. 
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This  amplitude  dispersion  will  have  a  direct  consequence  upon  receiver 
performance  for  a  fixed  radar  design  based  upon  deterministic 
amplitudes . 

Performance 

Rayleigh  Amplitude  Samples 

To  investigate  the  performance  aspects  of  fluctuating  targets,  it 
is  necessary  to  modify  the  simulation  program  of  Appendix  B  to  account 
for  "slow"  and  "ri  d"  fluctuations  on  signal  amplitude.  The  IMSL 
routine  GGWIB  provides  the  necessary  tool  for  developing  Rayleigh 
distributed  random  samples.  GGWIB  is  a  Weibull  random  deviate 
generator,  and  the  Rayleigh  density  is  a  special  case  of  the  more 
general  Weibull  density .  The  general  form  of  the  Weibull  probability 
density  function,  pdf,  of  the  random  variable  Y  is  given  by  [7:229]: 

fy(y)  -  (A/B) [ (y  -  C)/B] (A“1)exp[-(y-C)/B]Af  y  >  C  (73) 

where: 

A  *  shape  parameter;  A  >  0 
B  «  scale  parameter;  B  >  0 
C  *  location  parameter 

Consider  the  following  parameter  values  in  equation  73: 

A  «  2;  B  «.  sqrt(2)cr;  C  »  0 

f?(y)  »  'i/o' exp(-y2/2er2 ;  y  £  0  (73A; 

Comparison  of  equation  73A  with  equation  71A  shows  that  with  the 
parameter  values  as  designated,  the  Weibull  density  is  transformed  to 
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the  desired  Rayleigh  density.  This  transformation  is  used  in  the  target 


fluctuation  program  to  create  Rayleigh  distributed  amplitude  samples. 

An  annotated  program  listing  is  provided  in  Appendix  C. 

Signal-to-Noise  Ratio  (SNR) 

Signal-to-noise  ratio  as  designated  up  to  this  point  has  been  equal 
to  the  maximum  instantaneous  signal  power  to  the  average  noise  power. 
With  a  fluctuating  amplitude,  the  peak  value  of  signal  power  is  no 
longer  applicable,  and  an  average  signal-to-noise  ratio  must  be  defined. 
In  general: 

Egy  *  /s2(t,dt  =CJ"PW  A2dt  =  A2PW 

An  average  value  of  Egy  would  be: 

ECEgy}  «  E(A2}PW  *  2ct2PW  , 

where : 

from  71C,  E{A2}  -  2or2 

Egy  represents  signal  energy 

An  average  signal-to-noise  ratio  is  then: 

E{SNR)  «  2E(Egy}/NQ  «  4<T2PW /NQ 


Taking  the  ratio  of  peak  to  average  SNR: 

SNRpeak/SNRavg  “  E3y/E(Egy}  *  A2/E(A2) 

Solving  for  SNR  .  : 


“".v,  •  3NRi,.ikE'A2>/‘2 


(74) 


In  the  fluctuating  target  simulation,  the  standard  deviation  of  the 


Rayleigh  distribution  was  adjusted  to  provide  a  mean  value  of  amplitude 
equal  to  the  nonfluctuating  amplitude  (A  =  5).  This  was  done  to  provide 
a  performance  comparison  between  fluctuating  and  nonfluctuating  models. 
Solving  equation  7 ID  for  the  standard  deviation  required  for  a  mean 


The  result  of  equation  74A  indicates  that  to  maintain  the  same  mean 
value  of  amplitude  as  the  nonfluctuating  model,  1.27  times  as  much  SNR 
was  required  under  fading  conditions.  This  is  equivalent  to  adding  1.04 
dB  more  to  the  signal-to-noise  ratio  level.  Another  way  of  interpreting 
the  result  of  equation  74A  is  that  for  a  fixed  average  SNR,  the  average 
and  standard  deviation  of  the  amplitude  would  change  considerably  from 
the  nonfluctuating  values  of  A  **  5,  and  er  =  0.  Because  simulation  SNR 
was  free,  the  average  SNR  could  be  raised.  If  there  were  an  associated 

cost  with  increasing  average  SNR,  performance  may  degrade  even  further 

\ 

than  what  the  results  that  follow  indicate. 

fluctuating  Target  Simulation 

The  modifications  necessary  to  Include  provisions  for  amplitude 
fluctuations  are  simple  and  straight-forward.  The  program  modification 
is  anuoted  in  Appendix  C  and  includes  the  following  steps: 

I )  The  user  i*  prompted  for  "slow"  or  "rapid”  fading  of  the 
amplitude. 
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2)  The  external  Weibull  deviate  generator  is  called. 

3)  For  "rapid"  fading,  the  pulse  amplitude  is  fluctuated  for  each 


of  the  50  tracked  pulses.  For  "slow"  fading,  the  amplitude  is 
held  constant  for  each  run  (50  pulses)  but  is  allowed  to 
fluctuate  for  successive  runs.  Each  fluctuation  is  Rayleigh 
distributed  by  appropriately  transforming  the  Weibull  deviates. 

4)  The  tracking  error,  and  performance  bounds,  are  then  computed. 


The  outputs  of  the  fluctuating  target  simulation  program  are 
plotted  for  a  stationary  target  in  the  following  diagrams.  A  stationary 
target  was  chosen  because  as  noted  in  the  previous  chapter  the  only 
difference  between  stationary  and  moving  targets  was  a  slight  increase 
in  the  average  mean-square  error  (MSE)  value  and  larger  dispersions  in 
the  error  (increased  tracking  difficulty).  The  first  plot.  Figure  33, 
shows  the  general  trend  of  the  MSE  at  high  SNR  for  the  nonfluctuating 
model,  the  "slowly"  fluctuating  model,  and  the  "rapidly"  fluctuating 
model:  As  Figure  33  indicates,  even  at  high  SNR  the  rapidly  fluctuating 
model  is  very  erratic  and  with  a  considerably  larger  tracking  error  than 
either  the  nonf luctuating  or  "slowly"  fluctuating  models.  This  is  an 
indication  of  the  boresight  bouncing  wildly  around  the  target  position 
rather  than  settling  to  a  steady-state  MSE  value.  The  "slowly* 
fluctuating  model  shows  a  moderate  increase  in  MSE  above  the 


nonfluctuating  model. 


Slow  Fluctuation 

Figure  34  is  a  stationary  target  composite  plot  for  the  slowly 
fluctuating  target  mean-square  tracking  error  (MSE) ,  the  Cramer-Rao 
bound  (CEB),  and  the  Ziv-Zakai  bound  (ZZB)  for  SNR  levels  of  15,  7,  -7, 
and  -20  dB.  Comparison  of  Figure  34  with  the  previous  nonfluctuating 
target  model  shown  by  Figure  30  demonstrate  that: 

1)  The  MSE  for  "slowly"  fading  targets  shows  more  degradation  at 
high  SNR  than  for  low.  This  is'  an  indication  that  when  the 
target  tracking  environment  is  poor  (low  SNR)  slow  target  fades 
don't  contribute  significantly  to  the  tracking  error. 

2)  The  CRB  was  largely  affected  by  the  presense  of  target 
fluctuations. 

3)  Although  the  error  is  greater  for  all  levels  of  SNR,  the  radar 
is  still  able  to  track  the  target. 

4)  As  the  -20  dB  plot  indicates,  the  ZZB  is  still  tight  at  the 
lowest  measured  SNR  level  and  the  CRB  has  lost  its  usefulness. 


Rapid  Fluctuation 

Figure  35  provides  the  rapid  fluctuation  composite  plot. 

Comparison  of  Figure  35  with  the  previous  nonfluctuating  target  model 
shown  in  Figure  30  demonstrate  that: 

1)  Because  of  the  very  large  excursions  in  both  MSE  and  bounds,  it 
is  very  difficult  to  track  the  target  in  a  "rapid"  amplitude 
fluctuation  environment. 


2)  The  MSE  shows  relatively  little  sensitivity  to  changing  SNR  and 
doesn't  really  indicate  tracking  at  all. 

3)  The  CRB  exceeds  the  MSE  at  a  higher  SNR  level. 

4)  The  ZZB  shows  more  variation  at  low  SNR  than  previously  seen, 
but  still  proves  to  be  a  tight  bound  at  low  SNR  levels. 

5)  Characteristic  of  the  "rapid"  fluctuating  model  are  deep  nulls 
and  peaks  from  pulse  to  pulse  in  the  MSE  and  CRB  curves.  This 
indicates  that  tracking,  which  should  be  based  upon  past 
history,  is  independent  from  pulse  to  pulse.  Clearly  this  is 
not  a  very  desireable  tracking  characteristic. 
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Figure  35.  Rapidly  Fluctuating  Target  Tracking  Error  vs  SNR  Level 


Smoothing  Effects 

In  an  attempt  to  smooth  out  amplitude  fluctuations,  the  fluctuating 
target  simulation  of  Appendix  C  was  modified  to  accumulate  the  error 
signal  on  N  pulses  before  moving  the  boresight.  The  annotated  smoothing 
program  is  provided  in  Appendix  0  and  includes  the  following  steps: 

1 )  The  user  is  prompted  for  the  number  of  pulses  to  be  observed 
before  boresight  movement  (N).  N  must  be  a  factor  of  50  to  fit 
into  the  previous  programming  structures. 

2}  The  error  for  the  designated  N  pulses  is  then  accumulated  and 
averaged.  This  averaging  injects  the  "smoothing"  of  amplitude 
fluctuations  over  N  pulses. 

3)  The  boresight  is  moved  at  the  end  of  the  Nth  pulse  based  upon 
the  average  measured  error. 

4 )  Tracking  error  and  bounds  are  then  computed  in  the  same  fashion 
as  the  previous  programs. 


Figure  36  provides  a  plot  of  the  MSE  at  high  SNR  for  a  stationary 
target  with  "rapid"  amplitude  fluctuations  and  pulse  smoothing.  Two  and 
five  pulses  were  observed  before  boresight  movement.  Also  included  in 
Figure  36  is  the  nonfluctuating  MSE.  Inspection  of  the  provided  plot 
clearly  shows  the  tendency  of  smoothing  to  reduce  the  tracking  error. 
Further,  as  the  number  of  pulses  increases,  the  MSE  approaches  the 
nonfluctuating  model  MSE.  It  is  worth  noting  that  this  increase  in 
tracking  reliability  in  a  fading  environment  comes  at  the  expense  of 
requiring  N  pulses  to  accurately  locate  target  pc  ion. 
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V.  MULTIPATH 


The  last  chapter  was  an  extension  to  the  initial  model  that 
considered  target  tracking  and  performance  estimates  in  the  presence  of 
a  fluctuating  target.  This  chapter  extends  the  initial  model  to 
determine  performance  estimates  and  target  tracking  capability  for  a 
more  complex  channel  model,  specifically,  tracking  in  the  presence  of  a 
multipath  environment. 

As  a  signal  is  transmitted  from  one  point  to  another  ~^ove  an 
irregular  surface  with  electrical  properties  differing  from  that  of  the 
propagation  medium,  reflected  signals  may  be  generated.  These  signals 
then  appear,  along  with  the  direct  path  signal,  at  the  antenna  of  the 
radar.  The  reflected  signals  are  commonly  called  multipath.  The 
concern  of  this  section  is  to  determine  the  effects  of  multipath  on 
azimuth  estimation  performance. 

Rather  than  a  complete  spatial  diffraction  description  of  specific 

environmental  structures,  the  attempt  will  be  to  develop  an  adequate 

model  to  characterize  and  simulate  the  multipath  channel.  One  specific 

type  of  multipath  will  be  focused  upon,  terrain  bounce.  Figure  37 

describes  the  geometry  of  the  terrain  bounce  environment.  Observe  from 

Figure  37  that  the  target  radiates  or  reflects  signals  in  all  directions 

and  the  radar . receives  a  direct  signal  from  the  angle  8^ ,  and  a 

t 

reflected  signal  from  the  angle  9  . 


Figure  37.  Terrain  Bounce  Environment: 

Source  [16:688] 

The  reflected  path,  or  multipath  signal  return,  will  depend  upon 
the  relection  characteristics  of  the  terrain.  Strong,  well  defined 
returns  are  called  specular  multipath  and  they  result  from  large  flat 
reflecting  surfaces.  A  second  kind  of  multipath,  diffuse  multipath, 
occurs  due  to  small  surface  irregularities  and  appear  to  come  from  all 
angles  instead  of  a  single  well  defined  direction.  In  this  section, 
both  specular  and  diffuse  multipath  effects  on  azimuth  estimation 
performance  are  of  interest.  Figure  38  describes  the  two  types  of 
multipath. 
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Figure  38.  Multipath  Types 


Tracking  in  a  multipath  environment  will  be  investigated  in  the 
following  steps t 

1)  A  model  will  be  generated  for  the  return  signal  in  the 
multipath  environment . 

2)  Azimuth  estimation  performance  will  be  simulated  in  the 
presence  of  specular  or  diffuse  multipath. 

3)  A  method  to  reduce  the  tracking  error  for  target  tracking  in 
the  presence  of  multipath  will  be  proposed  and  investigated. 

Signal  Representation 

Characteristic  of  the  multipath  medium  is  the  time  spread 
introduced  in  the  signal  which  is  transmitted  through  the  channel 
[17:455].  A  short  transmitted  pulse  over  a  multipath  channel  might  be 
received  as  a  train  of  pulses.  Figure  38  shows  a  possible  response  to  a 
transmitted  pulse  through  the  multipath  channel. 
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Figure  39.  Possible  Response  to  a  Transmitted  Pulse 
Source:  [17:455] 


Another  characteristic  inherent  to  the  radar  problem  vs  that  the 
received  direction  of  the  individual  reflected  signals  will  differ  from 
each  other  and  from  the  direct  signal  received  from  the  true  target 
location. 

The  objective  of  the  radar  signal  model,  as  developed  here,  is  to 
obtain  a  received  observation  that  is  a  function  of  both  time  and 
azimuthal  direction  that  accounts  for  the  following  effects: 

1)  Direct  returns  from  the  target. 

2)  Reflected  returns  (either  specular  or  diffuse)  from  the 
terrain. 

3)  The  possibility  of  the  superposition  of  returns  (pulse  overlap) 
that  may  add  in  a  constructive  or  destructive  manner. 
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General  Signal 


Most  narrowband  radar  signals  can  be  expressed  as  [18:298]: 

s { t )  -  A(t) cos [wQ ( t)  +  0<t)3  (75) 

where : 

wQ  =  carrier  frequency  (rad/sec) 

A(t)  *  instantaneous  amplitude 

0 ( t)  *  instantaneous  phase 

A(t)  &  0(t)  vary  slowly  compared  to  wQ 

To  characterize  both  direct  pulses  and  a  random  number  of  reflected 
pulses  with  both  specular  and  diffuse  multipath  environments,  equation 
75  must  be  extended  to  a  more  gerrral  descriptive  form  of  s(t).  For  any 
arbitrary  pulse,  a  general  representation  of  the  signal  received  at  the 
radar  can  be  expressed  as: 

s(t)  *  Re(u(t)exp[ jWg(t) ) }  (75A) 

where: 

u(t)  is  the  complex  envelope  and  can  be  expressed  as 
u(t)  »  TA<t)exp[j0(t)J 
r  *  complex  reflection  coefficient 
Aft)  «  discrete  random  process  that  describes  the  received 
pulses  in  both  quantity  and  time  of  occurance 
0 ( t )  *  arbitrary  phase  relative  to  some  reference 
complex  representation  is  used  to  allow  combining  signals 
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Complex  Reflection  Coefficient 

The  reflection  coefficient  of  the  terrain  surface  is  generally  a 
complex  quantity,  T  =  pexp[ j^l .  The  real  part  of  r  (p)  describes  the 
amplitude  changes  on  reflection,  while  the  imaginary  part  of  r  (^) 
describes  the  phase  changes  on  reflection.  To  assign  appropriate 
magnitude  coefficients  for  the  different  multipath  cases  (spectilar  or 
diffuse),  an  average  value  of  surface  roughness  was  employed.  For 
average  surface  roughness,  the  magnitude  (p)  of  the  complex  reflection 
coefficients  for  both  specular  and  diffuse  multipath  is  approximately 
0.4  [16:689].  To  assign  the  appropriate  phase  after  reflection,  a  phase 
angle . of  pi  radians  applies  to  a  smooth  surface  with  good  reflecting 
properties  if  the  radiation  is  horizontally  (tangential  electric  fields 
across  a  boundary)  polarized  and  the  angle  of  incidence  is  small 
(specular)  [15:444],  and  diffuse  reflection  is  assumed  uniform  in  phase 
from  0  to  2pi.  The  complex  reflection  coefficients  are: 

specular:  T  *  0.4/pi  (76) 

diffuse:,  T  *  0.4/^>  (77) 

where: 

*  is  U(0,2pi) 

Envelope  A(t) 

The  changing  terrain  will  cause  independent  path  lengths  and  a 
random  number  of  pulse  returns,  with  unique  arrival  times,  at  the  radar 
receiving  antenna.  Defining  a  randan  variable  A  such  that  A(f)  equals 
the  total  number  of  returns  within  the  radar  pulse  repetition  interval 
(PRZ),  A  only  takes  on  discrete  integer  values  0,  1;  2,  ...  and  the 


event  {A=k}  is  the  event  k  returns  occur  in  the  time  interval  [0,PRI>. 
The  probability  of  k  returns  within  the  PRI  is  given  by  [11:285]: 

P{A  =  k)  -  exp[-\PRI]  (XPRDVkJ  k  *  0,  1,  ...  (78) 

where: 

X  =  average  number  of  occurances  per  unit  time 

Let  a  *  XPRX 

»  rate'time  *  avg  occurances  in  the  time  interval  (O.PRI) 

P{A  «  k)  «  exp[-a] (a)*/kl  k  «  0,  1,  ...  (78A) 

Equation  78A  gives  the  probabiltiy  mass  function  of  A(t),  the 
number  of  pulse  arrivals  during  the  PRI.  It  describes  the  Poisson 
distribution  with  parameter  a  [7:171].  At  this  point,  it  is  necessary 
to  determine  the  placement  of  c&ch  returned  pulse  on  the  time  axis. 

Each  time  interval  between  successive  pulses  will  be  a  random  variable 
and  can  be  termed  as  the  interarrival  time  for  the  Poisson  distribution. 
To  find  the  distribution  for  the  interarrival  times  assume: 

1 )  all  events  are  equiprobable 

2)  the  first  event  (To)  must  occur  after  some  time  t 
(ie  A(t)  *  0  for  t  <  To) 

P(A  -  0)  «  P(To  >  t)  «  exp[-a] (a)°/0l, 

P(To  >  t)  •  exp [-a] 

P{To  <  t)  •  P(To  <.  t)  «  1  -  expl-a]  (79) 

The  probability  distribution  of  the  interarrival  times  is  expressed 
in  equation  79,  and  its  corresponding  density  function  is  that  of  the 
Exponential  distribution  [7:209].  The  Exponential  density  function  was 


previously  defined  by  equation  69.  Thus,  the  envelope  of  the  general 
representation  of  the  signal  received  at  the  radar  will  take  the  form  of 
a  Poisson  distribution  with  Exponential  interarrival  times. 

Instantaneous  Phase 

To  account  for  the  phase  of  the  complex  envelope,  all  instantaneous 
phases  are  assumed  as  time  independent  and  equiprobable  for  any  phase 
between  0  and  2pi.  This  establishes  phases  as  being  uniformly 
distributed  from  0  to  2pi  and  constant.  From  equation  75A  and  the 
preceding  development,  the  received  signal  at  the  radar  vill  take  on 
three  seperate  forms: 

1)  For  the  direct  pulse:  +  -  0,  0(t)  is  U(0,2pi),  p  *  l 

s(t).  *  A(t)cos[wQ(t)  +  0]  (80) 

2)  For  specular  reflections;  ^  *  pi,  0(t)  is  U(0,2pi),  p  «  0.4 

s(t)  ■  0 . 4 A ( t ) cos [ Wg ( t )  +  0  +  pij  ( 80A) 

3)  For  diffuse  reflections;  ^*U(0,2pi),  0(t)  is  O(0,2pi),  p  -0.4 

s(t)  «  0.4A(t)cos(wQ(t)  +  0  +  <f>]  (80B) 

Observations  , 

The  received  observations  will  in  general  bet 

r(t,9)  »  s(t,8)V(0)  +  N(t) 

where: 

s(t,0)  is  described  by  the  appropriate  form  of  equation  80 

N(t)  is  additive  white  Gaussian  noise 

V( 0 )  is  an  antenna  weighting  dud  to  the  angle  of  arrival 


The  model  for  the  antenna  weighting,  V(8),  will  take  on  one  of 
three  seperate  forms.  The  direct  pulse  will  have  an  angle  of  arrival  as 
described  by  Figure  37.  The  two  forms  of  reflection  will  have  angle  of 
arrivals  as  described  by  Figure  38.  Therefore! 

1)  For  the  direct  pulse;  0  *  target  angle  of  arrival 

2)  For  specular  reflections;  with  respect  to  the  boresight,  0 
will  be  uniformly  distributed  from  0  «  0  to  8  «  BW/2.  All 
reflected  pulses  will  be  allowed  to  form  observations. 

3)  For  diffuse  reflections;  0  will  be  uniformly  distributed  from 
6  *  0  to  8  *  pi.  Only  reflected  pulses  that  fall  within  the 
interval  (0,BW/2)  will  be  allowed  to  form  observations,  as  all 
other  returns  will  not  get  into  the  received  antenna J 

Pulse  Overlap 

With  Exponential  interarrival  times,  provisions  must  be  made  for 
the  possibility  of  two  pulses  arriving  into  the  receiving  antenna  during 
the  same  time  interval.  Such  an  occurance  would  be  pulse  overlap,  and 
for  simplification  this  investigation  considers  any  ovelap  to  be 
complete  overlap  (ie  for  any  portion  of  overlap  in  the  pulsewidth, 
overlap  is  considered  to  be  of  length  equal  to  the  pulsewidth)  with  a 
resultant  received  pulse  equal  to  the  superposition  of  the  amplitudes 
and  phases  of  the  two  affected  pulses.  This  assumption  is  reasonably 
valid  provided  no  doppler  shifts  of  either  of  the  two  signals  occurs  (ie 
the  sum  of  two  sinusoids  at  the  same  frequency  is  sinusoidal  at  the  same 
frequency). 

recall  u(t)  *  A(t)exp( j0(t)  1 
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in  polar  form:  u(t)  =  Um/gm 


(81) 


where: 

Um  =  |u(t)|,  =  |r|A(t)  =•  pA(  t) 

0m  =  /u(t)  =  0  +  ^ 

in  rectangular  form:  u(t)  =  X  +  jY  (81A) 

where : 

X  *  Umcos[0m]  »  Re{u(t)} 

Y  *  Umsin(gm)  *  Im(u(t)} 

Using  the  rectangular  form  of  u(t)  as  described  by  equation  81A, 
two  overlapping  pulses  can  be  combined  by  adding  the  real  and  imaginary 
components  of  each  pulse.  For  the  two  overlapping  pulses  at  the  same 
frequency: 

Re{u^  ( t) }  *  Um^ costing] ;  ImtCu^t}}  *  Um^sintim^ 

Re(u2(t)j  *  Um2cos[0m2];  la{{u2(t)}  =*  Um2sin{gm2] 
where: 

Um  and  0m  are  described  by  equation  81 

The  resultant  magnitude  and  phase  of  the  sum  of  the  two  pulses  are: 

|u(sum)j  -  sqrt t (Um^cos( 0m  ^ )  +  Ura2 jos( 0o2 ) )2  ♦ 

(Um1sin(gm1)  +  Um2sin(0m2) )2]  (82) 

/g(sum)  «  [Urn  sin(0a  )  +  um  sin(0ra  )' 

tan*  . .  2 . .  (83) 

Um,)coB(0a1 )  +  Um2cos(0ra2)_ 

Only  the  magnitude  (equation  32)  is  of  interest  in  the  combined 
resultant  signal.  Equation  82  will  be  used  in  the  multipath  simulation 
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program  to  combine  two  pulses,  the  resultant  signal  {s^ft}]  is: 

sr(t)  *  Re( lu(sum) |/0(sum)exp[ jwQ(t)]  ) 

*  |u(sum) jcos(w0(t)  +  D(sum)] 

Performance 

To  investigate  the  performance  aspects  of  target  tracking  in  the 
presence  of  multipath,  extensive  modifications  are  required  in  the 
simulation  program  of  Appendix  B.  Extension  to  a  multipath  environment 
requires  a  random  number  of  additional  pulses  that  will  form 
observations  during  each  pulse  repetition  interval  ( PRI ) .  An  annotated 
listing  of  the  simulation  program  in  the  presence  of  multipath  is 
provided  in  Appendix  E.  In  the  multipath  simulation,  each  received 
pulse  is  allowed  to  form  an  observation,  compute  an  error  estimate,  and 
move  the  boresight  in  a  direction  to  reduce  the  error  estimate.  To 
facilitate  comparison  with  previous  results  the  simulation  outputs,  the 
tracking  error  and  performance  bounds,  are  again  computed  after  each 
pulse  repetition  interval.  Both  diffuse  and  specular  multipath 

1  i 

environments  are  considered. 

Poisson,  Distributed  Returns 

As  stated  in  equation  78A,  the  model  used  to  simulate  a  random 
number  of  arrivals  during  the  PRI  is  a  discrete  Poisson  distribution  on 
the  envelope  of  the  signal,  A(t).  The  IMSL  routine  GGPOS  provides  the 
necest ary  tool  for  developing  the  Poisson  deviates.  The  Poisson 
parameter,  a  of  equation  78A,  is  an  input  parameter,  RLAM,  of  the 
simulation  program  of  Appendix  E,  and  was  adjusted  during  performance 
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verification  to  simulate  different  multipath  environments. 

Specifically,  RLAM  values  of  0.5,  1.0,  and  15,  respectively,  were  used 
to  model  increasing  multipath  scenarios.  The  GGPOS  output  vector,  IR, 
contains  the  50  Poisson  deviates  describing  the  number  of  random 
reflections  for  each  PRI  of  the  50  pulse  run. 

Exponential  Interarrival  Times 

For  each  of  the  50  direct  pulses  of  a  run,  the  Poisson  generated 
number,  1R( Pulse)  of  the  simulation,  describes  the  number  of  reflected 
pulses.  As  shown  in  equation  79,  the  time  of  arrival  for  each  of  these 
pulses  will  be  Exponentially  distributed.  The  XMSX,  routine  GGEXP 
provides  the  necessary  tool  for  developing  pulse  interarrival  times 
within  the  PRI.  The  number  of  required  Exponential  interarrivals,  NEX 
of  the  simulation,  equals  the  random  GGPOS  output  for  a  specific  PRZ. 
The  generated  random  arrival  times,  the  vector  E  of  the  simulation,  are 
sorted  to  provide  sequential  arrivals  in  time.  The  difference  between 
arrival  times,  DTIME  of  the  simulation,  establishes  the  response  time 
that  the  servo-loop  will  have  to  each  generated  observation. 

Pulse  Overlap 

As  mentioned  earlier,  pulses  are  considered  to  be  totally 
overlapped  if  any  portion  of  the  pulses  overlap  ( ie  if  the  time 
difference  between  pulses  is  less  than  twice  the  pulse  width) .  Both 
direct  and  reflected  pulses  are  considered  for  the  possibility  of 
overlap,  if  overlap  occurs,  the  resultant  pulse  magnitude  is  computed 
by  equation  82 .  Magnitudes  are  assigned  according  to  the  appropriate 
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reflection  coefficient  times  the  pulse  amplitude  for  each  considered 
pulse.  All  phases  are  assigned  a  random,  uniform,  number  between  0  and 
2pi.  The  IMSL  routine  GGUBS  provides  the  necessary  tool  for  developing 
uniform  deviates  from  0  to  1,  and  an  appropriate  scaling  of  the 
generated  vector  ot  uniform  deviates  yields  the  desired  random  phases. 
Upon  computation  of  the  resultant  magnitude  and  phase  of  two  overlapping 
pulses,  a  program  flag  is  set  which  disregards  the  next  incoming  pulse 
and  subsequentialy  reset  in  anticipation  of  another  overlap. 

Angle  of  Arrival 

The  1KSL  routine  GGUBS  also  provides  provisions  for  random  angle  of 
arrivals  (AOA) .  For  diffuie  multipath,  IR( PULSE)  (Poisson  number  of 
generated  reflections)  uniform  deviates  are  generated  by  an  external 
call  to  GGUBS  and  scaled  to  be  uniform  from  0  to  (pi).  Each  reflected 
pulse  is  then  assigned  a  unique  AOA.  likewise,  specular  multipath 
uniquely  assigns  IR( PULSE)  AOA's  but  the  deviates  are  scaled  to  be 
uniform  from  0  to  BW/2.  No  AOA  is  assigned  to  the  direct  pulse,  as  it 
locates  the  true  target  location.  In  the  simulation  program,  no 
observations  are  formed  for  AOA's  that  are  greats,  than  BW/2  (diffuse 
multipath) .  Therefore,  the  difference  in  the  simulation  program  between 
specular  end  diffuse  multipath  is  the  AOA's  that  are  allowed  to  form 
observations  (all  for  specular;  only  those  between  0  and  BW/2  for 
diffuse). 

The  outputs  of  the  multipath  simulation  program  are  plotted  for  a 
stationary  target  in  the  following  diagrams. 
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Diffuse  Multipath 


Figure  40  is  a  composite  plot  showing  diffuse  multipath  effects  on 
the  mean  square  tracking  error  for  a  stationary  target,  at  high  SNR, 
with  varying  degrees  of  multipath  ( parang t^r  lambda) .  The  upper  left 
corner  shows  the  MSE  for  the  direct  pulse  with  no  multipath,  and  is 
provided  for  comparison.  Comparison  of  the  plots  with  parameter  lambda 
unequal  to  zero  with  the  direct  pulse  only  plot  demonstrate  that: 

1)  In  all  cases,  the  MSE  starts  at  some  large  value  and  settles  to 
a  smaller  value.  This  indicates  good  tracking  capability. 

2)  Although  the  actual  MSE  changes  from  plot  to  plot  vMSE  is  a 
random  variable),  diffuse  multipath  does  not  significantly 
degrade  tracking  capability.  This  result  verifies  a  statement 
made  in  the  literature  [19:1121  where  diffuse  multipath  is 
claimed  to  present  no  serious  limitation  in  direction  finding 


capabilities 


Figure  40.  Diffuse  Multipath  Tracking  Error 


Specular  Multipath 

Figure  41  is  the  corresponding  composite  plot  for  specular 
multipath  effects  on  the  MSE  for  a  stationary  target,  at  high  SNR,  with 
varying  degrees  of  multipath  (parameter  lambda).  Comparison  of  the 
three  multipath  plots  Wxth  the  direct  pulse  only  plot  demonstrate  that: 

1)  Serious  monopulse  tracking  errors  occur  in  the  presence  of 
specular  multipath. 

2)  At  an  average  one-half  additional  pulse  per  PRI,  MSE  has 
increased  considerably  but  still  shows  tracking  capability. 

3)  At  an  average  of  one  extra  random  pulse  per  PRI,  the  monopulse 
radar  tracks  for  approximately  5  pulses  and  then  the  tracking 
error  starts  to  run  away.  This  is  an  indication  of  not  enough 
SNR  to  perform  the  tracking  function  (breaklock) .  Observe  that 
there  are  a  few  points  where  the  error  dips  down,  or  decreases, 
in  the  lambda  “  1  plot.  The  bores ight  has  in  these  few 
instances  moved  to  reduce  the  error  and  for  the  purposes  of 
this  simulation  lambda  «  1  is  the  tracking  threshold  for  the 
radar. 

4)  At  an  average  of  one  and  one-half  additional  pulses  per  PRI, 
the  tracking  error  increases  monotonically.  Clearly,  the 
monopulse  tracker  has  lost  its  usefulness.  The  monotonic 
increase  in  tracking  error  may  be  useful  though  as  a  multipath 


detector  (19:61). 


I 


a 

l 

ii 

4 

1 

*  • 

i  W 


Specular  Multipath  vs  SNR 

In  general,  tracking  capability  was  indicated  for  specular 
multipath  only  when  the  parameter  lambda  (average  number  of  returns  in 

the  PRX)  was  less  than  1.0.  Of  concern  in  this  section  is  how  tracking 

\ 

and  performance  bounds  are  affected  by  changing  SNR  levels.  Specular 
multipath  was  chosen  for  observation  because  the  previous  results 
indicate  that  no  serious  tracking  limitations  are  introduced  by  diffuse 
multipath  enviornments .  Also,  since  tracking  capability  was  indicated 
when  the  parameter  lambda  was  less  than  1.0,  the  Poisson  paramenter 
lambda  whil>  'onitoring  tracking  potential  will  be  established  as  lambda 
*  0.5.  Figtu,  42  is  a  stationary  target  composite  plot  for  MSE,  the 
CRB,  and  the  ZZB  for  SNR  levels  of  15,  7,  -7,  -20  dB  in  the  presence  of 
specular  multipath.  Comparison  of  Figure  42  with  the  previous  direct 
pulse  only  target  model  shown  by  Figure  30  demonstrate  that: 

1)  Even  for  high  SNR  levels,  the  MSE  does  not  settle  to  a  steady 
state  value  but  is  larger  at  the  end  of  the  run  (50  pulse)  thah 
an  average  value  of  MSE.  This  is  indication  of  a  general 
tendency  of  the  MSE  to  be  increasing  in  the  multipath 
environment. 

2)  At  -7dB  SNR  level,  although  tracking  is  indicated  the  monopulse 
tracker  is  clearly  in  a  threshold  situation  where  there  is  not 
an  adequate  amount  of  signal  power  to  determine  target 
direction. 

3)  Once  again  the  CRB  exceeds  the  MSE  at  low  SNR  levels  while  in 
contrast  the  ZZB  lower  bounds  the  MSE,  even  at  -20dB. 
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In  an  attempt  to  overcome  the  tracking  difficulties  introduced  by  a 
multipath  environment,  the  multipath  simulation  program  of  Appendix  E 
was  modified  to  "lock-out"  all  undesired  pulses  within  the  radar's  PRI. 
The  "lock-out"  technique  requires  that  the  first  4  PRI's  be  tracked  in  a 
multipath  environment,  thus  the  radar  must  be  able  to  sacrifice  accurate 
tracking  in  the  first  4  PRI's  in  favor  of  increased  tracking- capability 
later  in  the  run  (the  remaining  46  pulses).  Additionally,  only  one 
target  is  assumed  fo  be  present.  Discarding  returns  (as  is  done  in  this 
technique)  with  multiple  targets  present  may  affect  target 
resolvability. 

The  objective  of  the  simulation  modification  is  to  look  at  the 
target  returns  over  the  first  4  PRI's  and  extract  some  form  of 
periodicity  in  the  target  returns.  The  characteristic  focused  upon  in 
this  investigation  was  the  time  of  arrival  of  pulse  returns.  Recall 
that  the  interarrival  times  are  random,  Exponentially  distributed,  over 
the  PRI.  The  probability  that  pulse  returns  would  consecutively  occur 
at  the  same  instant  in  time  for  the  first  4  PRI's  and  that  the  pulse 
would  be  due  to  a  random  terrain  reflection  is  quite  small  (considered 
equal  to  0  in  this  simulation).  Using  this  reasoning,  the  pulse  arrival 
time  (arrival  time  »  0  for  the  direct  pulse,  and  repetitive)  within  the 
first  4  PRI's  can  be  used  to  "lock-in"  the  desired  target  return  and 
"lock-out"  all  other  returns  not  oceuring  at  the  proper  time  instants. 
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The  annotated  "lock-out"  program  listing  is  provided  in  Appendix  F  and 
includes  the  following  steps: 

1)  A  (4X25)  arrival  time  storage  matrix  is  initialized  to  a  value 
of  1.0.  The  4  rows  represent  the  first  4  PRI's,  and  the  25 
columns  represent  the  direct  pulse  and  a  possible  maximum  of  24 
reflected  pulses  within  each  PRI. 

2}  As  the  simulation  program  proceeds,  the  arrival  time  for  each 
pulse  (direct  and  reflected)  is  stored  in  the  proper  matrix 
element  position. 

3)  At  the  beginning  of  the  5th  PR2,  each  of  the  25  columns  are 
averaged  over  the  4  PRI's.  This  column  average  is  then 
compared  to  each  corresponding  column  element.  If  total 
equality  for  all  elements  is  obtained  with  the  average  and  the 
average  is  not  equal  to  the  initialized  value  (1.0),  or  the 
next  PRI,  then  the  equality  designates  the  desired  pulse 
arrival  time  for  "lock-in".  All  other  received  pulses  will  be 
disregarded  in  the  5th  and  subsequent.  PRI's  except  the 
designated  "lock"  pulse. 

Figure  43  provides  a  plot  of  the  MSE  at  high  SNR  for  a  stationary 
target  in  a  specular  multipath  environment  and  pulse  "lock-out" 
employed.  For  comparison  with  the  previous  performance  curves  of  Figure 
41,  the  same  values  of  lambda  (Poisson  parameter)  were  used.  Comparison 
of  Figure  43  and  Figure  41  show  that  the  "lock-out"  technique  has 
removed  the  tracking  degradation  introduced,  by  all  of  the  multipath 
environments.  All  of  the  tracking  errors  settle  in  to  a  final  value 
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that  is  representative  of  stationary  target  tracking  at  high  SNR 
(reference  upper  left  corner  plot  of  Figure  43).  Zt  is  well  worth 
noting  that  this  increase  in  tracking  reliability  in  a  specular 
multipath  environment  is  for  the  particular  model  of  multipath,  and  form 
of  radar  simulation,  that  has  been  defined  herein  and  may  require 
extensions  to  combat  an  actual  multipath  situation. 
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Figure  43.  "Lock-out"  Effects  Upon  Multipath  Tracking 


VI .  RESULTS  AND  RECOMMENDATIONS 


This  chapter  summarizes  the  .devious  results  for  the  various 
tracking  capabilities  and  the  applicability  of  the  performance  bounds 
for  the  amplitude-comparison  monopulse  radar.  In  addition,  a  section  is 
provided  for  recommendations  for  continued  study  along  the  lines 
followed  by  this  thesis. 


Stationary  Target  Tracking 

As  Figure  29  indicates,  the  general  trend  of  the  mean-square 
tracking  error  (MSE)  for  a  stationary  target  is  inversely  proportional 
to  the  available  signal-to-noise  ratio  (SNR)  level.  There  is  also  a 
general  tendency  for  the  MSE  to  be  spread  over  a  larger  range  of  values 
for  the  smaller  values  of  SNR. 


Moving  Target  Tracking 

Figure  31  indicated  two  significant  results  for  the  MSE  for  a 
moving  target  relative  to  the  previous  stationary  target  tracking 

I 

trends:  I 


1)  For  high  SNR,  the  moving  target  MSE  settles  out  to  the  same 
value  as  for  stationary  target  tracking. 

2)  For  low  SNR  levels,  the  moving  target  MSE  is  slightly  higher 
than  the  corresponding  stationary  values,  and  the  MSE  spreading 
is  noticeably  larger  (increased  tracking  difficulty). 


Cramer-Rao  Bound 


Present  in  all  of  the  obtained  results  is  the  inability  of  the 
Cramer-Rao  bound  (CRB)  to  predict  tracking  capability  in  low  SNR 
environments.  Performance  results  indicate  that  the  applicability  of 
the  CRB  as  a  performance  indicator  is  limited  to  prior  knowledge  that 
the  MSE  exceeds  the  bound  (ie  the  SNR  is  above  a  threshold  situation). 
Further,  as  Figure  35  indicates,  the  CRB  is  extremely  susceptable  to 
target  amplitude  fluctuations  (producing  a  much  higher  error  with 
correspondingly  larger  valleys  and  peaks  in  a  fluctuating  target 
environment) . 

Ziv-Zakai  Bound 

The  Ziv-Zakai  bound  (ZZB)  in  all  observed  tracking  environments  was 
approximately  constant  and  no  greater  than  the  average  MSE  excursions. 
For  high  SNR  levels;  a  large  magnitude  difference  exists  between  the  ZZB 
and  MSE  (See  Figure  30).  Thus,  the  ZZB  is  more  applicable  at  lower  SNR 
levels.  For  low  SNR  levels,  the  ZZB  proved  to  be  very  tight  to  actual 
performance  values  in  that  the  bound  always  existed  at,  or  slightly 
below,  the  arithmetic  mean  of  the  MSE  excursions  (see  the  -20dB  »NR 
level  plots  of  Figures  30,  32,  34,  35,  and  42).  Thus,  the  ZZB  w.is  found 
to  be  a  very  good  performance  indicator  for  low  SNR  levels. 


Target  Fluctuations 


Significant  results  of  tracking  in  a  fluctuating  target  amplitude 
environment  are: 

1)  The  maximum-likelihood  receiver  structure  was  unaffected  by  the 
nature  of  the  fluctuating  target  amplitude. 

2)  "Slow"  amplitude  dispersions  only  moderately  deg. Aded  tracking 
capability  (see  Figure  34). 

3)  "Rapid"  amplitude  dispersions  imposed  serious  tracking 
limitations  (tracking  degradation)  upon  a  fixed  radar  design 
based  upon  deterministic  amplitudes  (see  Figure  35). 

Multipath 

Significant  results  of  tracking  in  a  multipath  environment  arc: 

1)  Diffuse  multipath  presents  no  serious  limitations  in  target 
tracking  capabilities  (see  Figure  40). 

2)  Serious  tracking  errors  occur  in  the  presence  of  specular 
multipath  (see  Figures  41  and  42)'. 

Recommendations 

In  the  way  of  recommendations  for  continued  study  along  the  lines 
established  by  this  thesis,  the  following  studies  art  proposed: 

Simultaneous  Parameter  Estimation 

One  of  the  possible  extensions  of  the  work  presented  here  would  be 
to  include  another  radar  parameter  in  the  estimation  process  and 
investigate  the  multiple  parameter  ambiguity  fuction.  Possible 
candidates  for  the  study  might  be  the  signal  delay  or  doppler. 


Two-Dimensional  Radar 


Appendix  A  provides  the  work  done  in  this  thesis  on  a  two- 
dimensional  target  tracking  radar.  The  mapping  from  probability  space 
to  a  moment  space  (in  the  estimation  error)  was  accomplished  via  a  two- 
dimensional  generalization  of  the  Tchebycheff  inequality.  The  method 
utilized  to  cone  up  with  a  two-dimensional  mean-square  estimation  error 
requires  close  scrutiny  as  to  its  correctness  and  applicability. 

Further,  the  error  probabilities  for  both  estimation  and  detection  as 
calculated  in  this  development  require  exhaustive  integrations  of 
standard  normal  density  functions  at  low  SNR  levels.  The  efficiency  of 
the  computational  algorithm  could  be  considerably  improved  with  a  Q 
function  "look-up”  table  rather  than  using  computer  integrations.  Such 
an  improvement  would  allow  runs  at  lower  SNR  levels  than  were  possible 
in  this  investigation  and  would  directly  provide  indication  about  the 
applicability  of  this  form  of  the  ZZB  to  the  two-dimensional  tracker. 

Simultaneous  Multipath  Components 

This  thesis  focused  in  on  either  specular  or  diffuse  multipath 
environments  when  in  a  practical  tracking  environment  both  components 
would  be  present  simultaneously.  Zt  is  well  worth  recommending  an 
investigation  that  considers  both  multipath  components  present  and  their 
resulting  consequences  upon  receiver  performance.  The  idea  would  be  to 
establish  a  probabilistic  description  of  the  received  envelope  by 
allowing  a  constant  specular  component  added  to  a  Gaussian  distributed 
diffuse  component  (Central  Limit  Theorem  applied  to  many  diffuse 
scatterers*.  The  constant  plus  a  Gaussian  should  result  in  a  Rician 
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envelope  distribution.  It  would  be  interesting  to  monitor  receiver 
performance  while  varying  the  specular  component  of  multipath  [ie 
varying  the  specular  component  from  zero  (diffuse  component  only),  to  an 
intermediate  value  (combined  specular  and  diffuse  components),  to  a 
large  value  ( specultir  component  only )  ]  . 


VS 
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Appendix  A 

TWO-DIMENSIONAL  RADAR  MODEL 

Introduction 

In  its  simplist  form,  amplitude-comparison  monopulse  radar 
determines  target  direction  by  comparing  a  single  pair  of  signals 
received  on  two  antenna  patterns  simultaneously.  This  is  sufficient  to 
determine  the  target  angle  of  arrival  in  a  single  plane  (azimuth  or 
elevation).  This  simple  form  of  monopulse  is  the  radar  model  thus  far 
developed.  Three  dimensional  tracking,  however,  requires  comparison  of 
two  pairs  of  signals  in  orthogonal  planes  (usually  one  in  azimuth  and 
the  other  in  elevation) . 

The  two-dimensional  development  that  follows  provides  only 
significant  features  and  results.  Extensive  derivations  are  omitted  as 
much  of  the  radar  development  parallels  previous  work,  and  a  rigorous 
development  is  lengthy.  The  scope  of  this  appendix  is  to  cover  the  two 
dimensional  radar  model  with  sufficient  detail  to  demonstrate  the  work 
done  in  this  thesis,  and  to  provide  an  avenue  for  subsequent  studies. 


Antenna  Functions 

Figure  44  shows  the  monopulse  antenna  signal  processing  circuitry 
for  a  four  horn  monopulse  radar.  Amplitude-comparison  monopulse  uses 
the  amplitude  difference  between  two  adjacent  horn  voltages  to  generate 
error  signals  (20t46) . 


Figure  44.  Monopulse  Antenna  Signal  Processing  Circuitry 
Source:  [20:46] 


Inspection  of  Figure  44  indicates  that  only  three  combinations  use 
adjacent  horns  (A+DorC+B  are  diagonal  combinations) .  The  three 
adjacent  horn  combinations  produce  a  sum  channel,  an  azimuthal 
difference  channel,  and  an  elevation  difference  channel.  The  diagonal 
combination  is  not  used,  and  is  simply  routed  to  a  load. 

Reference  4  derived  the  three  necessary  antenna  patterns  required 
for  further  processing  of  the  two-dimensional  monopulse  receive  signals 
Gaussian  beams  squinted  off  the  tracking  axis  in  the  two -orthogonal 
coordinate  directions,  azimuth  and  elevation,  were  used.  The  Gaussian 
pattern  is  almost  identical  (out  to  one-half  the  beamwidth)  to  the  Sa( ) 
pattern  resulting  from  the  uniformly  illuminated  rectangular  aperature 
that  was  developed  in  the  second  chapter  [4:268].  Since  all  angles  of 
interest  are  within  one-half  the  beamwidth,  the  squinted  Gaussian  beams 
can  be  successfuly  applied  to  this  two-dimensional  development.  The 


derived  sum  and  difference  patterns  of  squinted  beams  (for  small  angular 
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displacements )  are  [4:306]: 

V  =  2g„ 

s  ’0 

V  =  2gnk  9/BW 

e  ’0  m 

V  *  2g  k  0/BH 

a  u  m 

where : 

*  sum  voltage  pattern 

*  elevation  difference  voltage  pattern 

■  azimuthal  difference  voltage  pattern 

k  -  2sqrt ( 2 ) ln( 2 ) 8  /BW;  2sqrt(2)ln(2)0  /BW 
m  s  s 

*  0.98  if  »  -  BW/2f  0  -  BW/2 

s  s 

gQ  *  pattern  gain  constant 


Let:  gQ  “  0.5 

V  »  1 
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v#  •  e/Bw 

V  •  0/BW 
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(84) 
(84  A) 
(84B) 
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The  three  forms  of  equation  84  are  the  desired  antenna  function 
approximations  for  the  two-dimensional  monopulse  radar  model. 

Estimation  Model 

Observation  Formulation 

For  the  specific  case  of  additive  white  noise,  the  received 


waveforms  for  the  sum  and  difference  channels  are: 


r.(t)  *  As(t)V 

+ 

n  (t) 

(85) 
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r.(t)  =  As ( t ) V 

2  e 
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n2(t) 

(85A) 

r  (t)  =  As ( t ) V 

w  a 

+ 

n3(t) 

( 85B) 

Applying  the  Gram-Schmidt  orthogonalization  procedure  as  in  the 
Chapter  2  development,  the  waveforms  can  be  replaced  by  finite 
dimensional  vectors.  Characterization  of  the  received  observation  by 


way  of  a  joint  density  function  can  then  be  made. 

The  vector 

- —  — 

observations  are: 

R1  *  AV^sqrt(E)  +  N1 

(86) 

R2  *  AVesqrt(E)  +  N2 

(86A) 

R  -  AV  sqrt(E)  +  N 

J  d  J 

(86B) 

The  means  and  variances  of  the  three  observation  vectors  are: 


EfR.j/A,  8,8} 

*  Asqrt(E)V  j 

3  ' 

vartR^A.B.B) 

•  V2 

(87) 

E{R2/A,8,8} 

*  Asqrt(E)V  ; 

3 

var{R2/A,9,0J 

•  V2 

(87A) 

E{R3/A.9.0) 

“  Asqrt(E)V  ; 

3 

var(R3/A,9,0) 

-  Nfl/2  ' 

(87B) 

The  conditional  covariances  are: 

COV(RrR2/A.e,0}  -  COV{R^  ,  R^/A, 8,8}  -  COV^^/A,®,*}  -  0  (88) 

Equation  88  shows  the  vector  observations  to  be  uncorrelated. 
Further,  with  Gaussian  noise  assumed  and  linear  operations  to  develop 
the  vector  observations,  R1  R2  and  R^  will  be  Gaussian  and  uncorrelated, 
therefore  independent.  The  receiver  observation  is  then  the  joint 
density  of  three  independent  Gaussian  random  variables*  Independence 
implies  that  the  joint  density  of  R1#  R2<  and  «3  is  the  product  of  the 
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marginal  distributions.  With  Gaussian  distributions  and  statistics 
given  by  equations  87,  87a,  87b,  and  88  the  joint  conditional  density  of 
the  receiver  observation  is: 

pR(r/A,8,|»)  *  p(r^/A,O,fi)p(r^/A,9,0)p(r^/A,9, 9 ) 

•  - 

*  l(pi)NQ3 1*5exp(-1/N0[(r1  -  AVgSqrttE))2  +  (r2  -  AV^sqrttE))2  + 

<r3  -  AVasqrt(E))2] )  (89) 

Maximum-Likelihood  (ml)  estimation 

The  log-likelihood  equation  for  a  parameter  A  was  given  in  equation 
2S  of  Chapter  2.  The  extension  to  a  three  dimensional  conditional 
density  described  by  the  vectors  R^,  F  ,  and  R^  is  straightforward. 
Applying  the  procedure  described  by  equation  25  to  the  joint  conditional 
density  of  equation  89  and  solving  for  the  unknown  parameter  A  yields: 


A  .  - 
ml 


sqrt(E) 


r,v  +  +  r„V  * 

Is  2  e  3  a 

(V  )2  (V  >2+  (V  )2 

S  C  A 


Taking  the  expected  value  of  equation  90  produces  the  mean  value  of 
the  amplitude  estimate.  Performing  the  expectation,  the  amplitude 
estimate,  as  in  the  one-dimensional  case,  proves  to  be  unbiased. 

Applying  the  maximum-likelihood  procedure  of  equation  25  to  an 
estimate  of  the  elevation  angle  off  boresight  (o')  and  substitution  of 
the  amplitude  estimate  of  equation  90  into  the  result  yields  the 
elevation  error  equation.  The  elevation  error  (e^)  is: 
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\  ■  V„1BH  *  V„l*  -  r2IB"2  *  ,z‘ 

where: 

The  antenna  functions  of  equations  84,  84A,  and  84B  are 
used  in  the  development  for  elevation  error 


(91} 


Taking  the  expected  value  of  equation  91  produces  the  average 
elevation  error.  Performing  the  expectation  operation,  the  average 
error  value  is  found  to  be  equal  to  zero  when  equals  the  parameter 
9.  Assuming  zero  error,  equation  91  can  be  solved  uniquely  for  the 
elevation  maximum-likelihood  estimate.  The  result  is: 


ml 


r2lBW2  t  »2) 
r1 BW  +  r3« 


(92) 


Manipulation  of  equation  92  and  taking  the  expectation  of  both 
sides  yields  the  average  estimate  value  of  the  elevation  angle  estimate. 
Performing  the  expectation  operation,  the  maximum-likelihood  estimate 
for  the  elevation  angle  off  boresight  is  found  to  be  unbiased  (ie  the 
average  value  of  the  estimate  equals  9)'.  An  operation  similar  to  that 
for  the  equation  91  development  performed  for  the  azimuthal  angle  off, 
boresight  produces  the  azimuth  error  equation  (e^),  And  an  unbiased 
maximum-likelihood  estimate  for  the  aziomthal  angle  off  boresight. 

aa  "  VmlBW  +  Vml®  "  r3lSW?  +  ^  (93) 

.  r  IBW2  +  92] 

•ml - —  (94) 

r^w  +  r20 

When  joint  azimuth  and  elevation  angle  estimates  are  made,  the 
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estimates  are,  in  general,  coupled.  Inspection  of  equations  92  and  94 
show  that  both  estimates  have  a  scaled  Gaussian  numerator  and  a  Gaussian 
denominator  (sum  of  scaled  Gaussian  random  variables).  The  ratio  of 
jointly  independent  Gaussian  random  variables  has  a  Cauchy  density 
centered  at  zero  [11:198].  The  Cauchy  density  is  like  the  Gaussian 
density,  but  the  tails  are  off-axis.  Moments  do  not  exist  for  the 
Cauchy  density,  thus  making  it  difficult  to  determine  the  coupling 
coefficient  for  two  Cauchy  random  variables  (azimuth  and  elevation 
estimates).  This  investigation  will  make  use  of  reference  21  which 
generated  combined  azimuth  and  elevation  maximum-likelihood  estimates  in 
a  white  noise  environment.  Reference  21  shows  that  for  a  narrow-band 
signal  the  maximum-likelihood  elevation  and  azimuth  estimates  will  be 
uncoupled  if  the  two-dimensional  illumination  function  is  separable  into 
the  product  of  one-dimensional  functions  [21:44).  Since,  the 
illumination  functions  used  in  the  two-dimensional  development  are 
Gaussian,  the  estimates  are  assumed  uncorrelated. 

To  move  the  boresight  in  both  azimuth  and  elevation,  the 
relationship  between  the  two  error  equations  (equations  91  and  93)  must 
be  determined.  The  covariance  between  azimuthal  and  elevation  errors 
after  lengthy  calculations  is: 

coyCe  ,e  )  -  E{  [e  -  E(e  }][«  -  E(e  )]) 

VS  w  V  •  A  m 

•  \V.2>  *  •«!»<’  *  V.J>  *  <95» 

*  Sk 

Evaluating  equation  95  at  the  average  estimate  values  for  6  , ,  0  , : 

•  mi  mi 

cov(e  ,e  }  «  N  BW28P/2[V  2  +  V  2  ♦  V  2]  (95A) 

"  •  ^  A  i  fl 


Examination  of  equation  95A  shows  that  in  the  no-noise  environment 
the  covariance  will  be  small  for  small  angular  displacements  0,0.  Since 
target  tracking  implies  that  both  0  and  0  are  less  than  one-half  the 
beaswidth  (1.5  degrees),  the  coupling  coefficient  between  azimuth  and 
elevation  error  is  assumed  negligible.  Neglecting  the  coupling,  allows 
the  error  voltages  to  be  used  to  drive  two  separate  antenna  servo  loops, 
similar  to  the  one-dimensional  servo  loop  of  Chapter  2,  to  maintain 
track-axis  alignment  with  the  target. 

Discriminators 

Zt  is  necessary  to  express  the  two  error  equations  (equations  91 
and  93)  as  a  function  of  the  difference  between  the  estimated  and  actual 
target  locations.  Functions  of  the  diffference  allow  projected 
bores ight  movement  along  the  two  orthogonal  coordinate  axes  proportional 
to  the  computed  difference.  With  additive  noise,  the  error  equations 
are  random  quantities  as  R^,  and  are  independent  Gaussian 
distributed  random  variables.  Assuming  no-noise,  the  observations  will 
be  equal  to  their  respective  mean  values.  Using  the  mean  values  of  the 
observations  R^ ,  R2>  and  R^  and  expressing  the  elevation  difference  as 
dO,  the  elevation  error  (equation  91)  can  be  expressed  set 
.  e  ■  -  Asqrt ( E ) BWdB  -  Asqrt(E)02d0/BW 
where : 

d#  -  0  -  0 


Solving  for  d0  and  approximating  the  result j 


(96) 


d8  -  -  e  /Asqrt ( E ) BW 
c 

where s 

(i/BW)2  <<  1 

A  similar  development  for  azimuthal  error  yields  an  approximate 
expression  for  df: 

di  -  -  e  /Asqrt (E)BW  (97) 

A 

where; 

(B/BW)2  «  1 

Examination  of  equations  96  and  97  indicates  that  the  mapping  from 
estimated  error  to  step  function  proportional  to  required  boresight 
movement  is  linear.  Equations  96  an  97  will  then  be  used  to  actuate 
separate  elevation  and  azimuthal  type  X,  improved,  servo-control  systems 
to  steer  the  beam-axis  on  target. 

Ziv-Zakai  Bound  (ZZB) 

The  Ziv-Zakai  two-dimensional  bound  is  derived  by  comparing  the 
estimation  problem  with  the  optimal  detection  problem.  The  first 
development  for  the  Ziv-Zakai  bound  is  to  determine  the  error 
probability  associated  with  estimating  a  two-dimensional  target 
position.  The  next  development  is  to  compute  the  probability  of  error 
associated  with  detection  of  all  possible  combinations  of  the  two 
angular  dimensions.  This  is  an  M-ary  detection  problem  where  N  equals  2 
squared,  or.  4.  Finally,  an  inequality  must  be  established  between  the 
estimation  and  detection  problems  that  accounts  for  all  possible  target 
positions  in  the  a  priori  interval. 
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Estimation  Error  Probabilit' 


Consider  an  estimation  technique  for  the  two-dimensional  angle  off 
boresight  3  when  it  is  known  that  the  angle  is  either  at  position  or 
position  •  Figure  45  shows  the  described  situation  where  the  two- 
dimensional  position  vector  S  has  azimuthal  (0)  and  elevation  (8) 
components . 


Figure  45.  Estimation  Geometry 


The  Ziv-Zakai  approach  is  to  compare  S  with  an  average  value  of  S, 
Sav^.  The  decision  space  for  the  two  possible  positions  and  their 
average  is  shown  in  Figure  46. 


Assuming  equal  a  priori  probabilities  and  rearranging  terms: 

a  a 

Pe  -  0.5 JPC ( »  <  0,/2  +  e2/2  or  0  <  0 ,/2  +  02/2)/S,}  + 

A  A 

P{(0  >  9^/2  +  02/2  or  0  >  0^2  +  02/2)/S215  - 

A  A 

0.5  [P{  (0  <  0,/2  +  02/2  and  0  <  0^2  +  02/2)/S,}  + 

A  A 

PC (8  >  8t/2  +  92/2  and  0  >  0^2  +  02/2)/S2>}  (98A) 

Consider  the  or  terms  of  equation  98A  and  let: 


de 

•  e  -  0.; 

0, 

-  d8 

♦  04» 

0 . 

-  0.  -  d0 

2  1 

2 

1 

1 

2 

d0 

»  0  •  0  1 

2  1 

*2 

-  d0 

+  V 

*1 

«  02  -  40 

P{ (0  <  0,/2  +  02/2  or  0  <  0^2  +  02/2)/S1)  - 

A  A 

P((0  -  01  <  d0/2  or  0  -  0,  <  d»/2)/S^) 

P{(0  >  0^/2  ♦  02/2  or  0  >  0,/2  +  02/2)/S2)  - 
P{(0  -  02  >  -d0/2  or  0  -  02  >  -d0/2 )/S2) 


(98B) 

(98C) 


Imposing  symmetric  target  positions  let  S,  *  -S^  and.  fix  d0,  d0  at  the 


maximum  usable  angle  off  boresight  (ie  d0  -  d0  ■  -BW) .  Equations  98B 
and  98C  become: 


P(<0  -  0J  <  d0/2  or  0  -  01  <  df/2)/S,}  - 

A  A 

P{  1 0  -  0j  _>  BW/2/0 1  or  1 0  -  0,|  I  BW/2/0, J  (98D) 

A  -  -  - - A 

P((0  -  02  >  -d0/2  or  0  -  02  >  -d0/2)/S2)  - 

A  A 

PC  I  o  -  0,1  >  BW/2/0,  or  |0  -  0,|  _>  BW/2/0,}  (98E) 

Substitution  of  0,  ■  0  and  0,-0  into  equations  980  and  98!i  and  scaling 
by  0.5,  the  terms  can  be  added  together.  Equation  98F  is  the  result, 
and  represents  the  or  term  contribution  to  the  equation  98A  estimation 
probability  of  error.  The  or  terms  can  now  be  expressed  as: 


i24 


(98F) 


PC  I  e  -  8 1  >  BW/2  or  1 0  -  0|  £  BW/2} 

To  proceed,  equation  98F  must  be  transformed  into  a  second  moment. 
Reference  22  provides  a  two-dimensional  generalization  of  Tchebychev's 
inequality  (22: 18] : 

A  A  AAA  A 

P{|e  -  E{9) |  £  Xsqrt(var(9))  or  |0  -  E{0)|  £ Xsqrt(var(0) ) }  £ 

(1  +  sqrtd  -  p2)]/X2  (99) 

where: 

A  A 

p  is  the  correlation  coefficient  between  8  and  0 

Following  a  partial  development  in  reference  22,  and  applying  a  Biename' 
inequality  of  the  nature  of  equation  38: 

A  A 

P{ | 8  -  e|  £  aw/2  or  |f  -  0|  £  BW/2)  £ 

E{aax{(8  -  8>2,(0  -  0)2]}(1  +  sqrt( 1  -  p2)]/(BW/2)2  (99A) 

A  A 

Assuming  0  and  0  are  uncorrelated: 

P(  1 8  -  8 1  £  BW/2  or  |0  -  0|  j>  BW/2)  £  8e  (S)/BW  (99B) 

where: 

e2(S)  -  E(max ( (0  -  8)2,(0  -  0)2)} 

■  two-dimensional  mean-square  estimation  error 


Now,  consider  the  and  terms  of  equation  98 A: 

A  .A 

P((8  <  8/2  ♦  8/2  and  0  <  0/2  +  0/21/8/  - 

A  A 

p(|e  -  8|  £  BW/2/0  and  |0  -  0j  £  BW/2/0)  (100) 

A  A 

P{<e  >  0-/2  .♦  0„/2  and  0  >  tJ2  ♦  f9/2)/89)  * 

•  A  ‘  l  4  4 

A  A 

PC  I •  -  i|  £  BW/2/0  and  j 0  -  t\  £  BW/2/0)  ( 100A) 


Scaling  by  -0.5  and  adding  equations  100  and  100A  results  in  equation 
100B,  which  represents  the  and  tern  contribution  to  the  equation  98A 
estimation  probability  of  error.  The  and  terms  can  now  be  expressed  as: 

-P{|e  -  a  |  _>  BW/2  and  |0  -  0j  >  BW/2}  (100B) 

A  A 

Assuming  S  and  0  are  independent,  equation  100B  beccjes : 

A  A 

-lp{|e  -  e|  2:  bw/2jp{ 1 0  -  0|  _>  bw/2jj  (1000 

From  the  one-dimensional  development,  equation  100C  becomes  the  product 
of  scaled  one-dimensional  error  probabilities: 

-[e2(e)/(BW/2)2e2(0)/(BW/2.”‘  (100D) 

where: 

e  (0)  “  Et|e  -  a j  }  ■  mean-square  el  error 
a  2  A  2 

e  (0)  «  E{ | §  -  0 |  }  »  mean-square  az  error 

Combining  the  results  from  the  or  terms  (equation  99B)  with  the 
results  obtained  from  the  and  terms  (equation  100D),  the  expression  for 
estimation  probability  of  error  (equation  98A)  becomes: 

Pe(est)  «  (e2(3)  -  2e2(0)e2(0)/BW2]8/BW2  (101) 

Detection  error  probability 

The  detection  error  probability  is  the  error  probability  for 
deciding  whether  a  target  is  at  or  when  it  is  known  to  b_  at  one 
of  these  two  positions  with  equal  probability.  The  decision  regions  for 
the  two-dimensional  case  are  shorn  in  Figure  47. 
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Figure  47 .  Detection  Decision  Regions 


four  detection  hypotheses  sret 

V  R1  *  As^rt<B)Va(+et,+ft)  +  Ht 

R2  -  Asqrt(E)Ve{+et,+|t)  +  M2 

R3  -■  Asqrt(E)Va(+8t#+0t)  ♦  N3 

H2s  R1  -  Asqrt<E)Va<-9t,+*t)  +  » 

R2  *  A^qrt(E)Ve(-0t,+flt)  ♦  n2 

R  -  Asqrt(E)V  (-«+»)  +  « 

J  a  t  t  3 

H3:  R1  »  Asqrt(E)Va(-9t,-0t)  +  M 
R2  -  Asqrt(E)Ve(-9t,-flt)  + 

R3  »  Asqrt(E)Va(-0t,-0t)  +  ^ 

h4!  R,  *  Asqrt<E)Va(+8t,-0t>  +  h 
R2  -  Asqrt(E)V  (+9  -f  )  +  N 
R3  -  Asqrt(E)V  (+9  .-0  )  +  M 


Due  to  the  symmetry  of  the  four  decision  regions  of  Figure  47,  and 


with  ecrual  a  priori  probabilities,  the  detection  error  will  be 
equivalent  to  the  error  in  a  single  region.  Using  hypothesis  H  : 
Pe(det}  =  P{112  <  O/H^Pll^J  +  P{113  <  O/HjJPU^}  + 

p{114  <  o/h1)p{114) 


where: 


natural  logarithm  of  the  likelihood  ratio 
lntp(r/Hi)/p(r/Hj)J 


(102) 


Consider  the  general  log-likelihood  ratio  (1^): 

l±j  *  *N0"1[(r1  “  E{R1/Hi))2  +  (r2  -  EtRj/I^})2  +  (r3  -  EtR^H^)2 
-(r1  -  EtR^H^))2  -(r2  -  E^/H^))2  -{r3  -  ECR^H^)2] 

-  -(2Asqrt(E)/BWN0J  (r2(9j  -  8^>  +  r;j(0j  -  0^  ♦ 

0 .5Asqrt(E) /BW( 9^2  -  9*  +  t*  -  t*)\ 

where • 

E{R^ )  *  equation  87,  and  is  independent  of  9,0 
E(R2)  «*  equation  87a,  and  is  independent  of  0 
E(R3)  “  equation  87B,  and  is  independent  of  9 

An  equivilant  decision  rule  considers  1’^  where: 

1'^  »  lij(BWNQ)/(2Asqrt(E)) 

The. general  conditional  mean  of  1*^  follows  as: 

B{1*  :/H,)  -  Asqrt ( E ) /BW (9.(9.  -  9,)  +0.5(9  2  -  9  2  +  0  2  -  0  2)  + 


A  similar  analysis  on  the  general  conditional  variance  of  1*  ^  yields: 


var{l*ij/H1}  -  O.5NQ[(0i  -  0^)2  +  (0i  -  fl^)2] 


(104) 


To  uniquely  determine  the  density  functions  of  the  various 
combinations  of  1'^  ,  the  statistics  must  be  derived  from  the  general 
formulas  given.  With  Gaussian  statistics  for  each  of  the  joint 
conditional  densities,  1'^  is  the  sum  of  independent  Gaussian  random 
variables,  therefore  Gaussian. 


Consider  the  combination  i  *  1,  j  *  2i  i  »  (+9t,+Bt),  j  *  (”®t*+®t)* 


Etl'l2/H1}  "  2Asqrt(E)»t2/BW}  vartl'^/H.,)  -  2*0»t2 


(10S) 


P<1’l2/H1)  “  lk12(Pi)1"‘5cxPl*(kl2)"1(1‘l2  ’  2Asqrt(E)9t2/BW)2} 


where 


k  •  4N  9 
12  0°t 


Consider  the  combination  i  »  i,  j  •  3»  i  »  ( +9  , +0  ) ,  j  «*  ( -9 . . -0. ) : 

t  U  t  t 

Etl'l3/H1)  "  2ASqrt(E)(9t2+Bt2)/BW»  varU'^/H^  -  2NQ(0t2+0t2)  (106) 

Ptl'^/H,)  «  [k13(pi)]**5exp{-(k13)”1(l,13  -  2Asqrt(E)(0t2  +  0t2)/BW)2l 


where 


‘d  •  «V»tJ  ♦  *t2) 


Consider  the  combination  i  *  1,  j  ■  4f  i  *  (+9^+0^),  j  *  (♦O^-B^t 


E(l'14/H1)  -  2Asqrt(E)0t2/BW>  varJl^/Hj)  *  2N0Bt2 


(107) 


P<1'14/Hj)  -  tk14(pi)]*,5exp(-(k14)"1(l'14  2Asqrt(E)Bt2/BW)2l 


where 


*14  *  4Vt 


Using  the  appropriate  change  of  variables  in  each  of  the  three  density 
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functions  described  by  equations  105,  106,  and  107,  the  detection  error 
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(equation  102)  can  now  be  expressed  as: 

Pe(det}  =  P(1'12)Q(E(1* 12/H1)/sqrt(var{l* 12/H1>)1  + 
P{1’13)Q(E{1' 13/B1)/aqrt(var(la ^/H,} )1  + 
P(l*  14)Q[E{1*  14/H1)/sqrt(var(l* 


where : 


Q(a]  -  sqrt(2(pi) )  V—  x2/2)dx 
statistics  of  1*  are  given  in  equation  105 

statistics  of  1*  are  given  in  equation  106 

statistics  of  1'  are  given  in  equation  107 


(108) 


Substituting  the  appropriate  statistics  into  equation  108  and 
simplifying: 

P  (detj  -  0.33 [Q(A0  sqrt(2SNR)/BW)  +  Q(A0  sqrt(2SNR)/BW)  + 

©  t  t 

Q(Asqrt(9t2  +  0t2)sqrt(2SNR)/BW) ]  { 108A) 

where : 

SNR  -  E/Nq 

p(i*12}  -  p{i’13)  -  pa*14)  -1/3 

Equation  108A  describes  the  P&  of  the  best  procedure  for  deciding 
whether  a  target  is  at  S1  or  -S^,  when  it  is  known  to  be  at  one  of  these 
two  locations . with  equal  probabality.  To  compare  with  the  results 
obtained  in  one-dimensional  tracking  consider  0*0-0: 

W  to 

P  (det(O.O)}  -  .333 [Q{ 0 )  +  Q(0)  +  Q(0)J  •  0.5 

© 

This  result  is  the  same  as  obtained  in  the  one-dimensional  case  where 
the  target  separation  angle  was  zero.  Here,  target  separation  in  two 
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orthogonal  planes  is  zero  and  once  again  the  detection  error  indicates 
that  the  targets  will  be  difficult  to  distinguish. 

Detection/Estimation  Inequality 

Following  the  Ziv-Zakai  method,  a  performance  bound  can  be  derived 
by  comparing  the  estimation  problem  with  the  optimal  detection  problem. 
The  estimation  decision  rule  error  will  be  lower  bounded  by  the  error 
associated  with  an  optimal  detection  scheme.  Establishing  an 
estimation/detection  inequality: 

P  {det}  <  P  {est}  (109) 

e  —  e 

Using  the  estimation  probability  of  error  as  described  by  equation  101 
and  solving  equation  109  for  the  two-dimensional  mean-square  estimation 
error: 

e2(S)  >  BM2P  (S.-SJ/8  +  2e2(8)e2(f )/BW2  (109A) 

where : 

P  (S,-S)  is  described  by  equation  108A 
0 

*2 

e  (a)  refers  to  the  mean-square  error  for  parameter  a  and 
given  by  the  expression  in  equation  40C 

Equation  109A  is  a  lower  bound  to  the  mean-square  estimation  error 
of  any  pair  of  values  of  the  parameter  S  which  are  3W  units  apart. 
Extension  of  equation  109A  to  consider  all  possible  positions  in  the  a 
priori  interval  of  the  parameter  S,  and  using  the  worst  case  error  of 
all  possible  a  priori  values  results  in: 


e  (S)  >  max 


0  <  0  <  0 

“  -max 
0  <  i  <  0 

t  ~  max 


[S2P  (S,-S}/8  +  2e2 ( 8 )  e2  ( 0 ) /S2]  (109B) 

e  t  t 


where: 


S 


Observe  in  equation  109B  that  the  lower  bound  does  not  exist  at  S  * 
0.  Thus  a  restriction  on  the  bound  as  described  in  references  9  and  12, 
and  previously  used  in  equation  40C  of  this  thesis,  must  be  imposed . 
Further,  because  of  the  symmetry  invoiced,  the  bound  only  considers 
positive  values  of  azimuth  and  elevation.  Positive  parameter  values 
allow  the  substitution  of  the  magnitude  squared  of  the  parameter  for  the 
square  of  the  parameter.  In  vector  notation: 
s2  (pos  s)  »  |sj2  *  ssT  -  Tel  ;8t  *tl 


Redefining  the  bound  of  equation  109B: 


«2(S>  l  max  lsin< |s| )2P  {S,-S}/8+2e2(0). e2(»).  /sin( |si )2] ( 1090 
0<8  <8  8  C  Z 

t-  p 

where: 


sin (8  1 

•  O.S[sin(0 

) 

+  sin(0>] 

P 

max 

sin{0  ) 

-  0.5 (flin(| 

) 

+  sin(0) ] 

P 

max 

Equation  109C  will  be  the  form  of  the  Ziv-Zakai  bound  used  to 
predict  the  mean-square  estimation  error  for  maximum-likelihood 
estimates  of  elevation  and  azimuthal  angles. 
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Radar  simulation 


To  compute  the  mean-square  tracking  error  and  associated 
performance  bounds  of  the  two-dimensional  amplitude-comparison  monopulse 
radar,  a  simulation  program  was  designed  for  the  radar  model  as 
developed  in  this  appendix.  An  annotated  listing  of  the  simulation 
program  is  provided  in  Appendix  G .  The  program  is  a  straightforward 
extension  of  the  one-dimensional  simulation  provided  in  Appendix  B.  The 
major  modifications  and  significant  features  are: 

Three  independent  observations  are  required. 

The  Cramer-Rao  bound  is  not  computed. 

Provisions  for  moving  targets  are  removed. 

Equation  109C  (ZZB)  is  never  evaluated  at  d  »  $  -  precisely 
zero.  The  slight  offset  from  zero  allows  computation  of  the 
second  term  of  the  described  equation. 

A  storage  matrix  (16X16)  is  established  for  the  evaluation  of 
equation  1090  (ZZB)  for  the  various  segmented  angle  intervals 
in  azimuth  and  elevation.  The  stored  values  are  then  sorted  by 
row  and  column,  and  the  maximum  is  picked  off. 

Performance  Verification 

Before  comparing  the  results  with  the-  one-dimensional  tracking 
results,  a  few  comment  about  the  simulation  is  in  order.  Por  each 
computed  bound  of  the  simulation,  the  Q  function  is  evaluated  by 
integrating  the  tails  of  a  standard  Normal  density  in  an  external  call 
to  the  I MSI.  subroutine  HONOR.  One  run  calls  the  MDNOR  subroutine 


2) 

3) 

4) 

5) 
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approximately  thirteen  hundred  times.  Each  SNR  level  conducts  fifteen 
runs,  therefore  approximately  nineteen  thousand  integrations.  For  lower 
SNR  levels,  the  area  integrated  increases,  and  the  time  consumed  for 
each  run  is  quite  lengthy.  Low  SNR  runs  were  not  possible  because  of 
time  limits  on  the  computer  system  used  (CYBER). 

Figure  48  is  a  stationary  target  composite  plot  for  the  two-, 
dimensional  target  mean-square  tracking  error  (MSE) ,  and  the  Ziv-Zakai 
bound  (ZZB)  for  SNR  levels  of  30,  25,  15,  and  7  dB.  Cooparision  of 
Figure  48  with  the  previous  stationary  one-dimensional  target  tracking 
results  demonstrate  that: 

1)  The  MSE  is  much  larger  for  the  two-dimensional  tracking  case 
(reference  approximately  3X10E-5  for  one-dimension  at  15dB  SNR 
(see  Figure  30)  and  approximately  20X1 0E-5  for  two-dimensions 
at  15dB  SNR  (see  Figure  48)). 

2)  At  7dB  SNR,  the  two-dimensional  MSE  grows  unbounded.  Breaklock 
is  assumed  to  occur  at  this  point.  This  was  not  encountered  in 
the  one-dimensional  tracker. 

3)  For  all  SNR  levels  considered,  the  two-dimensional  ZZB  lower 
bounds  the  MSE  curves.  Although  not  as  tight  as  in  the  one¬ 
dimensional  tracker,  the  two-dimensional  ZZB  is  an  adequate 
performance  indicator,  and  merits  further  investigation. 
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■Figure  48.  Two-Dimensional  Target  Tracking  Error  vs  SNR  Level 


Appendix  B 
Program  Listing 


»*••*»««««**«*«*«*»»*••«»••*««««**««**«*«•**««««*»#****««»***»*******«** 


PROGRAM:  SIMULATION  OF  AMPLITUDE-COMPARISON  MONOPULSE  RADAR 

* 

PURPOSE:  COMPUTE  THE  TRACKING  ERROR  AND 

PERFORMANCE  BOUNDS 

* 

INPUT:  DESIRED  SNR  LEVEL(dB),  TARGET  MOVING  OR  STATI QNARY 

* 

OUTPUT:  3<1  BY  50)  VECTORS  CONTAINING  AVG  ERROR,  CRB,  22B 

» 

COMPILER:  FORTRAN,  VERSION  5 

* 

EXT  CALLS:  IMSL  ROUTINES  GGNML,  MDNOR 

« 

•ftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftft** 

PROGRAM  ACNON 

§ 

>  EBWS.UARlAftg 

DDOSI04  flW6<50)/W61<5l)^We2(50),PRE(14>,K30),B(3P) 

• 

DIMB6IW  SNRV<15),SER<15,58),S£RCR<13,5I),SERZ2<15,50) 

REAL  NM,NllN2,NmPlfCALC 

•  INTEGER  TO  REAL 

INTEGER  PULSE, ANGLE 

•  REAL  TO  INTEGER 

DOUBLE  PRECISION  DSEED1  ,DSEED2 

i  DOUBLE  PRECISION  SEEDS 

»  INPUT  STATEMENTS 

DATA  SWM/33.,30.,25.,21.,13.,11.,7.,5.,3.,1., 

• 

*“2.  ,-5.  ,-7.  ,-H. 

•  POSSIBLE  SNR  VALUES  <<*)• 

3 

CONTINUE 

• 

WITEO.IOOI) 

• 

READ*,  INS 

•  TARGET  MOVING  OR  QUIT 

IF(INS.EB.0)6O  TO  3100 

i 

WITE(*,I0H) 

» 

URITE(*,1028) 

• 

REA0«,I1,I2 

•  START/STOP  SNR  VALUES 

•  RADAR  PARAMETERS 

DATA  0SEE01 ,DSEED2/123457.000 ,329017.000/ 

DATA  OPAAO,R9,VEl,PU,GU/37. 3. 4900., 257.4, l.E-4, 3.0/ 

• 

DATA  PRI , SIS. ONEGA, AMP, SWR/10.E-3, 434. ,434. ,3.9,1  .E-3/ 

E6Y»W«ANPm2 

»  ENERGY 

WWQ 

i  30  GAUSSIAN,  SAMPLES 

NB*50 

•  30  GAUSSIAN  SAMPLES 

DO  309  IP>I1,I2 

t  SNR  DO  LOOP 

SNB»9UW(IP) 

• 

9MNB>t0.*«9Rt/10.) 

•  CONVERT  TO  NWERIC 

,  ROLO-M 

*  INITIAL  RANGE 

DO  200  >1,13  • 

•  >R1N  NUMBER 

CALL  69tf.(0SEEDl  ,Nft,R) 

»  NGMAL  DEVIATE  GENERATOR! 

CALL  6GW1<DSEE02,N8,8> 

•  NORM.  DEVIATE  GENERATOR! 

•  PIACE  TARGET  IN  BEAMU1DTM 

BS1TE-0.8 

•  STATIONARY 

TR6T*9.8t 

•  STATIONARY 

IF( IMS.EO. 1)8SITE*1 .5 

•  MOVING 

If(INS.E8.1)TRGT«0.8 

»  MOVING 

•  TRACKING  CALCINATIONS 

TKETA-TRGT-BS1TE 

•  TARGET  ANGLE  OF  ARRIVAL  • 

I-B 

• 

• 

00  191  PULSE-1,31 

•  SB  PULSES 

• 

PULSEC-flQAT(PULSE) 

•  CONVERT  TO  NWERIC 

• 

»e«amwM2*<va»wi*w.sEC)M2) 

IF(IMS.EB.l.)!l4MEI*<<R0Lfl/ll€W)M4)*SNSNEW 

i  mi  wk-rfi  mi  fnflri 
iWFvwjim 

NAMNSflRT(NMl) 

1*1+1 

NHK1HMAMP 

KM<1)«W* 

VD*$D«1.88*Tt£TA/BW) 

VS»<SBRT<2J/2J*<1.+COS<2.094*THETA/8W)) 

RJ*AMP*$8RT<ESY)*VS«N1 

R2-HMP»S«T(E5Y)*WW<2 

£RINll*VD-R2*VS 

!F(ERR/<SBRT<2.*EGT)*AMP)  .GE.1.99)  EJBW.99*S0RT(2.»EGY)*AMP 

lF<EHR/<SflBT<2.*£SY)*ANP).LE.-4.9?>  EBR-4.99*S8RT(2.*£Gr)*AMP 

DTIfT**-<BU/l  .8I)MSIN<  EM/(S9RT(2 .  <€6T)  «W) ) 

S£1MW.*^<-SI6*PRI)*<C0S<(XEGA»PR1>^IN<(XE6A*PRI>> 

IF<VD*ERR.6TJ.0)OTHET*-OTHETA 

BSITWSm+SERVO*OTHETA 

lF<JNS.E8.1)7R6faATAN<(VEL*PULSEC*PRl)/RI)»0PRA0 

THCTA-TRST-BSITE 

BMUrtPRAP 


COMPUTE  CHAN6JN6  RAN6E  • 
OXPUTE  CHANGING  SNR  < 
NOISE  POWER  • 
NOISE  AMPLITUDE  • 

t 

U91  SAMPLE  OBS  1  * 
WON  SAMPLE  OBS  2  » 
AhTBMA  BIFF  WEIGHTING  • 
ANTBMA  SIX  WEIGHTING  * 
SIX  CHAMCL  OBSEWATION  * 
OIFF  CHANEL  08SEWATION* 
ERROR  S1BML  * 
LIMIT  POSITIVE  ERROR  • 
LIMIT  NEGATIVE  ERROR  • 
STEP  PROPORTIONAL  TO  ERR* 
SEJM)  RESPONSE  • 
TRACK  PROPERLY  • 
MOVE  BORESIGHT  • 
OWNING  TARGET  POSITION* 
NBA  ANGLE  OF  ARRIVAL  • 


THETfrTKTA/DPIM 

SEKJ,PULSE)»T)CTA**2 

OOOH2.*SMNBMW**2)*<2.8A5+1 .77*C0S<3.74*THETA/HW> 
*-l  .l9S*C0S<4.18t*TICTArtW)) 

QNNH**2)/DMM 

SERONJ,PULSEMIi 

T1CTN(ASIN«SIN<BH^)*SIN<ABS{THETA) )  )«fl  .5))*0PRA» 

BMMPRAO 

THETA*THETA*OPRAO 

OATA  PRE/ 14*4.0/ 

ANGLE*1.+ABS<  IMT(10.*THETP)) 

00  90  W1  .ANGLE 
NCALC*<  FLOAT (AN6LE>/1 00 . )  <N 
IA»5./2*2.*CQS(2.0941CAlC/BW)+l.5*CQS(4.18filCALC/8W) 
taC8S(3i7MCAlC/Bi) 

RHO»<l.5+2.»COS(2.094KALC/BW)+4.5*CflS(4.1881CALC/8W) 
*+COS(  3 . 74NCALC/BW)/W 
AX*<S0RT((1.-RH0)W*SNRNEU/2.))«AN» 

CALL  M0N0R(AXfPX) 

ax*i.-px 

NCALONCALC/OPRAO 
PRE(N)*flX»<  SIN<  1 .  ACALC) )  **2 
OXTINUE 
00  94  JK*I,13 

01  92  KJ*JK1,14 
IF(PRE<JK).LE.PRE<KJ))60  TO  92 
STORE*PRE(JX) 

PRE(JX)"PRE<KJ) 


TRACT IN6  ERROR 
sawffo  ERROR 


OTBL 


*?*  a-.1 


jm 


DENOMINATOR 

BOIM 


CIXPUTE  2IV-ZAKAI  BOIXO 
ARITHMETIC  MEAN  ANGLE 


INITIALIZE 

RANGE  XII  FOR  DO  LOOP 
CHECK  FROM  0-THETP 


ANTBMA  WEIGHTING 

CORRELATION  COEFFICIENT 
INPUT  PARAMETER  AX 
COMPUTE  AREA  -INF  TO  AX 
0  FIXCTIBi  <P»> 

SOM 

SORT  ROUTINE  PREO 


Ptt(KJ)«STORE 
92  0NT1NUE 
94  CONTINUE 
ZZWREUO) 
SERZZCJ.PUISEWZB 


PRE(IO)  IS  MAXIMUM 


100  CONTINUE 
201  CONTINUE 

U8ITE(f, 1030) SNR 
URITE(i,ll40) 

00  400  INI, 30 
SUH.I 

stm««.o 

SIMM. I 
DO  300  L»l,15 
SUNSUNSEBd,*) 

SlNl*$IMl*SEBCRa,IO 

S(M2>StM2«SERZZa,X) 

301  CONTINUE 
4W6<KV>StH/l3. 

4M6KK)«SIM1/13. 

AW2«)«SII12/13. 

UUTE(*,1030MM6<U^G1(X) /WG2(K) 

400  CWTINUE 
3N  CONTINUE 
60  TO  3 

1000  FORNAT( 1X,'H(WIN6  TGT.STAT  TST,  OR  OUIT<  1  :WWE,2:STAt ,0:OUIT>') 
101*  FOamT<lX,'9NNW/33i30,23.20,l3,li,7t5.3fl,-2,-5t-?,-10,-2a/'^> 
1020  FONIATUX, 'BITER  MW  START/STOP  POINTSdE  START-1, STOP-15)') 
1030  F0WAT(/,1X,'SNR  LWEL  (OB)  »  ',£11.5,/) 

1040  FONIAKW/rtUSGER'.W/FWSSERCR'.OX.'flMGSERZZ') 

1030  F0I91AT(1QX,3(3X.E11 .5)) 

3000  STOP 
ENO 


HEADER  * 

AUERA6E  ERRORS  3UER  13  RIMS' 


OUTPUT  RESULTS 


ONLY  MAY  OUT  IS  0  IN  IHS< 
FOtHAT  STATEMENTS  IN/OUT  \ 


Appendix  C 
Program  Listing 


ft********************************************************************** 


PROGRAM:  SIMULATION  OF  "SLOW*  AND  "RAPID"  AMPLITUDE  FADING  * 

PURPOSE:  COMPUTE  THE  TRACKING  ERROR  AND  PERFORMANCE  BOUNDS  * 

INPUT:  DESIRED  SNR  LEVEL<dB),  TARGET  MOVING  OR  STATIONARY,  * 

TYPE  OF  FADE  CHARACTERISTIC  DESIRED  * 

OUTPUT:  3<1  BY  50)  VECTORS  CONTAINING  AVG  ERROR,  CRB,  ZZB  * 

COMPILER:  FORTRAN,  VERSION  5  * 

EXT  CALLS:  IMSL  ROUTINES  SGML,  MDNOR,  GGWIB  * 


ft************************************************ 


ft******************** 


program  fade 

0IMB6I W  AW 6(51)  ,AW61  (51)  /U62(S8)  ,PRE< 14)  ,R<50)  ,B(50)  ,A<50) 
D1HBCI0N  9llV<15),SER<15,5t),SERCR(15,59),S£R2Z(lS,58) 

REAL  NNR,N1  ,N2,NAHP,NCALC 
INTEGER  PULSE ^N6LE,SUERLS 
DOUBLE  PRECISION  DSEED1 ,DSEE92 

DATA  SNRV/33.,3l,,25.,;£^,15.,10.t7.,3.,3.,l., 
»-2.,-5.,-7.,-ll.,-2l./ 

3  CONTINUE 
URITE(»,1000) 

READ*, INS 

IF< INS. EQ. 1)60  TO  SIM 
URITE(*,1I1I) 

URITEl*,1020) 

KAO*, II, 12 
WRITE<*,1825> 

REA0*,SUERL6 

DATA  DSEED1  ,D$EED2/123457.0D0, 325017. 000/ 

DATA  DPRAO,RO,VEL,W, 8U/37.3, 4990. ,237.4, l.E-4, 3.0/ 

OATA  PRI,S16*d€M,S(W18.E-3,434..434.,l.E-3/ 

ALPHA" 3. 9994 
S»2 

T«SQRT<2.)*ALPHA 

NR-50 

NO*  50 

IF  (SUERL6.EQ.1)NA"I3 
DO  300  IP*I1,I2 
SNB*SNIW<1P) 

SNRNEU"I0.**<  SNR/10.) 

ROOM 

00  200  >1,13 

CALL  6®M.<0SEE01,NR,I> 

CALL  G9W.(DS£ID2,NB,B> 

CALL  6flWB<0$EE01,S,NM> 


DECLARE  VARIABLES 
SO  AMPLITUDE  SAMPLES 

FADE  TYPE:  SUERL6 
INPUT  STATEMENTS 


INPUT  FACE  TYPE 
RADAR  PARAMETERS 


RAYLE1SH  PARANETER-STOEV 
SHAPE  PARAMETER 
SCALE  PAMMETER 


58  RAYLEI6H  SAMPLES 
IS  SAMPLES  SUJU  FADE 


UEIBULL  DEVIATES 


PLACE  TAR6ET  IN  BEAHWIDTH 


Bsrrw.8 

TR6M.80 

1F(1NS.E0.1)BSITE»1.5 

IF<INS.£Q.1)TRGT*0.0 

THETA*TRGT-fiSI7E 

I«B 

CO  101  PULSE*1 ,50 
PULSEOFLQAT  ( PULSE) 

MU^ORT(ROH24<va*Pin>PaSEC)H2) 

IF(lMS.EQ.l.)SNfNEkH<R0LD/MEU)**4)*SffiNEW 

yaj^gtm/enmaj  - 
i TWFgfw  IWVp 

NM^SORT(NPUR) 

W*l 

NMKDMHP 

N24(1)4MV 

IF(SUERL6JC.1)AN^T<A<1) 

1F(SUERL6.EQ.DAW»T<A<J) 

EW«PU*AN?**2 

V0-SIN<1.88*THETA^L-) 

V*<S0RT<2.)/2.)»<l.*COS<2.0M»THETA/BU)) 

R1*AMP*S0RT  <  ESV)  *JS*N1 
R>AMP*S8RT(E6Y)<M>+N2 
EJBMM«AHtt*VS 

IF<ER*/<S0RT<2.*£6r>«AHP>.GE.I.?9>  ERJM).»«S0RT<2.*E6Y)«AMP 
IF(ERR/<SQRT(2.*EGY)*AHP) .LE.~4.99)  E»M).99*S*T<2.*I6T>«AHP 
OTHETAMOU/1.88>*ASIN<EJ8WSORT<2.*£GY)«AMP>> 

SEJW0»1 .B-£XP<-S16*PRI)*(C0S(QNEGA*PRI)-SIN<Q€GA*PR!)) 

IF(VO*€RR.ST.O.O)DTHETA«-0THETA 

BS1TE«8SITE+SEJW0*DTHETA 

IFdKS.EQ.l  )TR6T«ATAN<  <V£L*PULSEC*PRI  )/R0 )  »0PAA0 

THETA»TR6T-BSITE 

aWW/OPRAO 

THEWTHETA/OPRAO 

SER( J , PULSE )«THETA*»2 

D0W^<2.»SNf*€W«AMP«2)i(2.8i5M.77»COS(3.74»THETA/0U) 

*-l .095*COS<4 . i88»THETA/0U>  > 

CRB»<BUf*2)/DB«N 
SERCR(J, PULSE )*CRB 

THETP»<ASIN<  (SIN(BU/2)+SIN<A8S(TH£TA)  ))»0.3))*DPWD 

O»8U*0Pf!A0 

THfc7A«Tt€TA»0PRA0 

DATA  PRE/1 4*0.0/ 

AN6LE-1 .  «ABSUNT<  It .  »THETP> ) 

00  WW1.AN6LE 
NCALO<  FLOAT  (ANGLE)/!  01 . )  44 
IM./2t2.*C0S(2.09MttALC/SU)+0.3*COS<4.188NtiALC/3U) 
*-C0S(3.74«CALC/M> 

RMO*(0.3*2.*C0S(2.0944CALC/BU)*l.3*COS(4.168*NCALC/BU) 


TRACKINS  CALCULATIONS 


‘RAPID*  FADE  AMP 
‘SLOW  FADE  AMP 
E6Y  CHANGES  WITH  A 


TRACKING  ERROR 

CfWHER-MO  BOUND 


■jLintUl 


SiukUl 


•  21V-ZAKA1  BOUND 


110* 


\ 


*+C0S(3.76*NCALC/BU)/U 
AX=(SaRT((l.HM))«U*SN»«U/2.))i<W 
CAU.  HDNOR<AX,PX) 

®e»i.-px 

NCALONCALC/DPRAO 
PRE<N)«QX»(S1N(1  .*NCALC))i«2 
90  CONTINUE 
DO  94  JK*1 , 1 5 
JKWK+1 
DO  92  KJ«JK1,14 
IF<PRE(JK).LE.PRE<KJ))60  TO  92 
ST0RE»PR£( JK) 

PRE(JX)*PRE<KJ) 

PRE(KJ)*STOBE 
92  CONTINUE 
94  CONTINUE 
ZZW»RE<1« 

SER2Z(J,PULSEWZB 

100  CIHTINUE 
200  CtMTINUE 

URITE(*,1030)SNR 

URITE(«,1040) 

00  400  K-1,30 
SUM.O 
swi-o.e 
suoo.o 

DO  300  L*l,15 
SU^$IH+SER<L,K) 

SUNl*SWl*SERCRa,K> 

SW2»SlH2+SERZZa,K> 

300  CONTINUE 
.  AYG<K)*SUN/l5, 

AW6KK)*SUN1/15. 

AWG2<K)«SUN2/13. 

URITE<*,10S0)AW6(K)  ,AU61(K) ,AWS2(K) 

400  CONTINUE 
500  CONTINUE 
60  TO  3 

1000  F0WAT<lX,'fWlN6  TST.STAT  tST,  OR  0UlT(i:N0VE,2jSTAT,6:0UIT)') 

1010  FMNAT<lX,'SNRW/35,30,25,20,13,10,7,3f3fl,-2,-5,-7.-t8»'2D/V> 

1020  FORMAT!  IX, 'BJTER  SNW  START/STOP  POINTS! IE  START-1, STQP*13)'> 
1023  FORMAT!  IX, 'ENTER  DESIRED  FADE  MODEL  <SL0U«l  ,FAST»2)') 

1030  F0WAT!/,1X,'SNR  LE^EL  <OB)  «  ',Ell.5,/> 

1040  FORMAT!  14X,'AW6S£R' ,8X,'AWGSERCR'  ^'AWfiSERZZ') 

1030  FOWAT!10X,3!3X,E11.5» 

5000  STOP 

mo 


AMEBAGE  ERRORS  WER  15  R1HS< 


OUTPUT  RESULTS 


FORMAT  STATEHBfTS  IN/OUT 


************************************************************************ 


*  PROGRAMS  SMOOTHING  OF  'RAPID*  AMPLITUDE  FADING  OVER  N  PULSES  * 

»  PURPOSE:  COMPUTE  THE  TRACKING  ERROR  AND  PERFORMANCE  BOUNDS  * 

*  INPUT:  DESIRED  SNR  LEVEL*  dB) ,  TARGET  MOVING  OR  STATIONARY,  * 

*  NUMBER  OF  INTEGRATED  PULSES  DESIRED  * 

*  OUTPUT:  3<1  BY  50)  VECTORS  CONTAINING  AVG  ERROR,  CRB,  22B  .  * 

*  COMPILER:  FORTRAN,  VERSION  5  * 

*  EXT  CALLS:  IMSL  ROUTINES  GGNML,  MDNOR,  GGUIB  * 


•*•«*««**#«***«*«**•*****««*«•***«**««««•«#**•*««•#***************«•**** 


PROGRAM  FAOFIX 

DIMB4SKM  AV6<50)  ,AV61 (50)  fAV62(50>  ,PRE<16)  ,R<  50)  ,B(50)  ,A<5B> 
DIMENSION  SNRU(15) ,SER<15,50) ,SERCR<15,50) ,SER22( 15,50) 

REAL  NPUR,N1,N2,NAMP,NGALC 
INTEGER  PULSE, ANGLE 
DOUBLE  PRECISION  DSEED1 ,DSEED2 
DATA  SNRU/35.,30.,25.,20.,1S.,18.,7MS.,3M1., 
*-2.,-5.,-7.,-10.,-20./ 

5  CONTINUE 
WRITE**, 1000) 

READ*,  IMS 

IF(IMS.EB.0)6O  TO  5000 
WRITE**, 1010) 

WRi  fEO,  1020) 

REA0*,I1,I2 
WRITE**, 1025) 

WRITE<» , 1027) 

READ*,M 

DATA  0SEED1 , DSEE02/ 1 23457 . 0 DO ,  3250 1 7 .  ODD/ 

DATA  DPRAD,R0,VEL,PW,BU/57.3,4?0fl,,257.4,l.E*4f3.fl/ 

DATA  PRI,SIG,01ffiGA,SPUR/10.E-3,634.,434Ml.E-3/ 

ALPKM.9894 

S«2 

T=S0RT*2.)*ALPHA 

NR=5Q 

NB-50 

NA»50 

DO  500  IP-11,12 
SNR*SNRV(IP) 

SNR€W«10.*»<  SNR/10.) 

ROLD*R0 

00  200  >1,15 

CALL  69H(DSEED1,NR,R) 

CALL  6GWL*DSEED2,N8,8) 

CALL  60UIB(DSEED1 ,S,NA,A> 

BS1TE*0.0 

TRffM.80 


/ 

I 


i 


f 

* 


i  V* 


*.•  v 


1F<IMS.EQ.1)B$ITE*1.5 

IF(1MS.EQ.1)TRGT*0.0 

ti€ta*trgt-bsite 

1*0 

NO*1NT(50./FLCY  <N)> 

DO  100  PULSE»1,MQ 

STSlfH.I 

PULSEOO 

00  150  1100*1,8 
PULSEWULSEC+l 

Ha^80RT(ROH24(VB.*PR!>PU.SEC)H2) 

IFdMS.EB.l.)SNRNEl^<<R0lD/f0€U)«4)*$N&€U 

lUkCBJP/ttAfil 

nn^wnpv 

NAMMGRTOWR) 

NMK!)*W* 

N2*8(l)«Wf 

MNT<A<1> 

E6Y*fWm*2 
VD*SIN(  1 .  B8*T)€WBM) 
V9*(SQRT(2.)/2.)«(l.tC0S(2.094»TliTA^U)) 
RMMP*S0RT<E6Y)iVS*« 

KMHP*S0KT(E6Y)XflHN2 

EMWlHW-KiVS 

1F(ERR/(SQRT(2.*£SY)IAHP).6E.Q.99)  ERM.W*SQFT<2.«€6Y>«AMP 
IF(ERR/(S0RT(2.*£6Y)<ANP).LE.~0.99)  ER«H-0.mSQfiT<2.*E6Y)«AHP 
0THETfr*-<BU/1.88)*ASIN<ERR/<SQRT(2.*£GY)*ANP)) 

SERWM .8-EXP<-SI6*PRI)i(COS(QMEGA*PRl)-SlN(OHEGAtPRl)) 
]F(UO«ERR.6T.O.O)OTNETAp-ffTlCTA 
OTOfMTSUMTHETA 
ISO  CONTINUE 

DnCTAMTTSW/FLOAT  (H) 

BSITE*BS1TE+SEJW0«DTHETA 

I  F<  IMS .  EQ .  1  )TWT*ATAN(  (VELtPUL5EC»PRI  )/R0 )  «0PAA0 

THETA*TWT-8S1TE 

BWU/DPRAO 

THETA«THETA/DPRAO 

SER<J,PULS£)*THETA*«2 

DeW»(2.«SMWa«AhP»»2)*(2.W5»1.77»C0S{3.74»THETA/ftl) 

*- 1  .OW*COS<  4 . 1 88*THETA/8U > ) 

CR9»<0U*«2)/D0*H 

SERCR<J,PULSE)*CRB 

THETP*(AS1N(  (SlN(8U/2)+SIN<AB$(THETA) ) )  XI  .3)  HOPRAD 

BLHWOPRAO 

Tl€TA*THETA»OPRAO 

DATA  PRE/1410.0/ 

AN6LE*l.+A8S<lNT(10.iTKETP)) 

DO  90  F*»l  tAN6L£ 

NCALO<  FLOAT (ANGLE)/!  00. 

1*5  ./2*2 .  *C0S<  2 .  094*CALC/BU>  *0 . 3«COS<  4 . 1 88*<CALC/8U> 
•~C0S(3.764ffiALC/BU) 


*143- 


TRACK1N6  CALCULATIONS 

MWWHBER  OF  BSITE  HOMES 
INITIALIZE  ERROR  SUM 
9tt0THIN6  DO  LOOP 


ACCUMULATE  ERROR  SUN 
AVERAGE  ERROR 

TRACKING  ERROR  k  BOUNDS 
TRACKING  ERROR 

CRANER-AAO  BOUNO 


RHO*(  0.5+2.  *C0S!2 .094«CALC/0UMO . 5*C0S(  4 . 1 88#!CALC/BU) 

i 

> 

wCOS!3.76«CALC/BU)/U 

> 

* 

«t*(S«T((l.-RHO)*#SNi»€U/2.))*flMP 

• 

CALL  HONOR!AX,PX) 

> 

* 

OW.-PX 

* 

• 

NCALWCALC/DPMD 

« 

• 

PRE!N)*fiX»!SINd  .iNCALC)  )«2 

• 

* 

91  CONTINUE 

• 

* 

DO  94  JK*1,15 

• 

* 

JK1-JX+1 

* 

* 

DO  92  KMK1,U 

» 

• 

IF<PRE(JK).L£.P8E<KJ))60  TO  92 

« 

STORE-PRE!#' 

• 

* 

Pft£(JK)3PftE(Xj) 

* 

* 

PRE!KJ)*STORE 

t 

« 

92  CONTINUE 

• 

i 

94  CONTINUE 

• 

« 

ZZWRE!14) 

i 

ZIV-ZAKAI  BOUND 

* 

SERZZ!J,PULSE)»ZZB 

i 

> 

ISO  CONTINUE 

• 

* 

211  CONTINUE 

• 

* 

WRITE!*, 1038)90 

§ 

WRITE!*, 1048) 

• 

> 

• 

AVERAGE  ERRORS  (WER  IS  RUNS* 

D0  400  R»1,58 

• 

S9W.9 

• 

SUIM.I 

• 

SON2M.O 

• 

DO  380  L*1 ,15 

S9MUN4SER!L,K) 

• 

SW1«SUM1*SERCR!L,K) 

• 

SIH2»SW2*SERZZ!L,I!) 

i 

300  CONTINUE 

• 

AW6!K)*SJIt/lS. 

i 

1 

AV61<K>«SIM/15. 

• 

AVG2!K)«SU12/15. 

* 

i 

OUTPUT  RESULTS 

URITE(»,1050)AU6(K),AW61(K),AU62(K) 

• 

, 

400  CONTINUE 

• 

300  CONTINUE 

» 

• 

60  TO  3 

• 

■ 

* 

i 

I 

iH 

1000  FORMATdX, 'MOVING  T6T,STAT  TGT,  OR  QUlTdiMOVE,2:STAT,0!BUIT)') 

• 

1010  FORMATdX, 'SNRV/35, 30, 23, 20,15, 10, 7, 5, 3, 1,-2, -5, -7, *10, -20/' ,/> 

• 

1020  FORMATdX, 'ENTER  9*W  START/STQP  POINTSdE  START»1,ST0P»13>') 

• 

1023  FORMAT! IX, 'PULSES  INTEGRATED/1 ,2,3,10,23,50/' ,/) 

• 

■ 

1027  FOWAT! IX, 'ENTER  DESIRED  NUM8ER  OF  INTEGRATED  PULSES') 

• 

1030  FORMAT!/, IX, 'SNR  LEVEL  <DB>  »  '(Ell,5,/) 

• 

1040  F0mAT(14X,'MGSER',SX,'AV6SERCR',4X,'AV6SERZZ') 

• 

1050  F0flHAT!10X,3!3C,E11.3» 

• 

■ 

3000  STOP 

90 

. 

/ 


Appendix  E 
Program  Listing 


it********************************************************************** 


» 

* 

» 

» 

* 

» 


PROGRAM: 

PURPOSE: 

INPUT: 

OUTPUT: 
COMPILER: 
EXT  CALLS: 


SIMULATION  OF  MULTIPATH  EFFECTS  (TERRAIN  BOUNCE) 
COMPUTE  THE  TRACKING  ERROR  AND  PERFORMANCE  BOUNDS 
DESIRED  SNR  LEVEL(dB),  TARGET  MOVING  OR  STATIONARY, 
AVG  OCCURANCES  IN  THE  PRI ,  SPECULAR  OR  DIFFUSE  MULTI 
3(1  BY  50)  VECTORS  CONTAINING  AUG  ERROR,  CRB,  ZZB 
FORTRAN,  VERSION  5 

IMSL  ROUTINES  GGN1L,  MDNOR,  GGPOS,  GGEXN,  6GUBS 


* 


****«#•»•#«*»«***«•»#*****•****•»•«**•*»*»**•***## 


***#*•*«#**•«*•«**«* 


PROGRAM  MULFTH 

DIM0JS1W  AV6(50)  ,AM61  (50)  ,AM62(30)  ,PRE<14)  ,R(12S0)  ,0(1250) 
DIMBISION  SNRK13),SER<13,5Q),SERCR<13,50),SER2Z(1S,56) 
DINBISION  1R(S0) ,E<24) ,DT1NE<24) ,U(25) ,ANS<25. 

REAL  NPWR  ,N1  ,N2  ,NAMP  ,NCALC ,  PULSEC 

1NTE6ER  PULSE, ANGLE, REF, FIA6 

DOUBLE  PRECISIW  DSEED1 , DSEED2 , DSEED3 , DSEED4 , DSEEDS 

DATA  9W/3S.,30.,23.,28.,15.,lflM7.,S.,3.,l., 
*-2.,-5.,-7.,-lO.,-20./ 

3  CONTINUE 
URITE«,100I) 

READ*, IMS 

IF< IMS.EQ.0)6Q  TO  3001 
URITEO.lttlO) 

URITE(*,1020) 

READ*, 11, 12 

DATA  DSEED1 , DSEED2/1234S7.0D0, 323017. ODQ/ 

DATA  DSEED3,DSEED4,DSEED5/123457.000,  1 23457 . ODD,  1 23457 . 000/ 
DATA  DPRAO,RO,UEL,PV,6U/57. 3, 4900., 257. 4,1. £-4,3.0/ 

DATA  PRI, SIG, ONEGA, AMP,SPWR/10;E-3, 434., 434., 5.0, l.E-3/ 

NR-1230 

N8-1230 

NP*50 

FLAG-i 

DO  300  IP»I1,I2 
SNR-SWWUP) 

SNWEU»10.**(9IR/18.) 

ROOM 

URITE«,1023) 

URITE(*,1023) 

READ*, RIAN 
URITE<*,1027) 

READ*, REF 

1F(RIAH.LE. 0.0)60  TO  3 
IF(RLAN.6T. 11.0)60  TO  3 
00  200  >1,13 
CALL  69tl(DSEEDl,NR,R> 


DECLARE  VARIABLES 
DIN  FOR  24  REFLECT,  NAX 

POIS,EXP,UilF  DEVIATES 

FLAG  MONITORS  OVERLAP 
ADOITIOML  DBLEPRE  SEEDS 
INPUT  STATBIBTTS 


ONLY  MAY  OUT  IS  0  IN  INS 


RADAR  PARAMETERS 


PRINARY  *24  REFLECTED, NAX* 
PRIMARY  *24  REFLECTED, NAX « 
50  POISSON  DEVIATES  • 
INITIALIZE  FLAG  TO  0  • 


INPUT  AMS  0CCUM4CES/PR!  < 
LAMBDA  UNITED  TO  10,  PO$< 
RUN*P01SStN  PARAMETER  i 

.  i 

SPECULAR  OR  DIFFUSE  REFl  i 
LOWER  LIMIT,  MUST  BE  POS  i 
UPPER  LIMIT,  AMS  {  10  < 
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CALL  8WKDSEED2,»fB> 

CALL  6GPQ$<RLAM,DSEED3,HP,lRfIER) 

BSITEM.I 

7R6TM.M 

1F<INS.E0.1)BSITE*1.5 

1F<INS.E0.1)7RST«I.I 

TlCTA»TRST-8SrTE 

M 

DO  IN  PULSM,5I 
LASWKFaSE)*l 
NEX»ll<  PULSE) 

*Hl.lf-2>/RQAT<LAST> 

CALL  66GBMDSEEM^t^CX,0 
WHO 

Ml  GGUBSLDSEEDSfNUytl) 

1FOCX.LE.D8C  TO  69 
DO  61  LL-I.NEX-1 

ipmi*i 

N  39  KJWPljNEX 
IF(E<LL).LE.E(KX))60  TO  33 
TBM(LL) 

E<LD-€<KK) 

EHKMBT 
S3  CONTINUE 
61  CJKTINUE 

63  N  83  NULTW,NEX 
IF(FIA6J€.1)60  TO  67 
FIAO-O 
30  TO  83 
67  CtNTINUE 

IF<HULTI.EB.0)E<NULTI)*«.8 

E<LAST)«PR1 

!F<c<NULTI*l).GE.PRI)GO  TO  83 
DT1NE(NULTI*I)«£(MULTM)-E(H0LTI) 
C0EFFM.4 
00  70  II»1,NU 
AN6<ID<U<I  1X6.283 
71  CWT1NUE 

IFINULTl .£B.O)COEFF»1 .1 

IF(0TINE(NULTl*l).fiT.2«fy)S0  TO  73 

UNWOEFF 

AN61<AN6(MULT1) 

VK2-0.4 

AN62"AN6(HULTI*1) 

UHT»S8BT<(VN1«C08(AN61)*LN2«C0S<AN62))M2 

**<VMI«SIN<AN61)*VM2»S1N<ANG2))m2) 

COEFF-UNT 

FLAW 

73  CONTINUE 


POISSON  DEVIATES 
flflg  TAR6ET  IN  BBH610TH 


TRACKING  CALCULATIONS 


TOTAL  PULSE-raiSS«DI  RECT 
•  RE8  INTERARRIVAL  TINES 
COVOTE  AW  >  PR1/PULSES 
G8CBATE  NEX  INTRARRIVlS 
m  8  IffJH  DB7IATES  REB 
GENERATE  NO  AM  PHASES 
DONT  NEED  TO  SORT 
SORT  INTERARRIVAL  TINES 


NULTIPATH  P8I  DO  LOOP 
DIRECT  4  REFLCT  DO  LOOP 
FLAW  ON  2ND  OVRLP  PULSE 
RESET  FLAG  VARIABLE 
OBSRV  I  PULSE  ON  tWERLAP 

DIRECT  PULSE  ARRIVAL  TINE 
LAST  ARR1VAL*4€XT  PRI 
UNIT  REF  FRAME  TO  PRI 
DTINE  BEIUEBt  AOJ  PULSES 
REF  COEFFICIBiT  SPEC.DIFF 
PHASES  UNIFORM  0-2PI 
2P1*U<8,1) 

DIRECT  REF  COEFF  MAO 
compute  PULSE  overlap 
IF  PULSE  OVRLP  00  ELSE  73 
FIRST  PULSE  NAG  COEFF 
FIRST  PULSE  PHASE  ANGLE 
SECONO  PULSE  NAG  COEFF 
SECONO  PULSE  PHASE  AN6LE 

RESULTANT  NAG  COEFF 

FLAG  SET  OMITS  NEXT  REF 
M  PULSE  OVERIAP  CALC 
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PULSEC-fUMT(  PULSE) 

T1NE«PRU<PUIS£C-1  .)4E(IUTI*1) 

■W»S«T(W^<Va«TlME)M2) 

IF(IHS.EB.l  .)SH#CI*,<(WLI)/1WEU)h4)»SHWBU 

UBJ^CflUD/QAJBi 
irWPrinRv  SnmDI 

NAMP-SORTONR*) 

W*I 

NMKD4MHP 

IF(NULT].EQ.0)6O  70  94 

1F(REF.EB.1)THET<H)0UTI)»1.5 

IF<REF.E0.2)TICTAMKI1UITI)«180. 

1F(THETA.6T. 1.3)60  TO  8S 
M  CONTINUE 

E6T»fV*<C0EFF«AHP)w2 

^SWl.SWTICTA^W) 

VXSOKT(2.)/2.)»<1.*COS(2.0?4*T1CTA/BM)) 

R1*C0EFF*AHP»$QRT<E6T)*VS*N1 

R>C0EFF*AHP«SQRT(EGY)*A>*N2 

ERfMtlMH2*r$ 

IF(ERR/(SQRT(2.<E6T)<C0EFFtAHP).G£.l.99) 
•ERM.m$0RT<2.*E6r>«C0EFF«ANP 
lF<ERR/<SQRT<2.»E6r)«C0EFF*AHP) .LE.-6 .99) 
»£RR»-4.9P»SQRT<2.*£ST)*CflEFF»AMP 
OTICTfrMaVl .88»AS1N<ERR/<SQRT(2.*E6Y)*C0EFF<AMP)) 
SElWW.I”EXP<“SI6*OTINE(NUL71*l))*<C0S((I€6A«07IHEO10L7I*I)) 
*-SII«(ICGA«OTlM(NULTI*l))) 

IF(VDfEtt.6T.0.0)0THETAHmCTA 
BSlTE-6SITE«SEiMHKmCTA 
03  CONTINUE 

IF(  INS  .EO.I  )TR6WTO*K  (VEL*TIME)/RO>  *0PRA0 

THETA*7WT*BSI7E 

BHM/0PAA0 

TtFApTHETA/OFMO 

SEWJ,PULSE)«'nCTA«»2 

DENOH2.*SNff€U*AHP**2)*<2.843*l  .77iCO$(3.74*THETA/0U) 
»-1.0W*COS(4.188*TJCTA/8W)) 

cm*<H«2)/oeot 

sERC?<j.pasf)«aa 

THETF><ASIN(  (SIN(OU^)«SIN(AOSrnCTA) ) )  H.S) )  lOFIM) 

BMWOPRAO 

THETA*T1CTA«0PRA0 

SATA  PRE/UM.O/ 

ANGLE^l.+ABSUNTUD.iTHETP)) 

DO  «  WANGLE 
NCALC»<FIQAT<AN6LE)/100.)P< 
IA*5./2t2.*COS(2.0944CALC/lM)*l.5*COS(4.188<MCALC/SH) 
•*£0S(3.7MCALC/BU) 

RM(WI.3*2.*C0S(2.l944ttALC/BU)t|.3*C0S(4'lM4CALC/BI) 


CLOCK  FOR  ENTIRE  RtM 
M0U1N6  T6T  RANGE  UPDATE 


DIRECT  ADA  IS  T6T  LOCAT  i 
SPEC,  AQA  IS  U<6,BU/2)  I 
DIFF,  AQA  IS  U<0,PI) 

CHIT  ALL  REF  OUTSIDE  ftl/2 

i 

EFT  IS  RUCTION  OF  COEFF  I 


SW  CWNCL  OBSERVATION 
diff  cmta  oesEMATim 


DTINE  IS  PULSE  TIME  DIFF  < 

i 

WE  B08ES1SHT  EACH  PULSE1 
LOOP  0®  FOR  REFLECTIONS 

i 

UPDATE  HOVE  TGT  LOCATION  i 
COTUTE  THETA  AFTER  PRI 


TAACK1N6  ERROR 
TRACKING  ERROR  AFTER  PRI 
COMPUTE  CRMER-RAO  BOUND 


CRO  AFTER  PRI 
COMPUTE  2IU-ZAKAI  801N0 


•+C0S!3.7MCALC/8U)/U 
AX*(S8RT!(1.-RH0)*HM»€U/2.)>*AMP 
CALL  HONOR!  AX,  PX) 

(K»l.-fX 

NCALWCALC/UWAD 
PRE!N)«®!»(  SIN!  1 .  NCALC) )  *<2 
90  CONTINUE 
00  94  JK*1,13 
JX1«JK*1 
01  92  KJsJKl.H 
1F<PRE(JK).LE.PRE(KJ))60  TO  92 
STORE-PRE(JK) 

P«E<JX)»PIE<U) 

PREdUWTORE 
92  QHTINUE 
94  CWTINUE 
22WRE!1« 

S£822(  J, PULSE  W2I 

1M  CONTINUE 
2H  CONTINUE 

WRITE!*, 1930)MI 
WRITE!*, 1941) 

DO  411  K-1,3* 

SUM.) 

SUfl-9.8 
SUOI.I 
D0  390  L-1,13 
SUMUHSERiL.K) 

SU11-StNl*SERCR!L,K) 

SUO«$UG4SER22!l,K) 

391  CONTINUE 
AU6<K)*SUV13. 

AW61!K)-SUM1/13. 

AV62(K)»SIK2/13. 

URITE<*(1030>AM6<K)  t^M61<X> 

499  CONTINUE 
399  CONTINUE 
60  TO  3 

1999  FORMAT!  IX, 'NOVINO  TBT.STAT  T6T,  OR  QUITU:NOVE,2:STAT,0:QU!T>'> 
1919  F09NA7(lX,'SNW/33,39,23,20,13,J9,7,3t3fl,-2,-5,-7,'19,-20/V> 
1929  FORMAT!  IX, 'ENTER  MW  START/STOP  POINTSdE  START-1, STOP-13)') 
1023  FORMATOX, 'INPUT  RATE  OF  OCCURANCES/PRI  (IAN80A)') 

1025  FORMAT!  IX, 'U4180A  MUST  BE  RSH.  (9-11.9)',/) 

1027  FORMAT!/, IX, 'SPECUWR  OR  OIFFUSE  REF?(l:SPEC,2s0IFF)V> 

1939  FORMAT!/, IX, 'Ml  LEVEL  (08)  •  ',E11.3,/) 

1949  F0RMAT(14X,'4M6SEt',8(,'AV6SERCR',4X,'AU6SERZ2') 

1939  F0RMAT(19X,3!3X,E11.3)) 

3999  STOP 
M 


22B  AFTER  ffl 


AUEBA6E  ERRORS  ODER  15  RWSi 


OUTPUT  RESULTS 


FOfrAT  STATEHENTSJN/QUT 


Appendix  F 


Program  Listinc 


*««•»«*«»«#*«•****«*«******»******«*»**•»•»**»*******«««••**•**»**»*** 


PROGRAM: 

PURPOSE: 

INPUT: 

OUTPUT: 

COMPILER: 


•LOCKOUT’  OF  UNDESIRED  PULSES  WITHIN  PRKREQ  4  PULSES) 
COMPUTE  THE  TRACKING  ERROR  AND  PERFORMANCE  BOUNDS 
DESIRED  SNR  LEVEL (dB) ,  TARGET  MOVING  OR  STATIONARY, 
AVG  OCCURANCES  IN  THE  PRI ,  SPECULAR  OR  DIFFUSE  MULTI' 
3<1  BY  50)  VECTORS  CONTAINING  AVG  ERROR.  CRB,  2ZB 
FORTRAN,  VERSION  5 


EXT  CALLS:  IMSL  ROUTINES  GGWL,  MDNQR,  GGPOS,  GGEXN,  GGUBS 


PR06MN  NULFIX 

DW06IDN  AU6(50)  ,AV61  (54)  ,AU62<50)  ,PRE( 14)  ,11(1250)  ,8(1250) 

0IMBIS1«SNHIK15),SER<15,50),S£RCS<I5,50),SERZZU5,5« 

DIMENSION  11(90)  ,E<24)  ,0TIME<24)  ,U(23)  ,AN6(2S)  ,AWAL<4,2A> 

REAL  NPWR.N1  ,N2,NAMP,NCALC,PULSEC 

INTEGER  PULSE /N6LE, REF, FLAG, FLA61,FLA£2,FLA63 

DOUBLE  PRECISION  DSEED1  ,DSEED2,D$EED3,DSEED4,DSEED5 

DATA  SMN/33.,30.,25.,20.,15.,10.,7.,5.,3.,1., 
«-2.,-5.,-7.,*!i.,-20./ 

DMTDUE 

URITE(*,1000) 

SEAO»,IMS 

IF(D1S.EQ.O)60  TO  9001 
URITE<»,1010) 

WUTE(»,1O20) 

READ*,  II,  1 2 

DATA  DSEE01  .DSEE02/223497.0M, 329017.000/ 

DATA  DSEED3,0SEED4 .0SEED5/123457 .000 , 1 23497.000 , 1 23497 . 000/ 
DATA  DPRAO, R0,UEL,FV, (M/37,3,4900,, 297. 4, l.E-6, 3.0/ 

DATA  PR1 ,SI6, OMEGA  ,AMP,Sftft/10,E’3, 434., 434., 3.0,1 .£*3/ 

NR-1250 

N0-123# 


FLAM 

00  300  IM1,I2 

SNft*SMWlP> 

9NMB4>1O.**(9NR/10.) 

ROUHI0 

WRITE!  *,1023) 

URITE(*,1023) 

REAiX.RU* 

WRITE! • i 1027) 

READ*, REF 

1F(RU*.LE. 0.0)60  TO  3 
IF<RU*.6T. 10.0)60  TO  3 
DO  200  J-1,13 
DATA  ARUAL/1Q4*!,0/ 


DIM  ARRIVAL  TIMES  4  PRl'S) 
THREE  ADDITIONAL  FLAGS  I 


ONLY  WAY  OUT  IS  0  IN  IMS 


INITIALIZE  AU.  A8LMLM.0 
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INIT  FIA63  (LOCK  INDICAT) 


FIA6>I 

CALL  66NKDSEED1  ,NR,R) 

CALL  G9H(OSEE02,N8,5) 

Ml  S6MS(RtA1,DS£ED3^P,IR,]ER) 

BSITM.I 

TR6M.8I 

1F(WS.EB.1)BSJTE*1.5 

IF<INS.EQ.1)T*ST«4.I 

TJ€TA»T1ST-BSITE 

M 

00  111  PULSE-1 ,31 
LAST— IR<PULSE)-*1 
NEX-IKPU.SE) 

Rt»<l.lE-2)/R0AT(LAST) 

CALL  66E»K0SEED4f»1,NEXf£) 

MHCX 

CALL  6aJBS(0SEED3^U,U) 

IF(NEX.LE.1)60  TO  43 
DO  M  LL-I  /CX-l 
IPI-LL*! 

N  93  KJMP1,NEX 
IF<E(LL).LE.E<KK))GO  TO  33 
TB^€(LU 
E<LL>— €(KK) 

EUKMBT 
33  CONTINUE 
41  CONTINUE 

43  01  83  NULTI>I(NEX 
IF<FIA6.NE.1)60  TO  67 
FLAM 
60  TO  83 
67  CONTINUE 

IF(NULTI.EO.O)E<NULT1)*I.I 

E(LAST)*PfU 

lF<E(HULTM).6E.Pil)60  TO  83 
PTIMEtNULTl+I  )-€(«ULTI*l  K<NULT1) 
IF<FlA63.EQ.l)DTINE<HULTI«t>-PRS 
C0EFW.4 
DO  71  II-I.NU 
ANG<IIMKII)*4.283 
71  CINTIHUE 

IF(NULTI.EO.O)COEFF«I.I 

1F<0TINE<HULTI*1>.8T.2*PU)60  TO  73 

VNI-COEFF 

ANSI— NN6CNULT1 ) 

VW-8.4 

AN62*N8<NULTIM) 

UNT*S0RT((UNHC0S<AN61)«VH2iC0S<AN62)>H2 
«<VN1  »8IN(AN61  )HH2»SIN(AN62) )  «»2) 


PUCE  TAB6ET  IN  8EANUIDTH 


HACK1N6  CALCULATIONS 


HULTIBATW  PW  DO  LOOP 


* 

ON  LOCK  OBS  PEXIOO  IS  P*I« 


•Au1!! 


l  FOLSE  OVERLAP 
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COEFF-UHT 

FLAW 

73  CWT1NUE 

PULSEC*FLQAT (PULSE) 

TINWRIKPULSEC-l  .)*£«JLTI*1) 

»BM«T<Mh24(VEUTIIC)  H2) 
lF<IMS.E0.1.)SNa€)K(R0LD/lt£U)<*4)<SNf#€U 

Ifly^Muft/CyfluQi 

IfUraflw  9HVM 

NAMNS*T<NPUR) 

I«I*I 

NMKIMWT 

NM<!)«MHP 

IFOWLTI.Ea.D60  TO  81 

1F<«EF.E0.1)TKTAMKMULTI)«1.S 

1F(REF.E0.2)T10HI<MULTI)»18I. 

IFOICTA.ST.  1.5)60  TO  85 
88  CONTINUE 

!F(PULSE.l£.4)ARtJAL(PULS£,MULTI«l)*E  (MULTI) 

1F(PULSEJC.S)60  TO  84 
00  82  JJ>1,23 

T0TALa<ANAL(llJJJ)+ARUAL(2|JJJ)+ARUAL<3,JJJ)tAMAL(4lJJJ))/4.B 
lF(T0TALJ€.AfNAL(l,JJJ>)SO  TO  82 
lF(T0tAL.NE.AMAL(2,JJJ))6Q  TO  82 
IF(TUTAL.NE.AWAL<3,JJJ))60  TO  82 
1F(T0TALJ«^W8L(4I4JJ))60  TO  82 
IFCTOTAL. £0.1.0)80  TO  82 
IF(TOTAL.EQ.PR1)60  TO  82 
F1A62-AM-S 
FIA6>1 
82  CKT1NUE 

84  CWTINUE 

E6Y-FW*<C0EFF<AnP)H2 

V0-S1N<1.88«THETA/8U) 

V>(S0*T(2.)/2.)t<l.»C0S(2.8?4«T»€TA/8y)) 

Rl*C0EFF«AMP»$0RT(E6Y)*VS<Nl 

R2-C0EFF*MPtS0RT<E6r)«U0«N2 

ERD-RlMWttiVS 

IF<ERR/(SQRT(2.*E6Y)<C0EFF*AMP)  .8E.0.99) 
*ER1N8.9?»SQRT<2.*E6Y)*C0£FF«PMP 
!F(ERt/(S0RT(2.*£6Y)*C0£FF«AMP).LE.H).??) 
•£8N,"8.9?*S0RT(2.»£6Y)»C0EFF*4NP 
DTJCTAMSU/1  .88)*ASIN<ERR/<S8RT(2.*E6Y)*C0EFF«PMP)> 
S£JM^I.O-£XP(-SI6«OT!hE<NU.TlM))»<COS(ON£»»OTlN£(NULTM)) 
t-$iN<QME6A*0TIME<NULTI*l))) 
lFtVO»ERR.ST.0.Q)DTHETAM)'n€TA 
1F(NULTI  .EQ.FIA62)8SITE*8S1T£*S£IMWTHETA 

85  CONTINUE 


PULSE  OUERIAP  MONITOR 
00  PULSE  OVERLAP  COMPUTE) 


STORE  ARUM.  TIMES  FIRST  4) 
LOOK  FOR  PERI001C  ON  5TH 
QMSIDEh  23  PULSES  POSSIBi 
AU6  25  PULSES  OUR  4  PHI'S 
DOES  AU9-PRI1,  PULSE  JJJ?  i 
DOES  AUWR12, PULSE  JJJ?  i 
DOES  AUWRI3, PULSE  JJJ?  I 
DOES  AU04H14, PULSE  JJJ?  i 
RECALL  ANAL  U4IT  TO  1.0  i 
NEXT  PR1  NOT  DESIRED 
JJJ-1  BECAUSE  DIRECT  IS  0< 
SET  k  LOCKED  TO  PULSE  FL2< 
LOOP  THRU  ALL  25  PULSES  < 
OUT  IF  ROU.COL  .NE.  AU6  i 
B#  PULSE  LOCK  COMPUTE  < 


LOCK  ON  FIA62,  OUT  OTHERi 
00  MULTIPATH  DO  LOOP  i 


lF(lMS.E0.1)TR8TaATAN<<VEL«TlME)/RI)H>PHA0 

THETA-TRST-BS1TE 


artu/DPMo 

■nCT^THETA/DmO 


SRTJ.PULSEWHETAmZ 

DENOTO<2.*SM8€WN81Pm2)«(2. 845+1. 77*C0S(3.74»T1€TA/Bd) 
•-1  .895»Ct5S<  4 .188+THETA/8U) ) 
d^<BU«2)/D0O1 
SERCR<J,PULSE)«CR8 

T)CTI*<ASIN<  ( S1N(8M/2)*SW(A8S(T1CTA) »«8 .5) )  +OPRAO 

nwwiDPwm 

THETA*THETA«OWAD 

DATA  m/UH.V 

ANGLE“1.*ABS(INT(18.»THETP)) 

00  98  101  .ANGIE 
NCAUMFlQATIANGtfVIGG.)* 

IO5./2+2.»Cfl$<2.894NCAlC/8U)+l.5«C0S<4.1 888CALC/W) 
*-C0S(3.74<NCAlC/BU) 

BM^<«.5+2.*C0S<2.8?4*CALC/BU)+I.5*C0S(4.188*<ALC/BU) 

•+C8S(3.744CALC/BU)/V 

AX*<S0RT(<1  .-fiHO)N+*9*NBl/2.))*AHP 

CALL  rCMM (AX, PX) 

CDM.-PX 

NCALMCALC/OFMO 
«£«)»«•<  S1HU.ACALO)  m2 
98  CONTINUE 
00  94  JK»1,15 
JM-JK+1 
D*  92  X>4K1,14 

JF<PRE(JK).lE.PRE<KJ»G0T0  92 
STORE«PRE<JK> 

PRE<JX»PRE<KJ> 

PRE(KJ)*ST0RE 
92  CWTINUE 
94  COfflNUE 
ZZWKE(U) 

SERZZU.PULSEWZB 


TRACKING  ERROR 
CtHftTTE  CMHER-SAO  BOUND 


Z1V-ZAXAI  BOUND 


181  OHT1NUE 
208  CONTINUE 

URITE(«,1830)M 

URITE(*,1040) 

DO  480  10*1,38 
SUOO.O 
SW1«8.0 
SltOO.I 
DO  388  1*1,13 
SlfOSW+SERtt.X) 
$Wl*SIM+SERCRa,K> 
SUQ*$W2+SERZZa,K> 
388  CONTINUE 


AUERASE' ERRORS  tVER  15  RWS< 


AV6<K)*S0V15.  > 

AV6KK>*Stl11/15.  » 

AW62(K>=SIH2/15.  • 

*  OUTPUT  RESULTS 

URITE(«,185C)Ayfi(K)^WGl(K)>WS2(K)  • 

400  CONTINUE  • 

SOI  C04TINUE  i 

60  TO  5  * 

*  FOUST  STATEMENTS  IN/OUT 

1000  F0WAT(1X,'N0V1N6  TCT.STAT  TE7,  OR  QUITO  M)UCt2«STATt8:QUIT)'>  * 

ill)  F0»*T(lX,'SNW/35,30,25,2l,15,ll,7,5,3,l,-2,-5,-7,-ll,-20/V>  * 

1021  F0RST(  IX, 'ENTER  SNW  STMT/STOP  POINTSdE  S7AHT-1 , ST0P»1 5) ' )  * 

1023  FORMAT (IX, 'INPUT  RATE  OF  OCOJWCES/PHI  (UOfiOA)')  * 

1025  F0*AT(1X, 'LAMBDA  MUST  BE  REAL  (0-10.0)'/)  • 

1027  F0ft(AT(/, IX, 'SPECULAR  OR  DIFFUSE  REF?U:SPEC,2:0IFF)V)  * 

1031  F0MAT(/,1X,'M  LEVEL  (OB)  >  '.E11.5/)  • 

114  F0fmT(l«,'AV6SER',8X,'AV6SERCX',A,'AV6SER22')  • 

1058  F0*AT(1BX,3(»,EU.3)>  • 

5000  STOP  • 

END  • 


Appendix  G 


Program  Listing 


************************************************************************ 


*  PROGRAM:  TWO-DIMENSIONAL  AC  MONOPULSE  RADAR  SIMULATION  # 

»  PURPOSE:  COMPUTE  THE  TRACKING  ERROR  AND  PERFORMANCE  BOUND  » 

*  INPUT:  DESIRED  SNR  LEVEL(dB)  * 

*  OUTPUT:  2(1  BY  50)  VECTORS  CONTAINING  AVG  ERROR  St  ZZB  » 

*  COMPILER:  FORTRAN,  VERSION  5  * 

*  EXT  CALLS:  IMSL  ROUTINES  GGNML  St  MDNOR  » 


***ftft***ft#*#**#****#***ft*#***#*ft*ft**«»*##*«*****#**#*ft*#***«***ft«ft*ft#*#* 


PROGRAM  TUOOIM 

DIMB4SI0N  AU6<50) ,AW62(59) ,PRE(14,16) ,R(50) ,8(58) ,C(50) 
DINB6IW  SOW<13),SER<15,50),SEB2Z(15,50> 

REAL  NFW,N1IN2,N3,NAW,NMLC1,NMLC2 
INTEGER  PtR  IE  ,ANGL£1  ,AN6LE2 
DOUBLE  PRECISION  DSEED1 .0SEED2.DSEED3 
DATA  SNRV/35.,30.,25.,20..15.,10.,7.,5.,3Ml.f 
e-2.,-5.,-7.,-10.,-20./ 

3  CCNTINUE 
URITE(*,100I) 

READ*,  INS 

IFUNS.EQ.DSO  TO  3001 
URITEO.ltll) 

URITE(»,H20) 

RBtte.il.I2 

DATA  DS£EDl,OSEEO2,DSEED3/123457.O0fl, 325017.000, 133547.000/ 
DATA  DMAO,Pg,8U/57.3,l.E-4,3.0/ 

DATA  PRI  ,S16, (MEGA, AMP, SAA/10.E'3, 634. ,636. ,5.0,1  .£-3/ 
E6Y*W*AMPh2 

NR*  50 
ND*50 
NC*SI 

DO  300  IP*!!, 12 
SN«*SNRV<IP) 

SMNB6*II.«t(9W/IO.) 

DO  200  3*1,13 

CALL  GMKOSEEOIMR)  ' 

CALL  G9tl(0SEED2,N8,B) 

CALL  Ga01(DSEE03,NC,C) 

8$ITEL*0.0 

BSITEAM.O 

TRGTEL*0.80 

TR6TA2*I.4I 

THETA*TR6TEL-BSITEl 

PH1*TR6TAZ*8SITEA 

1*0 

00  100  PULSE-1,31 


DECLARE  VARIABLES 


RADAR  MttAMETEBS 


PLACE  TARGET  IN  gEANMIDTH 


TRACKING  CALCULATIONS 
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PULSEC=n.OAT  ( PULSE) 

NPUR=SWR/SN«NEU 

NAMP*SQRT(NPUR) 

1*1*1 

Ni=Rd)#m> 

N2*B(I)#W1P 

N3*C(1)4M«> 

VDEl*THETA/BU 

VDA2=PHI/0U 

vs»i. 

R1*AMP»SQ8T ( ESY) *VS*N1 

R2»AMP*SCRT<E6Y)*V0£UN2 

R>ANP*S9RT(E6Y)<V0A2*N3 

ERREL*R1  *BM*THETA*R3*THETA*PH  I  -R2«(  0Uw2*PHI  *«2) 

ERRA2=81  *8U*PH1  *R2*T1€TA*PHI“R3*<  BUm2*THETA**2) 
DTHETA-ERRELy<AHf>»SffiT(E6Y)*BU) 

DPHI»-ERRAZ/(flNP*$ORT  ( E6Y)  *8U) 

SERVE!*! . 0-£XP( -SI 6*PRI ) •( COS( 0HEGA*PRI ) -SIN( 0MEGA*PR1 ) ) 

SERMAZ*1 . 0-EXP(-Sl S*PR1)*( COS( 01EGA*PRI ) -SIN< DMEGA*PRI ) )  ' 

1F<VDEL*ERREL.CT.O.O)DTH£TA«-OTHETA 

IF(VDAZ*ERRAZ.6T.8.8)DPH1*-0PH1 

BSlTtt*BSlTEL*SEf?;a«OTKETA 

BSITEA*8SITEA+SERVAZ*0PH1 

THETA*TRSTEL-BSITEL 

PHI*TRfiTAZ-8SITEA 

BMU/DPRAO 

THETA*THETA/DPRAD 

PHI*PHI/DPRAD 

SER(  J  ,PULSE)*OI€TA«2*PHl  *2) 

THETP*(ASIN((SIN(BU/2)*SIN<A8S(TH£TA) )  )*0 .5) )  tOPRAO 

PHI P*(ASIN( ( S1N< BU/2)*SIH<  A8S( PHI ) ) ) *8 .3) ) *OPRAO 

BU*8U<DPRA0 

THETA»THETA*DPRAO 

PHI*PHi»OPRAO 

DATA  PRE/254»0.0/ 

ANSLEW  .*A8S<INT(1D.*THETP)) 

AN6LE2*1 .  *A8S<  INTU8 .  »PHIP) ) 

DO  98  N*1  .AN6LE1 
DO  83  M*1,ANGIE2 
NCALC1*<FIQAT(AN$IE1)/108.)<N 
NCALC2*(FIMT(AN6LE2)/100 . 

UI*5./2*2 .  *COS<  2.0W«NCALC1/BU)*0 .5*C0S<  4 . 1 88WCAIC1/8U) 
*-COS(3.74«CALCl/BU) 

U2*5./2*2.*COS<2.0P4<4CAIC2/BU)*0.5*COS<4.188«NCALC2/BU) 

*-COS(3.76«CALC2/0U) 

RHOl*<  0 .5*2. «COS<  2 .0?4*NCALC1/BU) *0 ,3*C0S( 4 . 1 88«NCALC1/GU> 
**C0S( 3 . 74*HCALC1/BU)/U1 

RH02»(  0 .3*2.*COS(2.094<MCALC2/8U)*8  .S«COS(  4.1 884KALC2/8U) 
**COS<  3 . 74«NCALC2/8U)/U2 
AX1*(SQRT<(1.-RH01)<U1*(SNRNEU/2.)))<ANP 
AX2*<S0RT((l.-RH02)AC»(SNWa/2.)))«AMP 
AX31*(AHP4CALC1/BU)  «SQRT  ( 2 .  *SNIMEU) 


a 


c 


>■. 

L‘. 

L*. 


ST*SWT  ( (NCALC1/DPRAD)  »*2*<NMLC2/DPfiAD)  «2) 
AX32«<flMP*ST*OWW)/8U)*SORT(2.*SNSN&) 
AX33=<AHP*NCALC2/BU)  nSORT(  2 .  *SNJWEW) 

MLL  N0NQR(AX1,PX1) 

MLL  NONORwv2fPX2) 

MU  MDN0R(AX31  ,PX31) 

MLL  raN0R(AX32(PX32) 

MLL  MCN0R(AX33,PX33> 

0X1*1. -PX1 

ox>i.-px2 

QX31*1.-PX31 
(K32*1.-PX32 
8X33*1. -PX33 
NMLC1-NMLC1/DPM0 
NMLC2-NMLC2/DPRAD 
PEST*.333*<8X31 *8X32*8X33) 
PRE<N,r*«*.125»<SlN(ST)H2)*PEST*<2^<SlN(ST)«#2)) 
*•<0X1  *SIN(  1 .  4CALC1  )<*2)  •<0X2>S1N<  1 .  JNMLC2)  **2) 
NMLCl^MLCUOPMD 
NMLC2»NMLC2«DPRA0 
80  CtNTJNUE 
90  CONTINUE 
DO  96  M*l,15 
DO  94  JK*1,15 
JK1«JK*1 
DO  92  KJ*JK1,14 

IF<PRE<JJ,JK).LE.PRE<JJ,KJ))60  TO  92 
ST0RE1*PRE(JJ,JK) 

PRE(JJ,JK)*PRE<JJ,KJ) 

PRE<JJ,KJ)«STORE! 

92  CONTINUE 
94  CONTINUE 
94  CONTINUE 
DO  98  JL*1|13 
JL1*JL*1 
DO  97  LJaJll.14 

IF(PRE(JL,14).LE.PRE(LJ,14))60  TO  97 
ST0RE2*PRE(JL,14) 

PRE<Jl,14)*PRE<LJ,14) 

PRE(LJ,14)=ST0RE2 

97  CONTINUE 

98  CONTINUE 
ZZWRE<14,14) 

SERZZ(J,PULS£)*ZZB 

100  CONTINUE 
200  CONTINUE 

URITEO, 1030)901 
UR1TE«,1040) 

DO  400  K>!,58 
SUH.O 
SJ12*Q.fl 
00  300  1*1,13 


ZIV-ZAKAI  BOOM) 


AUEMBE  ERRORS  WED  15  R1HS) 


156 


Stft*SUHS£R<L,K> 

SUM2=Sm2+S€R2Z(L,K) 

308  CONTINUE 
AM6(K)»SUN/13. 

AW62(K)«SUN2/13. 

WITE(«,1050)AVG<IC),AVG2<I<) 

430  C04TINUE 
300  CONTINUE 
SOTO  3 

1000  FOHWT(lXf'QUIT?<ENTER  1  TO  PSOCEED,  8  TO  QUIT)') 

ioio  F0iitATdx/SNW/3Sf30, 23,20,13,10, 7, 5,3,1,-2,-s, -7, -io,-28/V> 

1820  FOItOTdX, 'ENTER  SWN  START/STOP  POINTS! IE  START-1, STOP-15)') 
1030  F0SHAT</,1X,'SNR  LEVEL  (OS)  >  ',£11.5/) 

1040  F0Wffd4X,'AUGSER',2SX,'AW6SERZ2') 

1030  F0RNATdflX,3X,E11.5,17X,E11.5) 

3000  STOP 
END 
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monopulse  radar,  generates  an  error  signal,  and  excites  a  servo  system 
to  reposition  the  boresight  on  target,  thus  providing  target  tracking. 
Cramer- Rao  (CR)  and  Ziv-Zakai  (ZZ)  bounds  are  derived  to  indicate  system 
performance  under  varying  signal-to-noise  ratio  (SNR)  conditions. 

Actual  tracking  error  is  approximated  from  a  computer  simulation  and 
compared  to  the  CR  and  ZZ  bounds  for  varying  SNR  levels,  for  tracking  in 
•slow"  and  "rapid*  target  fluctuation  environments,  and  for  tracking  in 
the  presence  of  specular  and  diffuse  multipath.  A  two-dimensional 
tracking  model  and  associated  performance  bound  are  aiso  presented.^ 

At  high  SNR  levels,  the  CR  bound  results  lower  bounded  the  mean- 
square  tracking  error  (MSE),  but  for  low  SNR  conditions  the  CR  bound 
results  exceeded  the  MSE.  The  ZZ  bound  results  indicate  a  tight  lower 
bound  for  the  MSE  at  low  SNR  levels  and  in  both  the  multipath  and  target 
fluctuation  environments .  '“Slow"  target  fluctuations  and  diffuse 
multipath  results  indicate  that  target  tracking  capability  is  not 
seriously  degraded  in  either  of  these  two  surroundings. .  Conversely, 

P Rapid*  target  fluctuations  and  specular  multipath  environment  results 
indicate  serious  tracking  degradation  is  introduced  in  tlie  monopulse 
tracker.  I 
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